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v. Borböly in Berlin. 


On non-stationary spatial low around a wing. The problem of plane flow around an 
(infinite) oscillating wing has already been solved. The author treats first with sufficient 
completeness the (spatial) case of a finite supporting line with an elliptical vortex distri- 
bution. Then he shows that these results may be combined in order to study the finite 
oseillating wing with elliptical vortex distribution whose depth is small in comparison with 
its span, and obtains approximate formulae for the flow in the neighbourhood of the wing 
and for the angle of attack. 


Sur l’6coulement spatial non-stationnaire autour d’une aile. Le probl&me de l’&coulement 
plan autour d’une aile (infinie) oseillante est d6j& r&solu. L’auteur traite d’abord dans 
des conditions suffisamment larges le cas (spatial) d’une ligne -portante finie entourde de 
tourbillons & distribution elliptique. Ensuite il sait combiner ces deux re&sultats pour 
ötudier l’aile oscillante finie & distribution elliptique des tourbillons, sous la condition que 
sa profondeur soit petite par rapport & son envergure. Il deduit des formules approch6es 
pour les vitesses au vc!sinage de l’aile et pour l’angle d’attaque. 


A. Kneschke in Dresden. 


On the motion of vortices. The paper deals with the motion of isolated vortices in any 
simply connected plane domain. The author establishes a relation between the kinetic 
‚energy and Kirchhoff’s „Bahnfunktion“ and deduces the formula for its transformation 
corresponding to conformal mapping of the domain. The motion may also be derived from 
a variational principle. 


Sur le mouvement des tourbillons. L’auteur &tudie le mouvement de tourbillons isol6&s 
dans un domaine plan simplement connexe quelconque, en calcule l’&nergie cinetique qu’il 
met en relation avec la „Bahnfunktion“ de Kirchhoff, et donne la formule pour la trans- 
formation de cette önergie correspondant & une repr&sentation conforme du domaine. Le 
mouvement peut aussi @tre deduit d’un principe de minimum. 


E. Reuß in Budapest. 


On Lüders-Hartmann lines. The paper is concerned with the formation of gliding surfaces 
in an elastic body at the yield point, as shown experimentally by the Lüders-Hartmann 
lines. The author considers, as a typical case, the twisting of a cylinder of any cross- 
section. At the yield point the material will still follow Hooke’s law nearly everywhere, 
but there must be some surfaces (gliding surfaces) where it behaves like a plastic body. 
Castigliano’s principle shows these surfaces to be portions of planes parallel to the axis 
of the cylinder. The strain and stress, and in particular their discontinuities on the gliding 
surfaces, may be calculated by means of the logarithmic potential of a simple distribution 
on the gliding surfaces, whereof the density is determined by an integral equation of the 
first kind. E3 


Sur les lignes de Lüders-Hartmann. Il s’agit de la formation des surfaces de glissement 
dans un corps ölastique & la limite d’&coulement, surfaces mises en &vidence experi- 
mentalement par les lignes de Lüders-Hartmann. L’auteur considöre comme cas typique la 
torsion d’un cylindre de section quelconque, A la limite d’&coulement le materiel ob6it 
encore presque partout & la loi de Hooke, mais il doit y avoir quelques surfaces (surfaces 
de glissement) oü il se comporte comme un corps plastique. En vertu du principe de 
Castigliano ces surfaces sont des portions de plans parallöles ä l’axe du cylindre. Pour 
calculer les d&formations et les efforts, et en particulier leurs discontinuit6s aux surfaces 
de glissement, l’auteur se sert du potentiel logarithmique d’une simple couche ötendue sur 
les surfaces de glissement. La densit& en est determinde par une &quation intögrale de 
premiere espece. 4 
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HAUPTAUFSÄTZE 


Über einen Grenzfall der instationären räumlichen 
Tragflügelströmung*). 
Aus dem Flugtechnischen Institut der T. H. Berlin. Leiter: Prof. Dr.-Ing. Herbert Wagner. 


Von ®. Borbely in Berlin. 


ie Grundaufgabe der Tragflügeltheorie, zu einem Flügel mit endlicher Spannweite die 
D zugehörigen Luftkräfte zu berechnen, ıst schon beim stationären Problem nicht einfach. 
Die strenge Lösung der entsprechenden Aufgabe für den nicht-stationären Fall erscheint noch 
wesentlich schwieriger. 


Den prinzipiellen Ansatz könnte man in Anlehnung an die Idee der tragenden Wirbel- 
fläche, die ım stationären räumlichen und ebenen instationären Fall durchgeführt worden ıst')?), 
folgendermaßen formulieren : 


Der Flügel wird als unendlich dünn vorausgesetzt und als Trennungsfläche (Diskontinuitäts- 
läche der Geschwindigkeit) aufgefaßt, also durch flächenhaft verteilte Wirbel ersetzt. Die 
/Zirkulationsdichte dieser Trennungsschicht ist zeitlich und örtlich variabel und fürs erste 
unbekannt. Der Gesamtänderung der Zirkulationsdichte der gebundenen Wirbel entsprechend 
entsteht nach dem Helmholtzschen Satz eine weitere, aus freien Wirbeln bestehende 
Trennungstläche, die von der Strömung mitgeführt wird. Für hinreichend kleinen Auftrieb 
ist die Deformation der abgehenden Wirbellläche zu vernachlässigen. Man kann dann nach 
dem Bıi1ot-Savartschen Gesetz aus der Zirkulationsdichte des gesamten Wirbelsystenis 
(las zugehörige induzierte Geschwindigkeitsfeld in allgemeiner Form angeben. Die Zirkulations- 
dichte der gebundenen Wirbelfläche ist nun so zu bestimmen, daß das resultierende Gesechwindig- 
keitsfeld an der Oberfläche des bewegten Flügels die geometrische Grenzbedingung erfüllt. 
Damit sind auch die am Flügel wirkenden Luftkräfte ohne weiteres bekannt. 


Die Lösung der obigen Fragestellung bietet jedoch fürs erste zu große Schwierigkeiten. 
Darum versuchen wir im folgenden den Einfluß der endlichen Spannweite auch beim nicht- 
stationären Fall wenigstens in erster Näherung analog dem Prandtlschen Gedanken durch 
eine Korrekturgröße zu erfassen. 

*) Auszug aus der von der Technischen Hochschule Berlin angenommenen Dr.-Ing.-Dissertation des Verfassers. 
Referenten : Prof. Dr. R. Rothe und Prof. Dr. G. Hamel. 

Herrn Prof. Rothe und seinem ehemaligen ständigen Assistenten Herrn Prof. Dr.-Ing. G. Grüß möchte ieh 


auch an dieser Stelle für die vielen Förderungen meiner Berliner Jahre aufriehtig danken. 


!) Birnbaum: Die tragende Wirbelfläche als Hilfsmittel zur Behandlung des ebenen Problems der Tragflügel 
theorie. ZAMM Bd. 3 (1923), S. 290. 


®) Birnbaum: Das ebene Problem des sehlagenden Flügels. ZAMM Bad. 4 (1924), S. 277. 
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320 v. Borbely, Grenzfall der instationären räumlichen Tragflügelströmung Math. 
In der x y-Ebane des »yz-Raumes liege im Bereich z=a...a+t, y b...+b der 

Die Zirkulationsdichte der entsprechenden 

Für 


KFlürel von der Tiefe f und Spannweite 2h(Abh. I 
eebundenen Wirbellläche sei zeitlich (r) mit der Frequenz o harmonisch veränderlich. 


jeden Zeitpunkt der Sehwingung seı in Richtung der Spannweite die Zirkulationsdichte will- 


kürlich als elliptisch vorausgesetzt. Die Zirkulationsverteilung in Richtung der Tiefe nehmen 





‘) 


Abb. 2 





Abb, 1. 

wir als reell an®), über ihre Größe y(z#) wollen wir dagegen keine Festsetzung treffen. 
V (parallel zur positiven «-Achse) sei die Anströmgeschwindigkeit der Flüssigkeit in großer 
Entfernung vom Flügel. Die Zirkulationsdiechte der gebundenen Wirbelfläche läßt sich also 


ın diesem Fall 
Il. (Flüeel. räumlich) in der Form 


I(#,y:rT) A J 
. m 


ar (‚r) e! 
‘ 


wir unter denselben Verhältnissen wie oben .hilfsweise folgende 


darstellen. Nun betrachten 
drei weitere Fälle der gebundenen Wirbelverteilune: 
2. (Faden, räumlieh). An der Stelle # = a + f* (f* = Schwerpunktsabszisse der y (x) Fläche) 

eehbundener Wirbelfaden (Abb. 2). dessen Zıirkulationsdiehte durch 


liege eın 


/ (y:rT) / A 1, 


sei so bestimmt, daß die Gesamtzirkulation des Falles 1 pro Längen- 
2 überein- 


eegeben sein möge. 1, : 
einheit der Spannweite mit der entsprechenden Gesamtzirkulation des Falles 


stimmen soll. 
a—+t liege eine in Richtung der Spannweite 


eben). Im Bereich z=«a... 
der a =T,Yyde’" und 


Wirbellläche mit Zırkulationsdichte 


3. (Flügel, 


unbegrenzte gebundene 
weiterhin 
eben) an der Stelle z=«a-+f* ein unendlich langer Wirbelfaden mit deı 

wobei wiederum die Gesamtzirkulationen pro Längeneinheit 


/ p ei® 


t. (Faden, 
Zirkulationsdichte /(r) 
in Richtung der Spannweite im Falle 3 und 4 übereinstimmen mögen, d. h. es seı 


+ f 
"\y@)de=I,IT,=T,. 


Aus den Zirkulationsdiehten der gebundenen Wirbel ermitteln wir im $ 1...3 für alle 
oben angeführten Fälle die Zirkulationsdichten der zugehörigen freien Wirbelllächen und aus 


der Gesamtzirkulation yermöge des Biot-Savartschen Gesetzes 


ware: 

öu, (le 
7 e|° 

Diese Geschwindiekeiten stehen zur 


die in der «= y-Ebene induzierten Geschwindigkeiten. 
r „Ebene senkrecht, wir geben daher ım folgenden nur ıhre Beträge an und nennen sie 


auch bei komplexer y/x) Verteilung. 


Intsprechendes gilt 
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Abwindgeschwindigkeit w. Die beigefügten Indizes entsprechen der Numerierung der ein- 
zelnen Fälle. 

Die Zirkulationsdiehten und somit auch die Abwindgeschwindigekeiten der beiden räum- 
liehen Fälle sind durch die obigen Festsetzungen eindeutie bestimmt, die der beiden ebenen 
Fälle enthalten dagegen noch einen willkürlichen reellen Faktor /,, über den wir nun folgen- 
dermaßen verfügen wollen: 

Wir betrachten längs der Geraden 2 =0, b<y=9,<-+b die Abwindgeschwindigkeit 
der beiden räumlichen Fälle 1, 2 und vergleichen sıe mıt den entsprechenden Geschwindig- 
keiten derjenigen beiden ebenen Fälle 3, 4, deren gebundene Zirkulationsdiehte pro Längen- 
einheit der Spannweite mit der gebundenen Zirkulationsdiehte der räumlichen Fälle ım be- 
r vb’ —y 

, b 
Festsetzung von /', ist nun längs jedes festen Schnittes y=y, (-b<y<-+b) zu der Ge- 
schwindigkeitsverteilung se, (#, 9; 7), mw, (x —- E*, y,;r) des räumlichen Flügel- und Fadenfalles 
die durch die Wahl von y, eindeutig bestimmte Geschwindigkeitsverteilung m, (#, yo; T), 


trachteten Schnitt y= y, übereinstimmt. Das heißt, es soll sein Durch diese 


w,(.e #*, 9:7) des vergleichsweise herangezogenen ebenen Flügel- und Fadenfalles zugeordnet. 
Im $ 4 wollen wir zeigen, daß in hinreichend großer, jedoch endlicher Entfernung x vor 
{ 
dem Flügel auch für endliches aY 
- U) 
| /e\l 
W,,3” W.,|1 I 0| — || 
L vr, 


i . a a 2 h \ er I)(E") 
eilt. Dabei wird wie üblich mit O (£*) eine Funktion bezeichnet, für die lim 
e>u €" 

Hinter dem Flügel ist dagegen eine solche Übereinstimmung auch in beliebig großer 
Entfernung nur in der Grenze kleiner Flügeltiefen im Verhältnis zur Spannweite, also mit 


konst.==0 gilt. 


‚> zu erzielen. 
Z)) 


Im $ 5 zeigen wir weiterhin, daß auch ın unmittelbarer Umgebung des Flügels (und 
für die Berechnung der Korrekturgröße ist dieser Bereich der ausschlagzebende) die Grenz- 
wertzleicehung 


” . (t ! 
lım tw, (2, 9; 7) — mw, (X, Y,5; Td)y we —E,y;5T)— wm, (2 —- 1, y;T7),]= lim 0 7 In Y V 
nu su) 
> ( > 
DD sb 
mit @a...2,1"...a-+t gilt, d. h. 
vw, (a,y;D)=%, (a, y,;tT) + lm fw, (x — FF, y;T) - w, (2 — 1", y,;Dd}. 
> 
2b 
Es ıst also für kleine ‘= — a in erster Näherung 
WET BE se et 7 


Wir beschränken uns daher im folgenden auf sehr kleine Flügeltiefen im Verhältnis zur 
h { 
Spannweite, d. h. auf den Grenzfall 5, >» 
Nun betrachten wir im ebenen Vergleichsfall 3 unter Berücksichtigung der Festsetzung von 
/‘, einen Flügelschnitt an der Stelle y, und versehen dieses Flügelelement mit der zur Anström- 
richtung senkrechten Zusatzgeschwindigkeit D(O,y,:r), oder (was auf dasselbe hinauskommt) 
D 


wir korrigieren den Anstellwinkel des betrachteten ebenen Flügels um das Zusatzgelied — , so 


\ 
erhalten wir bei varıablem ,, einen räumlich und zeitlich verwunden schwinzgenden Flügel, 
der in der Grenze für kleine Flügeltiefe mit dem schwingenden räumlichen Flügel gleich- 
wertig ist. 

Wir können also die bekannte Lösung des harmonisch bewegten ebenen Flürels*)’) als 
Ausgangswert nelımen und brauchen nur die Neigung des Flügelprofils im ebenen Vergleichhs- 
fall 3 (unter Berücksichtigung der Festsetzung von //,) an jeder Stelle der Spannweite um 
das rechtsstehende Korrekturglied zu vergrößern, um in erster Näherung den Einfluß der end- 
lichen Spannweite in diesem speziellen Fall der Wirbelverteilung zu erhalten, wobei wir zur 
Bestimmung des Korrekturgliedes nur die rechnerisch wesentlich einfacher zu behandelnden 
Fälle der instationären Einzelfäden zu betrachten brauchen 

4) Küßner: Schwingungen von Flugzeugflügeln. Luftf.-Forschg., Bd. 4 (1929). 


5) Kaßner u. Fingado: Das ebene Problem der Flügelschwingung. Luftf.-Forschg., Bd. 13 (1936), 8. 374. 
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I. Das induzierte Geschwindigkeitsfeld eines unendlich langen gebundenen Wirbelfadens 
mit zeitlich harmonisch veränderlicher Zirkulation. In der x=»y-Ebene des rechtwinkligen 
xy 2z-Raumes liege an der Stelle #»=a«a parallel zur „-Achse eın unendlicher gebundener 
Wirbelfaden (Abb. 3). Die den zebundenen Wirbel umgebende Flüssigkeit ströme mit der 
Geschwindigkeit V parallel in Richtung der positiven «w-Achse. Die Wirbelstärke des ge- 
bundenen Wirbelfadens sei /'(r)=T', ei”", Nach den Helmholtzschen Wirbelsätzen wird 








Abb. 5. 


sieh von dem zebundenen Wirbelfaden der (zeitlichen) Änderung seiner Wirbelstärke ent- 
sprechend ein System von entgegengesetzt eeriehteten freien Wiırbeln ablösen. welche als 
Trennunesschieht von der Strömung mitgeführt werden. Zur Zeit 7, hat der gebundene Wirbel 
die momentane Wirbelstärke /'(r,). Sehen wir in erster Näherung von der Deformation der 
abeehenden Trennungsscehicht ab, so wird in der Zeit 7, — r, die Änderune der zebundenen 
Wirbelstärke zur Zeit r, von der Strömung an die Stelle «, a— V(r, — r,) mitgeführt und 
bildet dort einen freien Wirbel von der Stärke 


Nach genügend langer Zeit wird sich also ın der Flüssigkeit folgendes Wirbelsystem befinden : 
\n der Stelle = «a der eebundene Wirbelfaden von der Stärke /'(r) [', e'°" und im Bereich 
a<x © die abeehende Wirbellläche mit der Wirbeldichte 


I.r » ! -. (dx. 


Durch das gesamte Wirbelsvstem wird an der Stelle P(£,») eine Abwärtsgeschwindigkeit 
induziert. deren Größe nach dem Bıot-Savartschen Gesetz 


H TI. TI. 
» » : n (x 
(di S)\ cd „ wi | | e Mn Add 71 
| TH / ‘ no 3 / ! ‚ E* ns f 93 
\ | lat =) (} N a | ’ R | Hr a > in NY 
H x 71 FH. 
oder 
I 
’ » 
ES - y i) f 
are | N Ti - da (1) 
Fi ee | ri ! = 
A 
- () () r j 
beträgt. Setzen wir --(r.- ID=u und -= tt, so ereibt sıch®) 


‚ I 
I ac N‘ ! inet s \ MET 
a Pi £ \ u 
j 
# ar ‚| | 
— [ iue! \ > rlıdla ee) — 1 > Sttlıula &) . 
(«Al & - | u | 
JMahnke-Emde: Funktionentafeln mit Formeln und Kurven. ?. Aufl. (1935). 
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2. Das induzierte Geschwindigkeitsfeld einer unendlich langen gebundenen Wirbelfläche von 
der Tiefe /, deren Zirkulation zeitlich harmonisch veränderlich ist. Setzen wir in die Gl. (1) des 
vorigen Paragraphen an Stelle von a x, und betrachten einen gebundenen Wirbelstreifen von 
der Breite dx, und der Zirkulationsdichte /,= IT’, y(#,). so liefert Gl. (1) die Abwind- 
eeschwindiekeit dieses Wirbelstreifens an der Stelle & Durch Summation aller Wirbelstreifen 
über die Tiefe # der gebundenen Wirbellläche (wobei nötigenfalls der Cauchysche Hauptwert 
der Integrale zu nehmen ist) erhalten wir also die zur gebundenen Wirbelfläche mit der Zir- 
kulationsdichte 7’, y (#,) ei" zugehörige Abwindgeschwindigkeit: 


3 dl l 


/ f Fi 


"iu\y(a,)eitxı\ sr ie 2. 5.4 Vor SE 


vl(r \d.r. 
/ 1 | 
Hr Bi . r 2 


anw,e-tet ER 
. | 2 2 > 
Kr Ci x Vj 
Vertauschen wir ım Doppelintegral dıe Reihenfolee der Inteerationen, so ereibt sich das von 
Prandtl und Biırnbaum!')?)”) auf andere Art hereeleitete Ereebnis der Abwindeesehwindie- 
keit einer unendlichen zebundenen Wirbelfläche: 


Ht et / \ 
»(z)da emiux( | 
017 , l 1 u 
3nw,e l\‘ I iu\ \yda jetnida,da 
2, Ss l Be, 
V \ X 
(2) 
4 / 
iur 
fi \ .de\yd(e,)e ud, 
ap & ‘ i i 
\ a+1 Yj 





Ks ist leicht ersichtlich, daß die erste der obigen Gl. auch in der Form 


! { Y 
Ylu)du u '  e-iug 
aunw,e-t En . / iu \ydu)eit"\ : Ivdu 
IM > en ö a > 1 
71 il 7 ı) f] 7 
und die zweite 
! f f 
‚Ivlu)du . | e-iuv ( 
unw,a tt T,\- . "iu\ - \v(uler"dude 
> Er u “ A Een N 
„H ı) / u „H / 
4 ! 
j p / / 2 
I iu\ dv\y(u)er"du 
+ “ ‘ 
Ad > 1 
ı 71 ( 


eeschrieben werden kann. 


3. Das induzierte Geschwindigkeitsfeld einer Wirbelfläche von endlicher Spannweite und 
Tiefe, mit zeitlich harmonisch und in Richtung der Spannweite elliptisch veränderlicher Zirkulation. 

Wir betrachten nun an der Stelle = «a einen gebundenen Wirbelfaden von der Spann- 
weite 2b, dessen Zirkulationsverteilung durch 


’ b) h? nn . 
del, Je 
. 2 /, 


dargestellt sein möge. Der Gesamtänderung von I" entsprechend wird von dem gebundenen 
Wirbelfaden eine Trennungsfläche abgehen, die von der Strömung in erster Näherung in der 
x y-Ebene mit der Geschwindigkeit V mitgeführt wird (Abb. 4). 


un 


-b 
._ /y,r 





\bb. R, 


') Prandtl: Über die Entstehung von Wirbeln in der idealen Flüssigkeit. Vorträge aus dem Gebiete der Hydro 
und Aerodynamik, S. 18. (Innsbrucker Tagung 1922). Berlin (1924.) 
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Um dıe Zirkulationsverteilune der Trennunesfläche zu bestimmen. betrachten wir als 


erstes einen festen Zeitpunkt r und schreiten längs der Spannweite von y nach „+ dy fort. 
age 
. ns . , ( . 7 v. 
Dabei ändert sich / um _—dy. Der durch diese Anderung erzeugte abgehende Wirbel- 
04 \ 


streifen wird dureh den Vektor der Zirkulationsdiehte 


a/' ! ! ] | 7 
I" (#.u:r\)ds duds . PL Dur 
‚yır)ds, Y7 yds, Ei Au | yiıdır) 
bestimmt. Betrachten wir nun zwei solche Wirbelstreifen an den Stellen y= + y,, so ist 
es klar, daß sie den Schließkungsbedinguneen des Helmholtzschen Satzes nur dann ent- 
sprechen, falls von y, nach y, eine in x dichte Folge von Wirbelstreifen mit dem Zirku- 
latıonsvektor 
a r” x 
( 7 : l „, 
x I: -® - (dr } 
( 7 e. 2 | \ b] I)? u 


liegt, entsprechend der örtlichen =-Änderung der Wirbelstärke u 
Zu der abzehenden Wirbelfläche, die aus der örtlichen y-Anderung der Wirbelstärke 
des gebundenen Wirbelfadens herrührt, gehört also eine zweite Wirbelfläche mit dem Zirku- 
latıonsvektor 
of a so \r— = | h? / 


ed) ge deijdy,) Mi re n I (ide). 


yu=—b 
Der Zirkulationsvektor /"’ wird aber nun «eerade dureh die zeitliche Anderung der Wirbel- 
stärke des gebundenen Wiırbelfadens erzeuet. 


e KL (A . . ER . r. 
In der Zeit rn, =r yo. nd dr ändert sich die Zirkulation des gebundenen Wirbel- 
' (\ l dr } . r* .. 7 
[adens um \ dr, die damit abgehende Wirbellläche hat also den Vektor der 
()7 r Ar 


/arkulationsdiehte 


I" (@,n;r)ds, iT5 e rIEr I deijdy). 
. . 2 2 \ 7 e 
Aus dem Vektor der Zirkulationsdiehte des gebundenen Wirbels /'ds, und dem der abgehen- 
den Wirbeltläche /’”ds, + I"ds, baut sich somit das Wirbelsystem des schwingenden end- 
lichen Wirbelfadens auf. 
Die von dem ganzen Wirbelsystem im Punkt P(£&, 7) induzierte Abwärtsgeschwindigkeit w 


_ ( 


wird sich dureh UÜberlaererune der von den einzelnen Wirbelelementen induzierten elementaren 
Abwärtseeschwindiekeiten ereeben. Es ıst somit 





4 I)? 7° (cd &) 
Te *, / J = ‚dy 
j h | ia &) (» ’) 
4 ! 
. (A ) 7 I) y\drdy we 
/ \ \ { r - (+3). 
by I,- „ | l.ır ei m 13 ’) 
x f 
(A | 5 „ (. Sidrdy 
iu \ \v —— 
h | (. a. ’) 
| i | 


Setzen wir nun statt a ax, und d/\,—=y(rw,)da&,, so ergibt die obige Gleichung die Abwärts- 


vb’— y’ ; 


veschwindigkeit eines Wirbelstreifens mit der Zirkulationsdichte y(.r,) e ". Durch 


h 
Summation über ., von a@...a-t erhalten wir die Geschwindigkeitsverteilung eines Itecht- 
! 
ecektlügels mit der Spannweite N 25, und Tiefe f \reziproke Flügelstreckung 4 R |, dessen 


Wirbeldichte in Riehtung der Spannweite einen elliptischen, in Richtung der Tiefe dareren 
den willkürlichen y(x) Verlauf besitzt. 
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f l, 
j | l,? Tu Kr e) dr, d 7 
Itw«e \ \ (2) N ar 3 
) I iur, er 7 ) 
x [7 7 
/ FT h 
i i TH ÄN 3 VAR YVdrdydı, | 
\ \ \ y(a,)e » 1/12 2 Tr (4) 
a yet 
\ x x, 3 h 
t F I) 
! \ c.) ] I)? Ti er >) drdyd, 
\ \ \ Y\X,)‘e = —— 
bi (x Fe Di 
x X x J 


oder dimensionslos geschrieben mit a — £E=bl und n|»b für die betrachteten vier Fälle 











pP 1 p N hı 
! „ . u f ) 
Ine N „| [,iue \ : 1 ee Er FR ir 
' ] 
) Fi i 
v(e) (de ‚ u | 
ne m 7,8: "’iub\y(uJ)e'’"\ drdu 
"— | | ;r | 
f] ) H ) / Hl 
| 2 i | 
\ v(e)de ‘ \ / 
ge - - / iuh (u)e "dudev >a) 
3 \ j | \ ı ] \ (aa), 
! () f] ı) 7 / 
J 
. / hı 
rsuhl e ! »’ MH 
!.iube \ dA \v(u)e du 
I‘ ’ 
} / 7} ı 
| 
/ | u 
ine’ ww, I\ J du 
) v’ (77 1) 
Hu ] 
% l 
»1 „abe | u” ; 
",iue Tu J „dud® ee iie: 
A | u” (N m) 
! / „ ] 
fi ] 
I ul [A ( 
ao ij ] ieh / tt) 
a E P e \e \ 2 2 „du dr 
) ‚yi u" yo +ln —u) 
N / „H | 
/ ] 
| " 
ine m \y(e)(e+-nD\ a „Adnude 
.Ye+l’+n—u) 
/ u H | 
) y; ] 
| PR ORG | iubı | et duder dw 
iuhbe \y(w)e \e 74 \ — (ta). 
- | ! iM M) 
u () N / H H ] 
4 ] 
r \ (1°) \e \ »3 
| | u vVvirln in)” | 
ı) l - „H 7 ] 


wobei von den rechtsstehenden Ausdrücken der erste den induzierten Geschwindiekeitsanteil 
der gebundenen Wirbel, der zweite den der Querwirbel und der dritte den der Randwirbel 
darstellt. 


4. Im folgenden wollen wir unter der Voraussetzung, daß der Abstand x des Auf: 
punktes P von den Flügelkanten größer ist als die Spannweite, die wiederum größer sei als 
die Flügeltiefe, untersuchen, mit welcher Genauigkeit wir das induzierte Geschwindiekeitsfeld 

endlichen 


” | 
er . eebundenen Wirbellläche dure das entspreche  (rese 2 iokeitsfe 
unendlichen B ırch da ent prechendi reschwindiekeitsfeld 
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endliehen = 2 R i E 
de8 nondlichen eebundenen Wirbelfadens ersetzen können. Dabeı wollen wir den Wirbel- 
ee] Ie11e 


faden in die Schwerpunktsabszisse der y(z#) Fläche gelegt denken und die (dimensionslose) 
Abszisse « des Aufpunktes von der Schwerpunktsabszisse der y(#) Fläche ab messen. Das 


” f 
heißt, wir setzen x !-+-t#* wobeı wır ff statt j schreiben wollen (Abb. 5). 
J 


1 
4 





4 
| 





N) 
1} 
N 
R' 
Y 
1 











Abb. >. 


Betrachten wir nun der Reihe nach die Anteile der induzierten Geschwindiekeiten der 
vebundenen und @Querwirbel der Fälle 3 und 4, so erhalten wir nach Gl. (la), (2a): 


i '(r)de # /1 
.) / B 
Iıe m w.)geb.w.—=/, O\— \yew) et’ dov), X ). 
} > 3 \ ’ tr r ! hr („3 \; ) 
\ () " 
benso eilt 
anne! (9, m ,)Querw 
.() L R 
. ; hr e n N hı 
/ j f 
rub|\l u)et"bu ‚dedu dı 
\ \ j 2 f# I, \ n Hr 
7 () fh Hl q v) 
Führen wır ım ersten Integral u« — f 0 of " ein, so ereibt sich für die eckigee Klammer 
() f % % „H SL 
a « | ub 
\ (, u ha) P u P una P ! ! 
/ ( ge) — \ gr dı \ vawor dve)du / \ dı 
eu . E 4 : ! 4 ’ ! A 
H f } ı) }* ) PR () ! u" 


Nun müssen wir aber #70 gesondert betrachten. Im ersten Fall (Aufpunkt liegt vor dem 
lügel) ist der Integrand beschränkt. Integrieren wir die inneren Integrale zweimal partiell 
nach » und führen das Ergebnis in die obige Differenz ein, so erhalten wir 





% =63 L 
v . ı , -) A N hı | 
i (ru I) ‚ | | y P 
/ \ — Il) (a — Pr du|- - tr — . ‚de| 
ll y. 27 rub.r (ruba)' (iub)’. ("+ x) 
) ] 7 7 1 [5 
() 
4 | | | | \] 
r v(H ) . pi! h IH vez - er e! I) #7 ! \ d H 
\ 4? be " f (ru be)? | u Ey | 
7 - t} | E> (1 I- » 
N A 
3 „H t 
) \ ] f’ > 1. / / 
2 : Z% I yl(u)ett ou dvd 
"6 1 I): I — ır)" 
N tl E 
4 ‚ () 4 (} 
(dub) e yeba | | 
2.2. \ v(u)(u FV du ‚dv+-l',\y(w) I\du 
ne y | \ \v ta)" \ ru br Mu ! 
) ’ 7 } i } „ ; ı I+ m 
“2 
() ® [Zi ! 
. ) | E | : ’ ee ! bh I’ 
iz ä 3 -l\du —- 20T, — |\y(u) eirb w—t*) \ ‚drdu 
3 (3 l: be)? | A / e Er u I)? / (" 4 pe}? 
{ ) | l+ j H I 1" u 
d 
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r 
= 


l 
u. 





Nun gilt aber für 1 5'=F=-0 und im vorliegenden Falle «+ a | 


I # 
eipı | "dl(e!?xe) > | 
dv (0) ’ 
For m ‚r)/ Bi (+ »)/ a \/ (=) 
. I . arl 
/ Mr 
x 
d. h. es ist 
E#, 
| (1ı\ =. 
ne m, W,)Querw. =0V|75 \ y(u) (u "rd Ta 
Zu 
0 


Wie man nun nach dem Obigen leicht einsieht. ist 


)() 


ERBE et | E \ vy(u)du |, 


KL 
() 


also zılt in hinreichender Entfernung vor dem Flügel 


j: 
W; w,(1+0( |]. 
\ “IE j 
Anders liegen dagegen die Verhältnisse für e<0 (Aufpunkt hinter dem Flügel). Für die Ge- 
schwindigkeitsdifferenz der gebundenen Wirbel behält die obize Abschätzung mit x\'x ıhre 
Gültiekeit. Beı Betrachtung der Querwirbel müssen wir dagegen den Hauptwert des Inte- 
erals (x > 0) 








L x 
» vr . 4 ıbr » uh» 
e l hı ‚ pP ! ) p l ‚ı 
\ dv: -e ı#bx Jim |\ dev +\ dev) 
i ! 1 >0\, (Bi h ı 
! () 1 x 
SH \ 4 
P iuba sın ubr E iuba 
P { l X \ d } - 24 pP ! \ \ : d "Bi B= P ı fi 7 x \ d Th 
Bi ? ! 
I X () \ 
’ ! | 
ı 7 e / a N = 
ubx (u hr)? 
einführen ®). Damit ereibt sich wie oben 
2 JE E . ( I. . m, IQuerw ' 
20 % 
L rt; 17 
iub(u — f*)] er / | | 
# \ vi) \ - d Ta intemt"bx—+4 r - 0 | ‚ 
uf %: ei Fr; ıb.x (uhr) “ 
„H () ı ) 
.) () 
| | | | | 
/ v(m) | - - ei uhin Er | - - rg p! biu 5 ZT 
\; iub.e| u—I (i uber)’ u 
„H ı | \ | e 1 
| X Mr 
..3 
.. » | Co ee 
| \ yiu)(m f w da [e, = = 
N 
0 
Nun ıst aber hinter dem Flügel 
2 (0) 
J ER 
w, Oo \yao)du|e, P=3i; 
\ | u 


I) 
womit wir also in hinreichender Entfernung hinter dem Flügel die Abschätzung 


| K,\l 
ww, [1+0|e(x,+ 2)|) 
| /] 
erhalten. 


Im räumlichen Fall (bei dem noch die Geschwindigkeitsdifferenz der Randwirbel neu 
hinzukommt) kann man die Richtigkeit der entsprechenden Abschätzung wörtlich wie oben 
zeigen. 














20) ) fol a2 . i Ztsehr. f. arg 
HyAe) \ 0) Ir \V\,„ trenzlall deı ııstatlonaren rauımlıchen lraetlürels romUl}T Math. und Me 


re\y, 
ch. 





5. In der Einleitung haben wir gesehen, wie zu einem jeden Schnitt y=y, des räum- 
lichen Falles durch passende Festlegung von /, ein entsprechender ebener Fall zugeordnet 
wurde. 

längs dieses Schnittes „, wollen wir nun zum Geschwindigkeitsfeld dieser ebenen ge- 
bundenen Wirbellläche eine Zusatzgeschwindiekeit so bestimmen, daß innerhalb der Flügeltiefe 
sich in erster Näherung das Geschwindigkeitsfeld der räumlichen gebundenen Wirbelfläche 
ereeben soll. 

Wir betrachten an Hand des Prandtlschen Gedankens für den stationären Fall ver- 


suchsweise die Gesehwindiekeitsdifferenz 


Diff (m) = (w, — w,) — (mw, — w,) 


der ın der Einleitune definierten vier Fälle der gebundenen Wirbelverteilune und zeieen. daß 


im Bereich der Flügeltiefe diese Differenz in der Grenze mit 3,> verschwindet, d. h. daß 
) 


ın erster Näherung 


IE u, “7 lım (ur, ı",) 
>\ 


ilt. Wir untersuchen wiederum der Reihe naeh die Beiträge der Geschwindiekeiten, die von 
den gebundenen, Rand- und Querwirbeln herrühren. Nach der Bezeichnung des vorigen $ 


ıst also 


B, 
Ine-to?! Diff (gen W, vr) (f" ))e,(vV (# a)y)d 
() 
\ v(n) / N, ) ' m; f’ )() f’ 
.) - £ div : BE IE. DI ZZ nn Li Re mr ha 
u L u (f” 21 I, / li 
[| u” i 2 
falls wir e, (,n)=\-‚; „du setzen. Es gilt nun, wie später gezeigt werden soll 
vVa+( u)" 
($ 6), für kleine Z und |» | die Entwicklung 
2 y1 n? , 2yi 7" ” . R 
e,IC.N) - 5 fe, (n) ö’In E rc,109)%" 


die. in das Inteeral einzesetzt. für die Differenz der gebundenen Anteile 


tIre-tor Diff (ln)gen. W. 


27 | n"\y ")le,(r (f" w))In!r {2 ‚r) Pc,1® (t" z)) +. ‚|d k 


liefert. Nun kommen die Anteile der Randwirbel an die Reihe. Es ist 


() 4 
te - tor Diff(w)Randw. tuhbei"hbix ff) (ir) ee bw\e-iubr ve,(tn)drd m 
11° ( MR u’ 2 Y) 
, % 2.7 
. p iubı . | 
iub/,e' bx\e ! beve (v,n) dev FHiubet bAA 23 | n? \ dı \y(u) e' budn 
g 1 
= ! t* x) Hu () 
() () 4 
Ei Pr N ı nv A p ! hr 
N «b2] | 7\- = \yce’ h “"Audn iul,2I | ), e! bx\ d D_ 
" (f Ce). g i . " 
} „H I‘ \ 


Wir betrachten vorerst nur die ersten beiden Glieder, die wir in der Form 
2) 4 TR (?* Y) 4 
» l * 
i yl 7 ' l } E“ er . Y I) Il? 2 
rube Zr \y(o)e ee Ei EN 
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bw nach Potenzen von iubmw, so ergibt sich bis auf 


schreiben können. Entwickeln wir e’ 


(‚lieder zweiter Ordnung ın n 
'() y ’() % 
inbeitbxle-tubt"( y(e)dı-\. iub (ee dm . 
ı) ! \ u l x 
’() Hu (/ x 4 
> 1 hit \ 5 
4 ar 
\ ,\ 
7 () ! X ! X 
30 
oder. falls wir /,=b|[y(w) dw setzen. bis auf Glieder zweiter Ordnung in f 
1 
 ( H Nr Y) 
iuhei"bia ’)\y(g)eivbw\,; sumı "e,(t,n)dedw 
u u ! x 
2 () (t” \ Re 7 (ft Y) 
" pP ! { I) I’ 
) il / t .) . ar | . ul ) . 
iub et me 2y1 7 \ n) e! w \ Ir dn \ \ 
un () N x H 0 ®v (f v) 
20 u (!* x 
-\y(mw)eitbw\enim bv(e vIin v|+c,r+...)drdw 
I f] * 
1! () 1‘ = 
Nun nehmen wir das dritte und das letzte Glied der Differenz und zerlegen 
X % Zi) 4 
a - ‚& 
. e ( a = A Pe ! ! 
iubeinb(x—-t)2 y1 7? | \ de+\-- \\y(a)einbudu er \ dev) 
l f t “ . « ! 
1° 4 „H " 1 x 


Das zweite und dritte Glied dieses Ausdruckes hebt sich nach Entwicklung von e’"P“ resp. 
ei«bt* bis auf Glieder zweiter Ordnung in ft weg, ebenso das erste Glied mit dem ersten Glied 
des vorigen Ausdruckes, so daß nun die ganze Differenz bis auf Glieder zweiter Ordnung in / 


in der Form 
2.83 u u v) 
B E p iuchh®v 
Lne-tor Diff (ie) randw. iubei"b(ix— tl") 2 j | 17" \ y(mw)e'“ bw \ Indy 
’ 
1° () ! (t* 4 
’() w— (t!’—Xx) 


\ y\w)e 
‘ 


P iuh(& 5 \ pi \ y (u) e' hu d H dl | 
2 ıE K) 
Z I 


u —lf#" — a) 


dargestellt werden kann. Führen wir in das innere Integral des ersten Ausdruckes ® 


SH, 


ein und vertauschen die Reihenfolge der Integrationen, so erhalten wir 
= 13 2 U 


2.623 u () 
B i P iu h „u 2 P iu h „H 
etub(iÜ — X) \ y (m) eiubw \ z dudw= e-twb(x -t”) \ = y(m)eirdwdwdu, 
y ‚w—ıI"— x) _ \w—(f e)\‘ 
w =—) u =Q) u 0 m! [7 
das sich aber gerade mit dem letzten Glied weghebt. Bis auf Glieder zweiter Ordnung in f 
haben wir also schließlich 
267 u (6 X) 
+6, ")dedıe. 


iu h et" bix ")2J | „\y(w)eirbwl\e iuh "(e,rln ) 


v=XxX 


Izre-ttt Diff (r)Randw. — 
1° . () 


Nun betrachten wir die von den Querwirbeln induzierte Geschwindigekeitsdifferenz 


2() L 


ei bix )|\y(m)ei"bw\e iuh u e,Iv, y).d "i d 11° 


‚or Dil don IQuerw. 


Irre 
» m 
" u if” ”) 


Zu 0 
FT. 
"einbt\e nbvo.(v,n)de|, 


h 
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Et 
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Er, 
K3*. 
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Br, 
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5. 
I: 
S 
AR 
- 
x 
ar 
ur | 
Fa 
BD . 
'e ’ 
Bi 
DH %.. 
ER 
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a EEE nmuntheenn Brnseseinesuuanrnnuewinntin ner nn sammen EN 
mit IS 6) 
u ) u)du | o(c’e,) ER " 
| \yi uy. ( ) L OL zyl- le inie ia 
\ I. 


) und entwickeln den Faktor ei«»« 


Wir trennen wiederum das innere Integral ın ( 
/ 7 (f x 1 \ 
so ergibt sich bis 


des zweiten Ausdruckes nach dem lxponenten, ebenso den Faktor e“bt 


auf Glieder zweiter Ordnune ın f 
\) 
Y hHIX [® (m) p! h u p N I vo, (", 7 ) d 5 Ay. 


Diff (1° )Querw. / 
Hu () Bi X 


Nun setzen wir die Entwicklung für e, ein und spalten vom inneren Integral den Anteil mit 
” 


ec, (1+In|r|)+e, ab. Integrieren wir diesen Teil partiell nach ®, so ergibt sich 
23 { I! Y) 
pe! I) IX [ } vi) pi 7 1° p 1 hı e, ( v. 7 ) d ’h d Th 
Th ı l \ 
2 
pP rubt C, ! In f) -_ r) '(ıWw) e! 2 w d Ih 
11° ı) 
Oo 
(m)le, (m (f" ))In (f a) + c, (m (f* )) dw 
1 " 
247 / N Tv) 
iubet"d (X —}) (m)ei"bw\ei"b "(e,e In ev +c,r)dedw 
u 


Damit erhalten wir also für die gesamte Differenz bis auf Glieder ? In |" 
{re-!vtDifflw) 


20Q 2 
iubv(p vlt c)dvdw. 


\un ist nach Definition von f* bis auf Glieder zweiter Ordnune ın n 


() 7 > Y) 
et BAR *!)\v(m)eirbw\ie,dn dw 
N 
QO 20O 
GC 4 ° . . ao 
Er \yar)dmw et bH\y(w)(l+iubw) dw 2.3.33. 
t ud : ‚ s a 


also bis auf Glieder zweiter Ordnung ın f 
20 . 


"Ditfor) 2yl—n’c er da-H)|ylw)eirdw\e 
Ih u ı 5 


i«bve]n/r! de dw 


I ıe 


oder bis auf Glieder von der Größenordnune (2 2)? In (2 2) 
0 
Iıe "Ditt (me) 3v1 n° c, \v(mw) !(w — (t"— x)) In | w 7 — 2) zlnzY dw. 
/ | / \ 1} 
77° () 


Bis auf Glieder von der Größenordnune (2 @) In (2 2) zılt also für kleine 2 2 die Näherungs- 
eleiehune Diff Gr) 0 oder in der Grenze exakt 
; lim [(w, — w,) — (mw, — w,)]=0, 


(Od) —>\ 
womit also in erster Näherune die Abwindgeschwindiekeit des räumlichen Flügels durch di« 
des ebenen Flügels und der vergeleiehsweise herangezorenen Fadenfälle ausgedrückt ist. 
auf die Betrachtung der Korrektur- 


beschränken wir uns also ledielich 


Im foleenden 


dıfferenz 
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D(x,n,T)=(n ,) 
| g | 
e' , ıdı (@, 7) rube! ba | EN erh I, (a 1) 
I ) Ä 
mit 
I | } 
DD .- / 
| eet#,]) 2 
\ 4 u 
.( e-indbvydv (e-inbe 
» \j | ur \ = „du ayi n? \ dr 
0 vv ru) v 
H | \ \ 
t- ] Y. 
H \ e-tuba (7 e)de 
ee) z 3 du, 
‚yi-wW) ven —u) 
[71 ] / X 


wobei nach dem Obigen für die zahlenmäßige Auswertung allein die Bestimmung von 
lim D(#,1,r)=D(0,n,r) von Interesse ist. Bevor wir aber dazu übergehen, wollen wir 


xoN 
noch eine Umformung vornehmen. Wir erhalten durch partielle Integration des inneren 


Integrals von 4, 


Nennen wir das Restintegral P (x, uw), so ist leicht ersichtlich, daß das innere Integral 


cd 
( M : r R r 
des Ausdruckes J, dureh e; A dargestellt wird. Setzen wir also das elliptische Inte- 
® or ! 
/ 
| 
| du ER 
gral e, (#,n) = \ & „, so erhalten wir schließlieh 
| JYe’+ln u) 
avi 
| re, (x, 9) 
L 
| x 
: op ! hh)ı 
=e-!rers (2) iub\ 1 WPD (&, v)du -- 2yl y”\ = dv 
u ] / \ 
+1 
| | e or 
du 
3 de: ' ’ 
y wol 1) 
„H 


zu deren Auswertune mit 2>0 wir nun übergehen wollen. 


6. Berechnung der elliptischen Integrale. Bei der Berechnung des Grenzwertes der 
Differenz D müssen wir als erstes die darin auftretenden elliptischen Integrale bestimmen. 
Wir transformieren die Integranden durch eine lineare Substitution auf symmetrisch um den 
Nullpunkt liegende Verzweigungspunkte, wobei wir die Punkte (1,0) und (-- 1,0) festhalten 
wollen. Das liefert den Ansatz 


| 7 3 nr | ze +) + ll /.) VA 
/. co u 
I+ I+x rl] ‚)+(Il1-+% N 
Damit wird also 
14 h/ | -Aa ( 
du NZ: dr, 1—uw’= y:! L— x’), +(n- u)’ N 
a C+mn—1)/ Pe 
Setzen wir nun A() 0, so gilt = _, | 2 mt 
”-+-(n+1)’ \fi 
E> u 4 7" | 3 E m. | 
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besonders einfach werden die Ausdrücke, wenn wir statt der Koordinaten (£,n) die Para: 

ter (VW. 7 ler Kreisbüsel ae mer ZRBRE Di; I p\2 4 
ıneter (VJ,/) deı \rTEeISDLISt \e] li = 9 ö Ih rcosı) 1&c)° 0.» CONSt: R / konst. : 

Li 
| l+4 24 

MW 0, Al, ‘() s E W/; z | 2 Te ' ere nF ’ »ı 2 ın 

„ to ee ‚|; PER Fi je (ın unserem Fall « I) einführen 
(Abb 4). 

ei 
43 Bi 























Mr Mn .) en 9 MM Mm u; 
a: 
\ 
Sn 7=) Be 
Abh. hy 
Wıe man leieht einsieht, wird dann 
ö R 7 
s=Brl cos W), y=2r’(l+cosd)), kr =. 
| JAvos» Vi 
i . > r j ri Bi sın ı) 
Aus dem Dreieck ( 1. P. I) folet mit dem obieen © = „ und es ist 
| 2 Acos u — 4° 
+4: 14 | fiel yi-a 
7 Ä 0 u; zu se tzen. Damit erhalten wir nun e, (£, 1) = 2 | | \ Ä — ‚dx 
i j / - T ‘ A 
0 
pP ee 4 sın U ” N a 
mit y" te (0/2), g. Im Einheitskreis ist y >1, wir führen also «—=1-- .«? 
Y 7 (1-+ eos») 
als neue Veränderliehe ein und erhalten die Normalform 
L /\® | "du 
e,(L£,n) 21 | Y mit A = — sın (#2) 
| alvizpiiyı w“"yli—k’u ı1-+y? ade 
ı) 
’ -) 
Somit haben wir schließlich e, (£, ») ven „LE (k) KK (k)], wobei, wie üblich, E.(k) 
; E(k,na/2), K(k F(k,na/2) und k’= cos (H#/2) gesetzt ist. 
| RR ’ | u? 
Beim Integral e, (£,»7) können wir auf Grund der Identität y 1 — \i „ schreiben 
u 
| | 1] 
du | uch, u)du \ udu 
t = — ] > — = ) : u 
7 "vet ”) 1 "It 1)” "Jyl uw VL (1 tt) 
| 
Andererseits ist aber auch 
I | 
| vl—wdu | Yl- Win u)’du 
€ 5 } .. 6 
x JyYsrin Br h yvrn tt) 
Ber | ) 
BR, 
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\ ) H d | | 
Fassen wir F Fr | „| zusammen und integrieren das zweite Integral 
+h 1)” du\yZ@+ln uw) 


1 


partiell, so erhalten wir 


H du 
“2 t \ ' 
1 u’ VL’ ) ) 
. s N ’ s ö * (ce, (£,n)) j 
Es ergibt sich also einerseits die Bestimmungsgleichung für e,: e.(E.n). 


andererseits aber, falls wir die obige Gleichung in die erste Umformung von e, einführen. 


u H In 


P Se, t| = = =. 
‚yi BY Hr 1)” ‚1 “"yv+t 1)’ 


aus der sich leicht die Zurückführung von e, auf eine Normalform herleiten läßt. 
Das erste dieser Integrale liefert, vermöge der anfangs eingeführten linearen Trans: 


j ai dr 
lormatıon =. - — , das letzte 
q | 1 | Ky’ır 
1 | 
7 en” l dr 
7 \ | es yaeyı ey va 
| | 
L/ | dx N da 
| \ (9 )] J 
ı)vi-eyltya q3U+nyi—aryl+r?a 
Demnach wird also 
ei 
L/ | | 7 l dr 
Ge, —+(l ) ——, | 


Da das Integral über den ungeraden Integranden verschwindet, erhalten wir nach Substitution 
von W =] er 


| du 
E Ce, +2 Bd \ De £ 
y |] I+y’ (1+-Wuw)yl--wyl Fu 
( 
mit 
., 1) > ( 
u? Es "= sın (9]2). 
i | 7" 
Führen wir nun «== sin # ein, erweitern mit yI - A? sin? und setzen für den Integranden die 
en A RB 
Partialbruchzerlegung "Der a ya >.:,,2 34 an"), So ergibt sıch 
1—k’sm?i 1-+u?sın?! 
i dt h? . u” | | k? sın? # 
\ Zn —— —K (kl) + ———— | il 
„A+w#sın’t)y1l—k’sın’t A’ u @+k?) 1+ u sın? 
ı 0 


Das letzte Integral ıst eine der gebräuchliehsten Normalformen des vollständigen elliptischen 
Integrals dritter Gattune”). Demnach haben wir also 


= kEI — u" o ge 
E, Ce +2 ie | ;K(k) 
2 (a Ik® Pi u _ IR a 5 
li u ’ k’? Ei r s) ‚ ) j \ ’ ’ ) r ' ’ ) 
> D N 5) (K (Fk) E (k , P) T lH; (k) F (k A P) IN (I) F (h r P) ‘ | 
ek RT w- 5 


‘s), Hamel: Berechnung des vollständigen elliptischen Integrals dritter Gattung. Sitzungesberieht der Berl. math. 
Ges. 31 (1932). 

»), Nyström: Praktische Auswertung von elliptischen Integralen dritter Gattung. Soeietas Seientiarum Fennica 
C.P.M. VII. 12 (1935). 
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1 
} 4 | 
7 eote NM A E s ] COS PD, 0 } Im 


| | ’ sın ) | | ) 
wofür wır auch zusammenfassend schreiben können 


ke h' | a nyl- 
e Ä P -. (1 7)” ] k) + > 
| | Ks > Vkötıy k'® 


| 1? IK (k) E(k', pP) -E (k) F(k’ pP) K (1) F(k’, B)\ 


’ 


Im Hinbliek auf den nachfolgenden Grenzübergang mit Z[>0 wollen wir für die Integrale e, 
und e, schon an dieser Stelle die Entwicklung für kleine Werte von £ vornehmen. Zu diesem 
/,weck betrachten wir erst den Zusammenhang zwischen k,%k’, {,n. Aus 
Ir - g” | i & 5 .) 
cos I >67 und li ”+(ıl+n)%, " + (l 7)" folgt 


vos ı) 
(1 1 ) | ar R | \j - 
"TH 1 
also ım Quadrat [ j | 
„ & a; | 
KIM, > 
| Ä Ä I+1 ER min 
| } | [ "> KiE® | 
' (1 ) (1 1°) 
/ Ä 
Dementsprechend Ist 
| } 43 
ge 4 - = 
cos| " { | gr (] n”)” | | , N | 
sın J 3 | r n° 2° 
(1 )Il—+ —(?ı ._0(d£°) 
(1 ı) (1 ) \* | 
oder 
" 1} ( j ı) ® , 
J: COS, ‚(1—0,({?)), k zu si + (1 — 0,(6Y)) für SV. 
u / pen 


Für e, erhalten wır demnach 
,=- TEC) KR (k)]. 


\lıt 0 eeht 27”’>0, k>1l,. und es eelten die Entwicklungen *) 


FH (4) I+1/2(4 ID) Er (K?A),. N (k A+((k’?A)). 


| | 
wobeı | In|,,) gesetzt Ist. 


Ks wırd also sein 


d.h. limöle, ©.n) =—-|=lım (( A))=0O. 
„>| - .—>|\ 


.) | 7° FR 
et CA-+((k'? A) 


Beim Grenzübereane im Inteeral e, betrachten wir erst den Ausdruck in der geschweiıften 
hlammer. Es ıst 
lımı +47 E (1. rz[/2) F(V,p) t- Jım F(k, 2/2)’ D(k',d), 


—>- | — | 


da Di, p) beschränkt ıst"), wird 


lım en) = ATeCOs lım ((£? .D)) = arccos N. 
<i>) ee 
Ks ıst also 
lim le, (£,n) +2 y1l—n’Inl 271-9? (l Int(l-y’))+2naresıny. 





Band 18, Heft 6 


Dezember 1938 


Vs 





Demnach ergibt sich also lım A, («, 
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)=V und wie aus der Entwicklung des Integralsinus 
x>\ 
resp. Cosinus ohne weiteres folgt 
R . - = E ve i bı 
lim Je-i"bxXe, (x N, 2] | 7” \ dv 
x—>0| 
ı „N 
lim (e, (@&,n)+2y1l- n’In,='+2yl- |lnuby-+inl2 (x)))) 
x > (|) 
2y1—1?(—-1+1n4 (1— 1?) uby) + 2n aresin „+ inayloı, 
wobei In y die Eulersche Konstante bedeutet. Somit haben wir also 
() U iorle ul 
D(O,y:r) [># 7 I-tubA,(0, 1) 1,0, »,)] 
mit 
1,0, ,)=2yl—y’( -I+ln4(l = mn’) uby) +2naresinn 
Hl 
+iryl- mn’ -tub\yl- "PO, u)du 
| 
Hl 
(On) | u PO, 1) ’ 
19 = ei; 
. ‚vl—-w _d(n—w) 
u 
P iuhv 
7. Die Ermittlung des Integrals ®. Führen wir in das Integral ? - \ vv? „A » die 
v”rinNn 1)“ 
! e \ 
neue Veränderliche w= ubr ein und setzen ub() - u)=a, so erhalten wir 
J- 
e pP im 
Dx—0:3 D. (a) N D, (a) \ — „dw. 
Jva’ rw 
1? 0 
Nun gilt aber bekanntlich für die Hankelsche Funktion nullter Ordnune erster Art die 
Integraldarstellung 
S- 
cosmw dm 7 
P.(a): - s=i— H'(ila|) 
() 
womit also zur Bestimmung von ® 
untersucht werden muß. 


nur die zahlenmäßige Darstellung des Imagınärteiles 
Bilden wir die erste und zweite Ableitung von ® nach dem Parameter a und führen 
den aus der ersten Ableitung gewonnenen Ausdruck 


Iddß 
ada 


FI. SI. 
ld® ‚[e-iw.wdw a ne” 
i\ 5 +3 \e-ie - „dw 
a da )Va’+w’ ya’ + w‘ 
u 0" 
über, so ist leicht ersichtlich, daß die Gleichung 
SI. 
ED Idd \ı | 
— nn i\ pP vw (d| = 2] 
da? a da \ Va’+w’ 
0 
gilt. Nach partieller Integration des rechtsstehenden Integrals erhalten wir für ® die 
Bestimmungsgleichung 
‚Frr | 4 ».. 4 .. n Z ! | ’ | 
DPD’ + ya -D- ‚ also für den Imaginärteil 2,’ + Ps D,= 
( ( 


ad 





EB nin 
Mn 
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mit der Nebenbedingung P, (a=0) =. Da die homogene Del. durch Besselfunktionen nullter 


Ordnung mit rein imaginärem Argument befriedigt wird, setzen wir die Lösung in der Form 
P,=e,(a)J,(ia)+c,(a)i H,' (ia) 


1 
„Variation der Konstanten“ muß dann sein 


an. Nach der Methode der 
4 ® m. 1 . . ! . f} 1 ° | 
C, J/, (ta) rt HH, (ra): U). ıc, J, (ta) > cC, HM, (ia)= 
| " dl 
») 
Mit Zuhilfenahme der Funktionalgleichung®) J,_, (a) H,' (a) J,H,_, ergibt sich also 
id 
y . f ’ . .T e 
C, - l > I, Ga), C, > J.(ta), 
d.h. 
#1 A 
j wi IT | . . ‚ ° . ‚ . „T . 
P,=\-\i5 Hl’ (ie)dc+ KN,)J, (ia) HS, Ga)dce+K,li- H,(ia). 
\ 7 \ () 
Da mit a >V das erste Integral wie a?’Ina, das zweite wie a? gegen Null geht, J, da) >| 
H,' (ia) dagegen logarıthmisch unendlich wird, muß A, =,, K,=V0 sein. 
1) 

















Abb. 7. 


Mit diesem Ereebnis können wir also schreiben 


, 6 nr 
D, ar MH (ub 7 ") 


dd 
} P JI . . x D D JT - 
$D, | \; Fn H'(ka)daid), (iu l 7 tt |) I\J,lia)dzliz Hin I, 7 |) 
() iu Y ii pr 
7 
oJ. ub 7 u), 
mithin 
Ö D, 7 
- tub— H,'(iub|n u) ub|+ Ye (lub 7 u )| 
Ev, ") ‚ re 
(dl 
\D g 7 7 
NY? = + ub|| \i5 H(kia)da)i I, (tubiy wm +(\9, Ge) dr)5 Hub 
l $ : 0 
. „7 . 
! J (tub u)=ub|+g (ub|n u )], 


») 


u. 


) 





ew. 
ch. 


ter 


rm 


so 
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wobei a für ub(n —u) steht und das obere resp. untere Vorzeichen zu wählen ist, je nachdem 
(n— u) =0O gilt (Abb. 7)'%). Auf Grund dieses Ergebnisses können wir nun eine für die nach- 
foleende Auswertung der Restintegrale sehr vorteilhafte Umformune vornehmen. Wir setzen 
nämlich 


+ ] 7 ıIl+n 
(1) yl TR! (ub|ı u) du \r\ \yl (n xy” D (uba)dı 
I l N Kl 
1-1 
\J | (7,4 " Pl(ubx) do (| |» I). 
xl 
Es ıst also 
| 7 | } 
d j ) 4 5 (v I 7 
\ = Du br) dr \ = D(u br) d.r 
dı vl—(n—x)’ ‚vi-(n+x)’ 
Nun können wir aber auch schreiben 
] N) ] 7 
d f . 
\Pl(uba)dy | (3 2 \Pfu bx)d] | (n + .r)* 
dn ’ ’ ' 4 / 
Y ı) X “ ( 
5 I+n = | n 
= | | ("7 x” D (u h .r) «>07 | | (N + .r)” D (u hr) TER 
1 ) 1 ! 
„d® „daD 
. \ ] | (3 2A 7 dx \ | l (n x)’ my dx. 
() 0) 


An den oberen Grenzen liefern die ausintegrierten Teile den Wert Null. Zur Ermittlune des 
(‚renzwertes an der Stelle #=0 entwickeln wir die beiden Wurzelausdrücke nach x. Es ist 


26) 5% da”). 


Da P bei 2=0 höchstens logarıthmiseh unendlich wird, ıst der Grenzwert des Produktes. 
somit also der ausinteerierten Teile an der unteren Grenze ebenfalls Null. Wir haben also 


l —? 
d? f \ 7 ge Ad» ] \ N % da» 
n - GT = Ar 
d 17“ I | (» av dx ‚I l (17 + x)? dx 
X u () 
| 
u dPDıiub| 1] m) f H dP(ub 7 ") 
\ ; / du \ = du 
yi u di - u) y1 u din- u) 
H | H ! 
l 
A u ü dP(ub N rt) 
F} (7) \ M) du 
‚yl u“ din ) 
l 
oder, falls wir /=ub, hn,4)=2yl—n?(—-1+ln4(1- n?)y4)+2naresin n setzen in 
endgültiger Schreibweise 
1 75% Ze 
MD) (o,n:r) eor -_(A+ib) 
ha h 
N 2 f ( } 7 - 
. o 1. . ( (1,4) 
A=inyli—n’— Afeln,A)+ — ee ee er ee m. . 


Q 7” 


0° fi. (7, /) 
0 7” 


RB —— /h (7,4) + 4 fs (77, 4) 





10), Die Zahlenwerte zur Abb. 7 wurden mir von Herrn R. Müller (FIB) freundlieherweise zur Verfügung 
gestellt. 
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8. Die numerische Auswertung der Phasenfaktoren. Bei der numerischen Auswertung ist 
es notwendig, erst den Variationsbereich der Parameter festzulegen. Statt der Parameter t, b, u 


| ! : 
führen wir die reziproke Flügelstreckung 4 57 und die Wellenlänge der Schwingung 
Pr SU 
In In 
BF \ ein. Auf Grund unserer Voraussetzungen für die Herleitung der 
oO [Ei - 
pr A £ ns 
Gleichung (A) setzen wir J=V... !Yı und L=x... S=2b. Dann wird / ] RR 
4 
Also müssen wir die Größen A,B der Gl. (5) in den Bereich für =0...3,7=0,0...1,0 


hinreichend genau bestimmen. Die Form der Darstellung von A und B läßt die graphische 
Auswertung vor allen anderen als sehr zweckmäßig erscheinen, vorausgesetzt, daß in dem 
interessierenden Bereich der Parameter eine hinreichende Genauigkeit erzielt werden kann. 

K"alls wir F, (7) mit dem Integraphen durch Umzeichnung des Integranden für feste 
‚und ») Werte ermitteln, kann der Mittelwert der Ablesungen im Durchschnitt mit ca. 0,15 '/o 
Teilstrichen falsch sein. Falls E,, E, die Einheiten des Integranden auf der &- und y-Achse 
in em bedeutet, y der in Teilstrichen abgelesene Integralwert ist und nach der Eichung p 


rın . . a . Due . rg‘ . . HK Y 
eilstriche — x em? entsprechen, haben wir für die erste Integralkurve y Teilstriche — — "7; 
pExEy 
* Iyı , \4p 
Kinheiten. Der relative Ordinatenfehler der ersten Integration ıst also : - 4 
p 


Nach einmaliger instrumenteller Integration von F}(n) haben wir den Fehler der dureh Ay 
und .Iy auch bei genauen Integrandenwerten bedingt ist (>03), und dazu den Fehler, 
der dureh die fehlerhaften Werte des Integranden herrührt (> 0,30). Das gleiche gilt auch 
bei der 2. Integration. Das Doppelintegral der F-Kurven ist also möglicherweise um ca. 0,9" 
falsch. Da die Bestimmung von f noch eine Integration erfordert (Festlegung der Integrations- 
konstante), kann f um ea. 1,2°/o, F um ca. 0,3°/o falsch sein. Es ist also der Gesamtfehler 
ungefähr mit 1,22°+0,3"/o proportional. Bei kleinen Parameterwerten (bis =1) ist also 
die zu erwartende Genauigkeit der graphischen Bestimmung befriedigend. Bei wachsendem 4 
steigt der Fehler quadratisch mit /, da aber der wesentliche Bereich für 4 zwischen 0...2 
(im extremen Fall bis 3) liegt und in diesem Bereich nach dem Obigen ungefähr mit einem 
Fehler von ea. 5°/o zu rechnen ist, wollen wir diesen Grad der Ungenauigkeit auf Grund 
der sonstigen wesentlichen Vorteile der graphischen Methode ın Kauf nehmen. Daß bei 
Werten >53 sehr starke Streuungen auftreten müssen, ist nach dem Obigen klar, trotzdem 
wollen wir bei der Bestimmung von A und B die Lösung nicht nur für die Werte = 0,2; 
0,5; 1,0: 2,0: sondern auch für 4 == 4,0: 6,0; 8,0: 10,0 ansetzen. 


/ 
glAX)-g, (x) / Bd 
1 7 
ei 





[A 


Als erstes betrachten wir F,(y)=4\ »9'n- u )du, wobei nach dem vorigen 


(Is T r , : en ä 
y | zahlenmäßig bekannt ist. Durch Umzeiehnung der Kurve 4 (x) auf die Einheit 
FA P7 . 
erhalten wir 9 (Z.#) = 9; (#) (Abb. 8) oder mit z=n— wu 9;,(n-— u). Nun führen wir u = cos t 
‚T 
ein. Damit geht das Integral über in 4] 9; (7 — cos PM) d (sin #), d.h. übertragen wir die Funktions- 


werte g;. die zur Abszisse »; - cost zugehören, auf die Abszisse =sin ft, so ergibt sich 


- 
c 


PD p(E)AE, wobeır g;,(£ die auf die Abszisse E umgezeichnete 9; (x) Kurve bedeutet. 


(0, 1,0) 





eW. 


ech. 


ist 


ng 


ler 


'n 
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Das mit dem Integraphen auszuwertende Umlaufintegral über 9; (£) liefert den Wert 
von F} (n) für ein festes 4 und n. Liegt für ein festes 4 die 9, (x) Kurve gezeichnet vor, 
so ıst für verschiedene mit der Einteilung der Abszisse  übereinstimmende Werte von n die 
Umzeiehnung der 9; (x) Kurve für alle diese 7 Werte nach Abb. 9 ohne weiteres durchführbar. 














Abb. 9. 


Im Falle y, ist die Umzeiehnung (da 9, beschränkt ist) ohne weiteres durchführbar, 
im Falle 9. müssen wir das Integral (wegen der Singularität des Integranden bei u „)) 
erst durch Abspaltung der Singularität der graphischen Behandlung zugänglich machen. 


Wir setzen 


| ’ 
+ Q.(Alı |) + R.(/ | ul), G “od 
JS / 1 (7 ) ( / / 
und betrachten 
E.7 | | 
4 \ " \ u du \ H n 
/ sg. du = + 4 (+ B.,)d4W. 
| ML | un er | u? i 
] 
| 
Das erste Integral ergibt nach Einführung von Zu =e'!—+te ‘t, eil 
r » 2 . ns y 
[A du | x 4] dz | >20: \ ER 
Jvil-w n-u 2:J #"—2nz+1l z 2: Br 
l (0) 


der Integrand AR. des zweiten Integrals ist beschränkt und wird auf zeiehnerischem Were 
wie g, ausgewertet. 


Der Verlauf von F ist in Abb. 10 und 11 dargestellt. Da nun aber nach dem vorieen 


dFf df 


; Fi) mit 0 
d 1)" AN d 7 n ) 
gılt, erhalten wır 
n 
fi) \A,)dar H2\yl-wWPoiwdu. 
0 0 


Das zweifache Integral wird für feste 4 Werte nach Aufzeichnung der F},(n) Kurven mit dem 
Integraphen bestimmt. 
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w) J; U DY/; F 1-0 
Yu u. All 
S Vf-ı ofn U) 
7 | 2 
, \ + + + In 
| 
\\ 
f: )® \ 
7 > 
| 
\bb. 10. Abb. 11. 
Im letzten Inteeral setzen wir u sın?t und erhalten 
12 
. r . 2 sin ?2/ 
2\y1i- uw Piin) du \ P(isindddlt +). 
N t 0 . 
= ‚sın2t ER ie Ai an 
Führen wir also &=t+ „als neue Veränderliche ein, so ergibt sich [ PD, (£)dE, wobei 
. () 


- 


PD; (E£) die auf die Abszisse £ umgezeichnete P (A sin ft) = P; (sin f) Kurve bedeutet. Im Falle 


D, lieet der Umzeichnung, da P®, beschränkt ist, 


MH ) 
’ 


nichts im Wege. @. wird dagegen bei 
En 


logarithmisch unendlich, wir zerlegen also das Integral in (+ f. Der Zahlenwert des 


0 € 


ersten Integrals wird durch Reihenentwicklung, der des zweiten durch Umzeiechnung vermöge 


des Integraphen gewonnen (Abb. 12 und 13). 








Abb. 12. 






+ A=05 


| JVru?-G-Am-u)du 


7 











\bb. 13. 
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9. Die Korrektur des Anstellwinkels. In der anfangs angeführten Arbeit von Kassner 
und Fingado°) ist für das ebene Problem des harmonisch schwingenden Flügels der An- 
stellwinkel als Funktion des Geschwindigkeitssprunges der Unstetigkeitsfläche berechnet. 
Bezeichnen wir die letzt angeführte Größe mit @=u, ei’, so gilt 


' 
I 


nVa=ü r. | - u a i| 3 +f; | 5 u | |, 


"\2 u 


> 


wobei die Bedeutung von f,,. aus der unter °) zitierten Arbeit zu entnehmen ist. 


Indem wir den Geschwindigkeitssprung durch die Gesamtzirkulation des ebenen Ver- 
eleichfalles ausdrücken. haben wir 


Er [A ’ u | 
i F = 1+utfl5 “t)\Hiutfels ut). 


Nun gilt aber weiterhin nach dem vorigen die Winkelkorrektur des instationären Falles an 
der Stelle 7 
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Abb. 14. Abb. 15. 











Abb. 16. 


(Abb. 14, 15 und 16), woraus mit /,=y1-—n?” und der in der Tragflügeltheorie üblichen 
Orientierung des Winkels gegen die Anströmrichtung der ungestörten Strömung 
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folgt. Führen wir statt «a die Wellenlänge Z ein und setzen 5, 2, 50 ergibt sich 
) 
/ 


JT 


u [> BL }© ' 0 1 
4 (1+2| AEIA L, 2)\)yl 7 zZ l a)! 


Apr: 


„(7 


+2) rl 2)vi-® 58 


Da mit a>V afsta), af.la), Bia,n)>V und Ala,n)>— 7 geht, wird für endliches 


S,Q und L>x 
I’ | : f 
ARäumlich» stationäı 7 vt| | | 7 ’ ) N 
Setzen wir noch elliptischen Umriß t=#,y1-—n? voraus und führen das Seitenverhältnis 
F I, 7 


—=A °” ein. so erhalten wir 
Pre Ya | N 


ARäumlich» statıonär 7 Vt 


d. h. das Prandtlsche Ergebnis des räumlichen stationären Falles. Unsere Gleichung stellt 
also für den instationären Fall in erster Näherung eine Erweiterung der Prandtlschen 
Gleichung für endliche Spannweite dar. 


In Abb. 17 ist der Anstellwinkel als Funktion der Zirkulation ke 7 - F(w,4%,n) 








ae ’ f Eur „ | Zu 
speziell für den Fall 7 =0 skizziert, wobei 2=,, und n die reziproke Geschwindig- 
SU J 
keitszahl’) bedeutet. 
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Abb. 18. 


Übernehmen wir aus der unter?) zitierten Arbeit in erster Näherung den Auftrieb als 
Funktion der Zirkulation, so können wir analog der Kassner-Fingadoschen Darstellung 
den Auftrieb (ebenfalls in erster Näherung) als Funktion des Anstellwinkels durch den 
Auftriebsvektor 7’ (@, 2,1) kennzeichnen. In Abb. 18 ist P speziell für den Fall 7=0 


skizziert. (2==0 entspricht beidemal dem Grenzfall der unendlichen Spannweite. 562 
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Zur T'heorie der Wirbelbewegung. 
Von A. Kneschke in Dresden (Sa.). 


ie Theorie der Bewegung isolierter Wirbel ın ebenen, einfach zusammenhängenden Flüssig- 

keitsgebieten ist durch Arbeiten von Routh, Lagally und Masotti') zu einem ge- 
wissen Abschluß gebracht worden. Bei allen diesen theoretischen Überlegungen wird von 
dem Helmholtz-Kirchhoffschen Grundsatz ausgegangen, daß die Bewegung jedes 
einzelnen Wirbels nieht von seinem eigenen Wirbelfeld beeinflußt wird, sondern daß der 
Wirbel sich vielmehr wie ein regulärer Punkt der nach Wegnahme seines eigenen Strömungs- 
feldes verbleibenden Strömung bewegt. Zahlreiche Abhandlungen beschäftigen sich mit der 
Auswertung der Routh-Lagallyschen Theorie und prüfen an Beispielen ihre Leistungs- 
fähigkeit?). Dabei hat sich ergeben, daß der Helmholtz-Kırchhoffsche Grundsatz ein 
geeignetes Hilfsmittel zur Untersuchung der Bewegung von Wirbeln in begrenzten Flüssig- 
keitsgebieten darstellt. Während die Kirehhoffschen Untersuchungen sich nur auf den 
Fall beschränken, daß die freien Wirbel sich in der unbegrenzten Ebene bewegen, so haben 
Routh und Lagallvy in den zitierten Arbeiten ein Verfahren angegeben, das ermöglicht, 
die Bewegung isolierter Wirbel in einfach-zusammenhängenden Gebieten zu beschreiben. 
Von der speziellen Wirbelbewegung in der Halbebene ausgehend, hat M. Lagallv mittels 
konformer Abbildung die Transformation auf solche Gebiete durchgeführt. Im ersten Teile 
der vorliegenden Untersuchungen sollen die Bewegungsgleichungen der isolierten Wirbel 
allgemein aufgestellt werden. Weiterhin wird das Strömungsfeld der Wirbel energetisch 
untersucht und die Transformationsvorsehrift der kinetischen Energie des Feldes aufgezeigt. 
Und im dritten Teile soll schließlich ein Minimalprinzip angegeben werden, das den Bewegungs: 
gleicehungen äquivalent ist und in energetischer Hinsicht Aufschlüsse über den Bewegungs 
voregang der freien Wirbel in abeeschlossenen Gebieten geben kann. 


Die Bewegungsgleichungen. In einem ebenen, einfach zusammenhängenden Flüssigkeits- 
gebiete @ mit der zeitlich unveränderten Umrandung $ sollen sich an den Stellen 2, =a,-+ib, 
n freie, isolierte Wirbel mit den Zirkulationen 7, befinden. Die Umrandung 5 des Gebietes G 
ist als Stromlinie Y/= const aufzufassen. 

Wie M. Lagally in der zitierten Arbeit erwähnt, lassen sich die Bewegungsgleichungen 
der Wirbel dadurch aufstellen, daß man ihnen in der unbegrenzten Ebene eine im Inneren 
von @ reguläre, zeitlich veränderliche Potentialströmung überlagert, die die Normalkomponente 
der Gesamtbewegung an der Begrenzung S zum Verschwinden bringt. Wir können demnach 
als Stromfunktion für das Gebiet @ ansetzen: 


I. a A # yr 
lg (a) + y— dr)? N Ile - a”? +(y- bi°+ \ 


4a u) EN — 


! l (dh) 
+flry: a; bj) | 


wobei also unter f(x, 9; a;,b;) der in @ reguläre Anteil der Stromfunktion zu verstehen ist, 
der durch das Vorhandensein des zen Wirbels bedingt ist und dazu beiträgt, die Umrandung S 
als Stromlinie auftreten zu lassen. Y ist — kurz gesagt Greensche Funktion des 
Differentialausdruckes 4 Y für @, woraus ihre Symmetrieeigenschaft f(#, 9: &»)=fl&, 1; :, 9) 
offensichtlich wird. Daß der dem Aten Wirbel entsprechende Anteil von der Summe ab- 
getrennt wird, ist analytisch, wie wir sehen werden, zweckmäßig. Nach dem Helmholtz- 
Kirchhoffschen Grundsatze werden die Komponenten der Geschwindigkeit des Akten Wirbels 
lauten: 


3 fly; a,b) | 


da; Of ur} \ gi | 
le (a, — a, ° + (b, — bb)’ + (a, b,; a; b;) 
At Q b;; | Te | JT r u + r ) u —— / ‘ . 
i i 
N 
+ f (2y: a; b,.) 
R „ P ix a] 
Yy bh; 
1) E.J. Routh: Some applieations of eonjugate functions. Proc. London Math. Soc. 12, 1880, 73.— M. Lagally: 
Über ein Verfahren zur Transformation ebener Wirbelprobleme. Math. Z. 10, 1921, 231. A. Masotti: Sul moto di 


un vortice rettilineo. Atti della Pont. Accad. d. Se. Nuov. Lincei 1931, 235. Sopra una proprieta energetica del moto 
di un vortice rettilineo. Atti della Pont. Accad. d. Se. Nuov. Lineei 1931, 464. Sulle azioni dinamiehe dovute ad un 
vortice rettilineo. Atti della Pont. Accad. d. Se. Nuov. Lincei 1931, 3. 


2), E. Paul: Bewegung eines Wirbels in geradlinig begrenzten Gebieten. ZAMM 14, 1934. A. Kneschke: 
Über die Ruhelagen von Wirbeln. ZAMM 14, 1954. W. Barth: Wirbelbahnen um Wände und Platten. ZAMM 10, 
1930. A. Kneschke und S. Matthes: Wirbelbewegung um einen Kreiszylinder. Ann. d. Phys. 9, 1931. Vom allge 


meinen Verfahren etwas abweichend behandelt G. Hamel in der ZAMM 13, 1953 die Wirbelbewegnung um eine Buhne, 
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Ab, Of u £ , i  Tgfe 
] \ lo (a, (t;)” | (hy b;)” + f (dt b,.: a; b;) 
dt od, | PER In — | 

! ! 
\ 
( » 
\ f Br „: (lı h;,) 
I A x da] 
y—=bx 


Da ferner als Folge der Symmetrie von f die Beziehungen 


‚ N 
\ 


/ N 
N = Ei. y,a,D) Dur „ub; dt, h\. 


») 
x d _ 0 h 
7) bh 


fix.y; a,b) fia,b: a,b) 
velten u können die 


Geschwindigkeiten dureh partielle Ableitungen der sogenannten Routh- 
schen Funktion 


", a2 | 
Ay \ le (a, — a,” + (b; b;” + \ fla;, bi; a;b)+ 5 flarbr; ab). . (ID 
ee T — _ 
ın der Form 
d (tl, Q .e d b, Q N, II] 
dt Ob, ’ lt 5. En tn | j 
dargestellt werden. 


Kinetische Energie des Strömungsgebietes A. Die Form dieser Bewegungsgleichungen 
läßt einen Zusammenhang der Routhschen Funktion mit der Strömungsenergie des Feldes 
vermuten. Schaltet man in bekannter Weise durch Kreise 0; 
eulären Stellen z, (k =1,2 


zu) 


mit den Radien e, die sin- 
.n) aus, so zerfällt das Gesamtgebiet @ in kreisförmige Wirbel- 
eebiete W, und in ein Restgebiet R (Abb. 1). Vorerst soll für das Strömungsgebiet R die 


Abh. 1. 





kinetische Energie der Flüssigkeit von der Dichte o berechnet werden; sie ist gegeben durch 


() 


Zi Pr’ + PN)dadn.:. . . > 2:2... 
ir 


woraus durch partielle Integration folgt: 


» 
( 


Oo he e 
a > \\vawaway- 5 Pd WiP,de-Pdy)+5 \ PiP,dae— Pedy). 
E R A S :: IR | . 


h 

Wirbelfreiheit von R aber verschwindet das Flächenintegral; das Integral über die 
Berandune S aber liefert einen Beitrag zur Energie, weil auf der Begrenzung des Gebietes 
VW = const ist. Nach Normierung auf Y: 

Berechnung von T, der Ausdruck 


Wegen der 


0 verschwindet dieser Integralwert, so daß zur 


" 
Le; 
[#) \ 


Zu - J . IV ‘ 
= Yv(Y,dc — Y,dy) .(la) 
— — « 
IR | 0 
k 
) N) un R (< 12.9: 5,9 TiE.n: 5, 
E lim n 
8a x=i JSi>0 u 
‚ HU J 
andererseits Ist: 
EEE TE. 1; Ifks+JI5n;:5 +4I5,n) TIE.M3-8, 9) 
= 1m £ 
ds JE> 0 2% 


woraus sich dureh Entwiekeln des Zählers ergibt: 


öfis,n;$5.N) fIlE+ JS,n: 5) (6,0; 5,9) orlx,9; 5, 
- > lim 2 


Us 1:20 45 ÖX 
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zur Verfügung steht. Nach (I) ıst demnach die Integralsumme, wenn man den hei 27 regu- 
lären Anteil der Stromfunktion einstweilen mit r (x, y) bezeichnet, 


a. Be A fr, dr.) 
1,= > P \’ (x, 9) ‚„'s (re — a,” + (y — bj.) r ay®? 3—d 7 
BER 0 


n 
0 \ , N, N | I, n ol (X a) d Y% (Y b; ‚dx 
nr: (2,9) le (@ — a,” + (y - b})” ar = 
De ER, ' 9 Br | IT d | (x a; + (y b;)’ 
ul, l 


auszuwerten. Setzt man 2 =a7;,+e;, eosy und y=b;-+ter,sing für die Integration längs 
der Kreise so, mit ihren Radien &),. so erkennt man, wenn nach dem Mittelwertsatze 


[Zi n 
y 
Be : y 
u > Ig la — a)? + (bu bi] + N fta,; br; a;b)+flar; br; a; bu) Wi; 
x Aa JI 
; . Fi rd 
7 bi ! 1 1 ) 


vor das Linienintegral gesetzt wird, daß das erste Integral für beliebiges &, verschwindet 
und das zweite Integral den Wert 
[7 


TR—2 \'ln | Pen I 


ETO 2 lo #;. 
JT 


u EN 2 
k—| 
oder, ausführlicher geschrieben, 
4 „ 
Fi ”. 
j . 2 
IE? = RR \ | \ le (a; (t;)“ t (bh; b;)”] 
2 Lu hr nn I 
k=1 i=1 


(IV) 


» k 
+ f (a; b,: a;b;) + f (a; br; a; bu) — 5. |g er), 





erhält. Um nun den Zusammenhang der kinetischen Energie mit der Routhschen Funktion 
aufzuzeigen, braucht man nur ihre partiellen Ableitungen unter Berücksichtigung der Symmetrie 
von f zu bilden; also: 


oT ‚0Xu f Abu OT f Au r des 
9) E 5) #) 
R . At Obu ia P O bu « e At 


Die Bewegungsgleichungen lauten dann 


N dd Aıyı 0 T7; r d Du ( T,; V 
Mr m . IV), 
A dt Obu ET Rz? 
wobei 
n 
\'lı Ekı’,\ 
"RR (ıy 
7} = 0 > | Pr I lg Er) 
l 1 


zu setzen ist. Multipliziert man die erste der Gleichungen (V) mit db.„, die zweite mit da, 
so folgt nach Summation beider als Integral der Differentialgleichungen der Wirbelbahnen 


a 


y 


worin der Satz von der Erhaltung der kinetischen Energie der Strömung während der Bewegung 
zum Ausdruck kommt. Schreibt man die Energie in der Form 


n 
Y '& 2 
en k \ 
ar it Var Eee \ ——]8 
t = \ 1 n 1 n Sr? | 
k=1 
„ 
h | Fi, u.4, ee u 
mit (2: 27 \ /,; U’, so erkennen wir in 2 die Kirchhoffsche Bahnfunktion'), 
IT u 
k= | 


4) Dieser Zusammenhang ist bereits von E. Trefftz erkannt worden und findet sieh ebenfalls in der an zweiter 
Stelle angegebenen Arbeit von Masotti. 
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die also im wesentlichen die kinetische Energie des Strömungsfeldes unter Ausschluß der 
Wirbelzentren darstellt. Sie ist abhäneize allein von der gegenseitiren Lare der Wirbel und 
von der Begrenzung des Flüssigkeitsgebietes. Wir wollen diese Energie T, deshalb als 
„Lageenergie* des Wirbelsystems bezeichnen. 
Es ist wichtig zu wissen, wie die Energie T,; im Gebiet @ mit der komplexen Ver- 
änderlichen 2=x2-+iy dureh konforme Abbildung auf ein Gebiet @, mit der Veränderlichen 
&+in dureh die Funktion z2=f({) sieh transformiert. Bei einer derartigen Abbildung 
werden auch die singulären Stellen in @ in solehe in @, übergzeführt: insbesondere «eeht ein 
in «an der Stelle 2, = a, + ib, gelegener Wirbel ın einen in @, an der Stelle £,= a; —+# ß}. 
E] dz; FR 
I (Sr) 


eeleeenen über. Wir erhalten demnach aus IV, weıl lım 


€] dl GC 


0 


[7 


Pe hr ,; N, .r - , 
IR, 2 JI d \ 4 Pin m lo / (Gy ) I. le Er] 


un” In 
I, | 
oder 
n 
Ur, l,; 0 P} N, le f’ (<;) m . . ° . . ° . . . . . (VI | R 
l | 


wenn man in dieser Transformationsvorschrift in 77, mittels der Gleichung a, + ib, =fla;, ti Pr) 
die Koordinaten a, und b, durch «a, und p, ersetzt. Für das Gebiet @, läßt sich demnach 
die Energie in der Form 


; N 
Tr, ) 7 OD \ Ih | P° u 5 lo € ko 
% u “T PER 
h l 
mııt 
I} | 
wo. y FAN = le / (£x) 


ansehreiben. 


Ein Variationsprinzip. Die meisten Differentialgleichungen der mathematischen Physik 
stehen in enger Beziehung zu Variationsprinzipen. Daher liegt die Vermutung nahe, daß auch 
die Bewegungsgleiehungen für freie, isolierte Wirbel sich als Eulersche Gleichungen aus 
einem Minimalprinzip mit Hilfe der Methoden der Variationsreehnung herleiten lassen. — Offen- 
bar lautet in unserem Falle dieses Grundgesetz der Wırbelbewegung, wenn man 


„ 


ö } B) \ Balcı da, db, \ op . y. 
einen Ausdruck ’) T, 5 I. \b, 14 1 | als „Energiefunktion* des Wirbelsvstems 
einführt, 
t, 
fo 


oder: IneinemeeerebenenebenenStrömungsfelde bewegen siıch freie, isolierte 
Wirbel so, daß das Wirkunesinteeral 


f, 7 

. ) O \ Y i d (ty db; ee. ) r 
/ -\12 l, Ib; dt (l); dt In, dt . A . j . . (IX) 

fü h l 


vegenüber allen benachbarten, virtuellen Bewegungen, die in demselben 
ZJeitintervallvon der Ausgangslage des Wirbelsystems zur Endlage führen, 
stationär wird. Kürzer formuliert: Für die wirkliche Bewegung der freien Wirbel ist der 
zeitliche Mittelwert der Differenz aus der „Energiefunktion* T, der Wirbel und der „Lage- 
energie“ 7, des wirbelfreien Strömungsgebietes stationär"). Nach den Vorschriften der 
Variationsreehnung ergeben sich hieraus sofort die allgemeinen Bewegungsgleichnungen 
2 z dd, Q Tl; a d hu OÖ I; 
ol dt Obu” ol dt Od 
die mit (V) übereinstimmen. Damit aber ist die Richtigkeit dieses Minimalprinzipes bereits 
bewiesen. 333 


5) Die hydrodynamische Deutung des Energieausdrnekes 7» ist bisher nieht einwandfrei gelungen. 
6), Für Ratschläge bei der Abfassung der Arbeit sei Herrn Prof. Prandt] herzlichst gedankt. 
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Über Lüders-Hartmannsche Linien. 
Von Endre Reuss in Budapest. 


1. Allgemeines über die Fließgrenze. Bei der Formänderung von Flußstahl und einigen 
anderen Materialien beobachtet man im Augenblick, wenn die Spannungen die Fließgrenze 
erreichen, an der Oberfläche das Erscheinen von Fließfiguren. Diese Fließfiguren, auch 
Lüders-Hartmannsche Linien genannt, bilden zwei nahezu orthogonale Liniensysteme, 
welche überall in die Richtung der Hauptschubspannung fallen. Ihre Anwesenheit kann an 
der ursprünglichen Oberfläche, aber auch auf Schliffen durch das Innere des Körpers durch 
die Frysche Atzmethode sichtbar gemacht werden: sie sind die Spuren von Fließschichten 
an der jederzeitigen Oberfläche. 

Die Frysche Ätzung greift vorzugsweise die am stärksten verformten Stellen des 
Schliffes an und zeigt, daß das Material der Fließschichten eine große bleibende Verformung 
erlitten hat, während das sie umgebende Material praktisch unverändert blieb. Man schließt 
hieraus, daß die Formänderunge an der Fließgrenze nieht stetig, sondern in den Fließschichten 
stattfindet. 

Dieser räumlichen Unstetiekeit der Formänderung kann eine zeitliche Unstetiekeit der 
Spannung gegenübergestellt werden, die Unstetigkeit an der Fließgrenze. 


= 


Wird ein Zerreißprobekörper aus Flußstahl in der Versuchsmaschine allmählich ver- 
formt, so steigt die Spannung nahezu proportional zur Formänderung an, bis sie an der 
„oberen Fließgrenze* plötzlich auf einen niedrigeren Wert, auf die „untere Fließgrenze*“ 
herunterfällt. 

Die beschriebene Unstetiekeit tritt nicht nur beim Zugversuch auf, sondern kann auch 
beim Biegungs- oder Verdrehungsversuch wahrgenommen werden, indem bei diesen Versuchen 
die Belastung an der Fließerenze oft etwas fällt. 

Betrachten wir den Verdrehungsversuch an einem kreisrunden Stab etwas näher. Der 
allgemeinen Gepflogenheit folgend nehmen wir an, daß die Spannung auf den zur Stabachse 
konzentrischen Zylinderflächen konstant ist. Dies folgt aus der axialen Symmetrie der 
Formänderung und aus der Annahme, daß die Spannung r(y) eine eindeutige Funktion der 
Schiebung y ist. Bezeichnet 9 die relative Verdrehung zweier Querschnitte der Entfernung 1, 
r den Halbmesser irgendeines der erwähnten Zylinder, R den Halbmesser des Stabes, so 
\ r 
ist y: und das Drehmoment der Spannungen 


I 


R 
M=\2nr’t(y)dr. 


Differentiation und darauffolgende partielle Integration liefert 


R 
dM 7 r R’t(y) M 
») so & a2 . .) / BD} 
1) ala, r? dı An y Pe 
woraus 
| M.,. 
(= H3M| a Ba GO Ge | 


folgt. 

Will man mit Hilfe dieser Gleichung aus dem Drehmoment-Verdrehungsdiagramm den 
für das Material charakteristischen Zusammenhang zwischen der Schubspannung und der 
Schiebung, 7 (y), berechnen, so erhält man an den Stellen, wo das Drehmoment fällt, für 7 so 
kleine Werte, daß diese mit anderen Versuchen nicht in Einklang gebracht werden können. 


} RN ‚.dM | er 
Ja. schon an einer Fließerenze mit 19 0, bis zu der sich der Körper ideal elastisch ver- 
j ! 
hält und dem Proportionalitätsgesetz folgt, müßte 7 (y) von seinem Maximalwert auf dessen 
« rg‘ . ” . “1. fi I M * 
a Teil herabfallen, da aus dem Proportionalitätsgesetz 15° M folgt. 
Fr 


Für den Biegungsversuch ergeben sich ähnliche Verhältnisse: die aus dem Biegungs- 
versuch berechneten o(e) Werte stehen mit den unmittelbar aus dem Zugversuch erhaltenen 
im krassen Widerspruch. Daß sie im späteren Verlauf der Formänderung miteinander gut 
übereinstimmen, deutet darauf hin, daß unsere Vorstellungen über die plastische Form- 
änderung im Grunde genommen richtig sind und nur für die Fließgrenze einer Berichtigung 
bedürfen. 
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Den Weg aus diesem Widerspruch weist uns die unmittelbare Beobachtung. Das Auf- 
treten der Lüders-Hartmannschen Linien zeigt, daß unsere Annahme betreffend der 
Gleichheit der Spannungszustände auf koaxialen Zylindern falsch ist. Auch bei Berechnung 
des Momentes beı der reinen Biegung eines rechteckigen Stabes wurde bisher immer an- 
genommen, daß in den zur neutralen Ebenen parallelen Ebenen überall derselbe Spannungs- 
zustand herrscht. Aus dem Auftreten von Fließfiguren folgt, daß man in diesen Fällen mit 
einem viel verwickelteren Spannungszustand rechnen muß. 

Ks ist allgemein bekannt, daß das Auftreten von Fließfiguren eine Materialeigenschaft 
ist, die nieht auf eine Beanspruchungsart beschränkt ist. Wir schließen hieraus, daß der 
Mechanismus der Entstehung von Fließfiguren von der Beanspruchungsart weitgehend un- 
abhängig ist. Gelingt es, für einen speziellen Fall die Entstehung der Gleitschichten, die 
Spannungsverteilung und das Fortschreiten der Formänderung in alle Einzelheiten zu ver- 
folgen, so glauben wir uns berechtigt anzunehmen, daß der so beschriebene Mechanismus 
auch für andere Belastungsfälle nahezu derselbe ist. 

Um überhaupt Berechnungen durchführen zu können, muß ein idealisierter Körper ins 
Auge gefaßt werden. 

Wir nehmen an, daß sich unser Körper für große Formänderungen ideal plastisch ver- 
hält, in dem Sinne, daß die größte Schubspannung immer gleich einer konstanten kritischen 
Schubspannung %k ist. Für kleine Formänderungen soll sich unser Körper ideal elastisch ver- 
halten und das Proportionalitätsgesetz befolgen. Der Erscheinung der „oberen“ und „unteren“ 
Klıeßerenze tragen wir dadurch Rechnung, daß wir dem Vorgehen von Nakanishi, Prager, 
Odqviıst') folgend annehmen, daß die Schubspannung im elastischen Gebiet auch über den 
kritischen Wert k steigen kann. Die Körperteile mit >% befinden sich in einem labilen 
Gleichgewichtszustand. 

Die Frage, unter welchen Umständen ein labiler Spannungszustand bestehen kann, er- 
scheint als eines der wichtigsten der noch ungelösten Probleme der Plastizitätstheorie. Am 
einfachsten ist die Annahme einer maximalen elastischen Hauptschubspannung. Nach Prager 
hängt dieser Höchstwert auch von der Körperform ab. Es ist möglich, daß für die Stabilität 
des Spannungszustandes außer den Spannungskomponenten deren räumliche Abgeleiteten 
maßgebend sind. 

Glücklicherweise spielt dieses Problem in unseren Betrachtungen keine Rolle, so daß wir 
diese Frage ruhig offen lassen können. Es genügt uns im folgenden die triviale Tatsache, 
daß an den Stellen, wo sich für die elastischen Spannungen unendliche Werte ergeben, kein 
elastischer, sondern nur mehr ein plastischer Zustand herrschen kann. 

Im UÜbergangszustand vom ideal elastischen zum ideal plastischen nehmen wir an, daß 
der Körper aus räumlichen rein elastischen Gebieten besteht, die von flächenhaften ideal 
plastischen Gebieten durchsetzt sınd. 

Die flächenhaften Gebiete, die Gleittlächen, welche die Fließschiehten endlicher Dicke 
idealisieren, sind Unstetigkeitsflächen, so daß Punkte, die auf beiden Seiten derselben zu- 
einander unendlich nahe waren, sich nachher um eine endliche Strecke entfernen. Die ganze 
plastische Formänderung ist also in diese Gleitflächen konzentriert. Auch die Spannungs- 
komponenten erleiden an den Gleittlächen Sprünge, aber so, daß sie die Gleichgewichts: 
bedineuneen überall erfüllen. Sie müssen in den Gleitflächen außerdem der Fließbedingung 


N , eenüeen. 


2. Qualitative Betrachtungen über den Verdrehungsversuch. Wir wenden uns nun dem spe- 
ziellen Beispiel zu, an dem wir alle Einzelheiten der plastischen Formänderung an der Fließgrenze 
verfolgen wollen: der Torsionsbeanspruchung eines Stabes von beliebigem Querschnitt. Im 
rein elastischen Gebiet wird der Spannungszustand durch eine einzige Spannungsfunktion 
erschöpfend beschrieben und durch das Prandtlsche Membrangleichnis veranschaulicht. Im 
ideal plastischen Gebiet läßt sich die Spannungsverteilung ganz ähnlich beschreiben, indem 
das Mempbrangleiehnis zum Sandhaufengleichnis ergänzt wird, so daß das Ausweichen der 
Membrane durch eine Böschungsfläche über dem Querschnitt begrenzt wird. 

Da die Schubspannungen in der Querschnittstläche Hauptschubspannungen sind, besteht 
das eine System von Gleitflächen aus Querschnittsebenen, das andere System steht hierzu 
senkrecht, miuß also aus Zylinderflächen bestehen, die zur Stabachse parallel sind. Die Quer- 
sehnittsebenen sind in diesem Falle miteinander offenbar nicht mehr gleichwertig, was eine 
beträchtliche Verwieklung des Problems mit sich bringt. Um dies zu vermeiden, lassen wir 
das System von Gleitflächen, welches aus lauter Querschnittsebenen besteht, ganz außer Acht 
und stellen so die Gleichwertigkeit aller Querschnitte wieder her. Ob dies zulässig ist, läßt 


1) Vel. W. Prager: Mecanique des solides isotropes au delä! du domaine «lastique. (Memorial des sciences 


age 


mathematiques.) Paris, Gauthier-Villars 1937, wo auch ausführlichere Literaturhinweise zu finden sind. 





ngeWw. 
Mech. 


Auf- 
der 
nung 
an- 
IngS- 
ı mit 


chaft 
; der 
| un- 
die 
ver- 
smluS 


Tr Ins 


| ver- 
schen 
| ver- 
eren” 
Der, 
r den 
ıbilen 


l, ©X° 

Am 
acer 
ılıtät 
iteten 


b wir 
sache, 
‚ kein 


ı. daß 
ideal 


Dicke 
n Zu- 
vanze 
1ungs- 
rıchts- 
neung 


11 Spe- 
renze 
t. Im 
nktıon 
t. Im 
indem 
7] der 


‚esteht 
hierzu 
(Juer- 
IS eine 
en wir 
r Acht 
t, läßt 


sciences 


Band 18, Heft 6 en E Enz 2 
Dezember 1938 Reuss,. Über Lüders-Hartmannsche Linien 349 





sich bestreiten. Es läßt sich dagegen einwenden, daß die einzelnen Systeme von Gleitflächen 
nicht einzeln untersucht werden können, weil plastische Spannungszustände nicht super: 
ponierbar sind; daß beim Verdrehungsversuch plastischer Körper der Bruch stets als Gleit- 
bruch in der Querschnittsebene stattfindet und daß durch die blockierende Wirkung der Ein- 
spannköpfe eben das Zustandekommen des zur Stabachse parallelen Systems gehemmt wird. 
Trotz dieser Einwände müssen wir das aus Querschnittsflächen bestehende System doch 
außer Acht lassen, da wir sonst überhaupt keine Berechnungen durchführen könnten. Die 
jerücksichtigung dieses Systems muß als eine weitere, noch zu lösende Aufgabe betrachtet 
werden. 

Die Spuren der Gleitflächen auf den Querschnittsebenen nennen wir im folgenden Gleit- 
linien. In den Gleitflächen ist die Schubspannung zur Stabachse parallel und hat den Wert k. 
3erücksichtigt man die Symmetrie des Spannungstensors, so schließt man daraus, daß die zur 
Gleitlinie normale Komponente der Schubspannung in den Querschnittsebenen ebenfalls den 
Wert k hat. Folglich werden die Gleitlinien im Membrangleichnis durch Kurven mit kon- 
stantem, der Spannung k entsprechendem Gefälle dargestellt. 

Bei starker Verdrehung, wo das Sandhaufengleichnis angewandt werden kann, sind die 
Linien mit dem Gefälle k die geraden Erzeugenden der Böschungstläche über dem (@uer- 
schnittsumriß. Die Gleitflächen sind also Ebenen, welche zur Mantellläche senkrecht stehen. 
Wir nehmen an, daß dies auch nahe zur Fließgrenze so ist. Diese Annahme wird auch durch 
ddas vorhandene reiche Versuchsmaterial?) gestützt. 

Auf Grund unserer Annahmen können wır nun folgenden Satz aussprechen: Die Mem- 
brane, welche den Spannungszustand nach Überschreitung der Fließgrenze veranschaulicht, 
schmiegt sich einem Skelett der Böschungsfläche an, welche aus lauter Erzeugenden besteht. 
Die Schubspannung wird hierbei stellenweise größer als k ausfallen. 

Wir sind nun in der Lage, die Ausbildung des plastischen Zustandes Schritt für Schritt 
zu verfolgen. Dynamische Einflüsse sollen nieht berücksichtigt werden. 

Laut unseren Annahmen kann die Schubspannung im rein elastischen Zustand die untere 
Fließspannung k überschreiten, so daß in Teilen des Querschnittes ein labiles Gleichgewicht 
besteht. Dieses labile Gleichgewicht soll an irgendeiner Stelle, vermutlich dort, wo die Schub- 
spannung am größten ist, aufhören. Einem bekannten Satze der Potentialtheorie zufolge tritt 
dıe größte Schubspannung am Umfang auf. An dieser Stelle fällt die Schubspannung auf %k 
herunter. Es entsteht eine Störungsstelle, die ähnlich wie eine Kerbe wirkt und zum Vor- 
dringen der Störung führt. 

Dieser Vorgang kann am Membrangleichnis veranschaulicht werden. Von der Stelle der 
größten Schubspannung an wird die Membrane durch eine sich nach dem Inneren verlängernde 
Erzeugende der Böschungsfläche niedergedrückt. Auf diese Weise entsteht die erste Gleitlinie. 

Wann kommt das Vordringen der Gleitlinie zum Stillstand? Denken wir uns die Er- 
zeugende der Böschungsfläche genügend verlängert, so gibt es einen Punkt, bis zu welchem 
die Membrane sich der Erzeugenden anschmiegt. Wir behaupten, daß sich die Gleitlinie 
nau bis zu diesem Punkte ausbreitet. Den strengen Beweis werden wir später bringen. Daß 
sich die Gleitlinie nicht weiter ausbreitet, ıst offenbar, da die Schubspannung über diesen 
Punkt hinaus den Wert k nicht erreicht. Man könnte aber denken, daß das Fortschreiten 
der Gleitlinie schon vor diesem Punkt haltmacht und hier eben auch ein labiler Spannungs- 
zustand bestehen bleibt. Wir werden später zeigen, daß dies nicht der Fall ist, da dann das 
Gefälle und mit demselben die Schubspannung im Endpunkt unendlich wird und deshalb die 
Störung des Gleichgewichtes weiter fortschreiten muß. Wir gelangen so zu einer Erklärung 
für das plötzliche Erscheinen der Gleitlinien. 

Durch das plötzliche Auftreten einer Gleitlinie wird das Gleichgewicht des ganzen 
Körpers gestört. Wır betrachten zwei extreme Fälle. 

Einmal halten wir den Verdrehungswinkel fest. Im Mempbrangleichnis bedeutet dies, 
daß der Druck unter der Membrane konstant gehalten wird. Durch das Auftreten einer Gleit- 
linie, was durch das Aufdrücken einer Erzeugenden der Böschungsfläche veranschaulicht wird, 
wird das Volumen unter der Membrane, also auch das mit demselben proportionale Drillungs- 
moment, verringert. Das Auftreten einer Gleitlinie entlastet daher plötzlieh den Körper. 

Das andere Mal halten wir das Drillungsmoment, also das Volumen, unter der Membrane 
konstant. Wird eine Erzeugende der Böschungsfläche auf die Membrane gedrückt, so muß 
das verdrängte Volumen anderswo Platz finden. Die Spannungen werden an anderen Stellen 
erhöht, das Gleichgewicht hört hierdurch auch an anderen Stellen auf. Es entstehen neue 
Gleitlinien, neue Erzeugende werden auf die Membrane gedrückt, der Druck unter der Mem- 
brane steigt, der Verdrehungswinkel nimmt zu, der Körper fließt. 


e- 


5) 


2, Nädai, A.: Der bildsame Zustand der Werkstoffe, Berlin 1927. — Mac Gregor, Ü.W.: Recent Investigations 
in Plastie Torsion. Journal of Applied Mechanics. 4 (1957), p. A—163. 





e 
we 
* 








LET Sn 





or ' z un Ztsehr. f. angew. 
O0) Reuss, Über Lüders-Hartmannsche Linien Math. und Mech. 








Die Erseheinung des Fließens von verdrehten Stäben, die bei Versuchen gut beobachtbar 
ist. hat so eine natürliche Erklärung gefunden und ist mit der Bildung von Fließschichten 
in einen engen Zusammenhang gebracht worden. 

Kehren wir zum ersten Belastungsfalle zurück. Wird der Körper weiter verdreht, so 
treten auch hier weitere Gleitlinien auf. Da die bereits vorhandenen auf ihre Umgebung 
entlastend wirken. treten die neuen immer womöglich entfernt von den schon vorhandenen 
auf, Ist das System der Gleitlinien schon genügend dicht, so treten die neuen Gleitlinien 
immer in der Mitte zwischen zwei vorhandenen auf. Die Gleitlinien zeigen ein Bestreben, 
sich äquidistant zu verteilen. Für den Kreisquerschnitt muß demnach die Anzahl der Gleit- 
linien eine Potenz von zwei sein. 

Bei fortschreitender Verdrehung wird das Svstem der Fließschichten, das Skelett von 
Erzeuzenden immer diehter, der Spannungszustand nähert sich asymptotisch dem aus der 
Plastizitätstheorie berechneten Spannungszustande. 

Bei der zweiten Belastungsart erfolgt die Verdrehung in Sprüngen, sonst verläuft alles 
ebenso wıe beı der ersten Belastungsart. 


3. Berechnung der Spannungen. Wir setzen die Lage der Gleitlinien und den Verdrehungs- 
winkel 9 des Stabes von der Länge eins als bekannt voraus und stellen die Aufgabe, hieraus 
die Spannungen zu berechnen. Bezeichnet F(x,,y) die Ausweichung der Membrane, so genügt 
sie im rein elastischen Gebiete der Differentialgleichung 


1? F a er 


mit der Randbedineune F=0 am Querschnittsumriß. In den Fließschichten ist überall =k, 
also berechnet sich daselbst F=ks, wo s die Bogenlänge der Gleitlinie, von dem Umriß an 
zerechnet, bezeichnen soll. Könnte die elastische Spannung unbegrenzt erhöht werden, so 
hätten wir auch bei der Verdrehung % einen rein elastischen Spannungszustand: die Aus- 
weichung der Membrane, die mit Hilfe bekannter Methoden der Elastizitätstheorie berechnet 
werden kann, würde dann F*(r, y) betragen. F* genügt ım Inneren des ganzen (Juerschnittes 
der Gleiehung 

j? F” 2G0 Er a Go ar Gar Sr Ve .. (>), 


am Querschnittsumriß der Randbedingung F*=(). 
Setzen wir 


Re. a er ae Ben u Si 


so genügt u der Laplaceschen Gleichung 


Li 
-_ 


7 a rt ne | (>), 


erfüllt am Querschnittsumfang die Randbedingung «—=0 und in den Gleitlinien die Rand. 


hedinzung 
en ee 7 te ee ET Aa: 


Unsere Aufeabe ist auf ein Problem des lograrıthmischen Potentials zurückeeführt. Es 
steht uns frei. die Gleitlinien mit zum Rande zu rechnen. oder sie als Unstetiekeitsstellen zu 
betrachten. Wiır wollen letzteres tun. 

Bei der Membrane verhält sieh F in der Gleitlinie stetie. Bezeichnet » die Normale 
\ F 
( 
RT 
daß in den Gleitlinien eine einfache Belerune » anzenommen werden muß. An allen anderen 


zur Gleitlinie, so weist da einen Sprung auf. Dasselbe eilt für u. Wir schließen hieraus. 


Stellen verhält sich die Potentialfunktion regulär. 

Für die nun folgende Untersuchung setzen wir die Greensche Funktion Kür, y: £&,n) 
als bekannt voraus. Wäre auch die Größe der Belegung in den Gleitlinien bekannt, so 
könnte man 

ul, y) riE, IR. .dEM 77): 7 A 
berechnen, wo die Integration bei konstantem “x und y über alle Gleitlinien auf die Variablen 
S(s), n{s) zu erstreeken ist. Da entlang der Gleitlinien der Wert von « aus (6) bekannt Ist, 
läßt sich für die Belegungsdichte die Integralgleichung erster Art 


u 77 E53 & | re Eee ee . (8 


anschreiben. Nach der Berechnung von » aus dieser Gleichung ergibt sich aus (7) und (4) 


die Gleiehunz der Membrane 
F—-F* I\rKds a a er 
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und weiter die Schubspannungskomponenten 
OF \ OK 
T T y’ { N 
07, . 07, 
(10) 
04 OK 
Ty z 7 \ y— ds 
VER j Vu 


mit Integralen, die im Stetiekeitsbereich gewiß konvergent sind. Die mit einem Stern be- 
zeichneten Größen beziehen sich auf den fiktiven ideal elastischen Zustand mit dem Ver- 
drehungeswinkel 9. 

Die Berechnung des Spannungszustandes ist also auf die Bestimmung der Greenschen 
Funktion und auf die Lösung von Integralgleichungen erster Art zurückgeführt. 


4. Die Gleitlinien. Wir wollen jetzt den Sprung, den die Spannung und die Verschiebung 
an den Gleitlinien erleiden, berechnen. Wir rechnen jetzt die Gleitlinien mit zum Rande und 
bezeichnen ihre beiden Seiten mit 1 und 2. Bezeichnet ferner » die innere Normale, so ist, da 


IF* aF 


0) 
On, On, 
| u O7, oF of " 
E m u |, 
on, On on, On, 
oder, da 
OH 
Ts 
on, 
Ts, + Ts, y Gr . (12) 


Sinngemäß ist der positive Sinn von s, und s, auf den beiden Seiten der Gleitlinien einander 
entgegengesetzt. 


Die Normalkomponente von r beträgt laut unseren Annahmen 


In, In. I: VE 5 A ern SE 


Es darf uns nieht befremden. daß der Absolutwert von r an der Gleitlinie größer als k aus 
fällt. Es herrscht da ebenso ein labiles Gleichgewicht, wie wir am (Quersehnittsumfang an- 
ceenommen haben. 

Zur Berechnung der Verschiebung in den Gleitebenen erinnern wir uns, daß dıe Ver- 
wölbung des Querschnittes durch die zu 


F’ r? 
7 ee Te GE (4) 
(11) 2 
konjugierte Funktion 9 angegeben wird. 9,  y, ergibt dann die relative Verschiebung ın 
\ N 
a . vg 01) en. 
der Gleitebene und kann durch Integration von | ‚ entlang der Gleitlinien um den 
Us un 


Endpunkt herum erhalten werden, da am Endpunkt der Gleitlinie g, 9,=0. Aus (14) folgt 


mit Rücksicht auf (4) 
Oy 1 /d F* O7 o /r 


F | 
on G3\on Ö | on\2 | 
Bei der Integration fallen die stetigen Glieder aus der Differenz 9, g, heraus, und es bleıbt 


l I 
07 | (du RT 
DD, ( un di DE wu 
di R. \ OS N On On, 


ds, 


wo im zweiten Integral nur mehr einfach bis zum Endpunkt zu integrieren ist. Beachtet 
man (11), so folgt weiter für die relative Verschiebung in der Gleitebene 
l 
ı f | 
G ’ q (; R \ vVitiS (5; y 


a ee eh . (15), 
I 
die also zum Integral über » proportional ist. 

Die Kraft. die von der Membrane auf die Erzeugende der Böschungsfläche ausgeübt 
wird, beträgt, auf die Längeneinheit bezogen, 


‚[d F 04 16 
g=BI7 Tau usa 4. al BE Sa Zur Sa Be Be Dee Ba ee au Ze Ze °___ 
1 2 


wo S die auf die Längeneinheit bezogene Spannung der Membrane ist. 
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5. Der Endpunkt der Gleitlinien. Wir sind noch den Beweis schuldig, daß jede einzelne 
Gleitlinie bis zum Punkte fortschreitet, bis zu welchem die Erzeugende der Böschungsfläche 
die im Mempbrangleiehnis die Gleitlinie veranschaulicht, sich der Membrane anschmiegt. Wir 
werden zeieen, daß das Fortschreiten vorher nieht zum Stillstand kommen kann, da sonst das 
(sefälle der Membrane, also die Spannung, im Endpunkt unendlich wäre. 


Im Punkte. bis zu welehem sieh die Erzeurende der Böschungsfläche der Membrane an- 
(\ F N 
schmiegt, ist, ebenso wie am Querschnittsumfang, g=%), FE To -V und daher auch v —=V0. 
( ! ( 
u 


In den Punkten, in denen die Erzeugende die Membrane niederdrückt, ist 49 >0, also auch 
y 0". 


Wir untersuchen nun, wie sich die Membrane in der Umgebung des Endpunktes einer 
Gleitlinie verhält, für welehe mit Ausnahme des Querschnittumfanges » >V ist. Zu diesem 
Zwecke legen wir den Anfangspunkt des Koordinatensystems in den Endpunkt der Gleitlinie, 
die positive -Achse in die Gleitlinie. Aus = 2 |vKds, wo die Summe über alle Gleit- 


Yu ’T OK 
= \ \ y- ds, was für die Punkte außerhalb 


(0) X (0) X 


linien zu erstrecken ist, berechnen wir 


der Gleitlinien zulässie ist. Wir zeigen nun, daß diese Größe über alle Grenzen wächst, 
wenn .e von den negativen Werten der «-Achse her dem Anfangspunkte zustrebt. Für K 
können wir 


| 
h u log r-+y1%2,95%,N) 
schreiben, wo r’=(.r — &)°+(y-— n)” und y mit allen Ableitungen bis zur zweiten Ordnung 


stetie sind. Für genügend kleine Werte von x sind mit Ausnahme des untersuchten Integrals 
alle Integrale der Summe gleichmäßig beschränkt. Das Integral, welches sich auf die unter- 
zu u 4 2 ee | Ö . 
suchte Gleitlinie bezieht, beträgt \v(&)|, 7 dz. Das Integral über das erste 
a J/TI& = 0) 31 
N / 
Glied zerlegen wir in eine Summe |=|-+]|. Wegen v(0) >O und der Stetiekeit von v kann 
ı ı 
man immer ein e finden, so daß für &<e vi£)> m >V ist. Der zweite Summand ist gleichmäßig 
beschränkt, für den ersten ergibt sıch 


m l € 
ri) - log 
\ inc —E 2 ‚| ) & en x 
() m 
ou 
läßt man . dem Werte Null zustreben, so strebt dieses Integral und damit \_ _, welches 
Ox 


davon nur um eleiehmäßie beschränkte Glieder verschieden ıst, dem Unendlicehen zu. 


6. Die Gestalt der Gleitlinien. Wir machten die Annahme, daß die Gleitflächen Ebenen 
sind, welehe je eine Erzeugende der Böschungstläche des Sandhaufengleichnisses enthalten. 
Mit Hilfe des Castiglianoschen Prinzips läßt sich eine Gleichung herleiten, welche an den 
(leitlinien erfüllt sein muß. 

Sind auf einem Teil der Oberfläche eines Körpers die Spannungen, auf einem anderen 
Teil die Verschiebungen gegeben, so sagt das Castiglianosche Prinzip aus, daß ım Falle 
des Gleichgewichtes die Ergänzungsarbeit 


A#=-A—-IJ)miä+tovHrnZ)do.:. . . : 2 2 502... MM 

N 
ein Minimum wird. Es bezeichnet hier 4 die gesamte Formänderungsarbeit des Körpers, 
x, », w die Verschiebungskomponenten, =, H, Z die Komponenten der auf die Oberfläche 


wirkenden Spannung, do das Oberflächenelement. Die Integration ist nur auf jenen Teil der 
Obertläche auszudehnen, wo die Verschiebungen gegeben sind. Bei der Variation sind nur 
solehe Spannungszustände zugelassen, welche die Gleichgewichtsbedingungen erfüllen und an 
den Teilen der Oberfläche, wo die Spannungen gegeben sind, mit diesen im Gleichgewicht sind. 
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Beachten wir, daß im Falle eines verdrehten Stabes die Spannungen auf der Mantel- 
fläche und in den Gleitflächen vorgeschrieben sind, daß für die Länge eins die Form- 
änderungsarbeit 


A: G \ \(Fy? FF )dxdy, 


und daß im zweiten Glied von (17), wenn wir den Stab an dem einen Ende eingespannt denken, 
am anderen Ende Z=F,, H Fr... 2=0, u y, e=„)%x ist, so erhalten wir nach 
partieller Integration und Berücksichtigung der Randbedingung F=0 


y rz \\l ä m i 2G0 Fidady, 


wo die Integration sich auf den ganzen Querschnitt zu erstrecken hat. Im Gleichgewicht 
soll A* ein Minimum haben, also 


B !- „ nn F,° 
o\\l 


.) 


G 6 4* 2G9F\dady=0 . ..... (8 


A) 

sein, wobei für all&@ zum Vergleich herangezogene Funktionen am Querschnittsumfang F=0, 
in den Gleitebenen F=ks ist. Diese Randbedingungen führen auf der Oberfläche zur 
Spannung Null, ın den Gleitflächen zur Grenzspannung k. 

Wir gehen nun zur Berechnung der Variation (18) über. Bei dieser Variation ist, da 
wir die Gleitlinien mit zum Rande rechnen, außer Fauch der Rand zu variieren. Es wird daher 

| It FE+ RR} 

GoA—=\\ö| B > n 


Ja NS 


I/F£+fF)/ 
260 F\dzay+\\\| - = =—2G9H \dedy |; ° 


Im ersten Glied erstreckt sich die Variation nieht mehr auf den Rand. Im zweiten Glied 
haben wir die von der Variation des Randes herrührende Variation besonders angeschrieben. 
Das erste Glied ergibt bekanntlich, wenn wir den Rand mit R bezeichnen, 
(I(F.. F,„t2@V)öFdxzdy-+ |öF-( F.dy F,d«). er 
S R 
im zweiten können wir die Inteeration nach der Normale ausführen, indem wir bedenken, 
daß wir die Unstetiekeiten der Ableitung von F zum Rande gerechnet haben, also F stetig 
differenzierbar ist. Wir führen Koordinaten ds und dn ein, wo ds das im positiven Umlauf- 
sinn gerechnete Bogenelement, dn das Element der inneren Normale bedeutet. An die Stelle 
des Flächenelementes dx dy tritt so dsdn. und es folet mit Rücksicht auf die Stetirkeit 
von F,, F,, F 
(Fe +Ff,y 
\| ; 9) - 2G0fF \ö nds ; \ 3 i 3 ; 21). 
R 
In den Randgliedern drücken wir alles durch s und n aus. Bezeichnet 9 den Winkel, 
den die in die positive Umlaufrichtung weisende Tangente des Randes mit der „-Achse ein- 
schließt, so ist in (20) 


F.dy F,de=(F,sing F,cosy)ds= F,ds. 


Da der Betrag des Gradienten von F gegenüber .der Drehung des Koordinatensystems Invarıant 
ist, läßt sieh in (21) für Fy?+F,? F2+ F,? schreiben. Wir erhalten also aus (19), (20), (21) 
endlich 


7 | Fe: 
(1. ö.4* \\(dr+26G9)öFdzdy—\F,öFds—\| R 


H 


Il) Fjoands 0°) (232). 


IV 


S Ik di 
Wegen (1) ist das erste Glied Null. Zerlegen wir den Rand # in den Umfang U und 
in die einzelnen Gleitlinien, so verschwinden alle Glieder, die sich auf den Umfang beziehen, 
da dort Ö8F=0, ön=0. Das Randintegral zerfällt in Glieder, die sich aus den einzelnen 
Gleitlinien ergeben. Diese sind voneinander unabhängig, müssen daher einzeln verschwinden. 
Wir fassen eine dieser Gleichungen näher ins Auge. Sıe lautet 
Ft 2G9Flönlds=0 . . . ... (823), 


/ 


w 


\IMÖöF-+| 
3) Die Herleitung einer Formel, von der (22) einen Spezialfall darstellt, ist bei Ch. Jordan, Cours d’Analyse 

Mathematique, 3. ed. t. III, p. 550—555 zu finden. 
5" 
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wo die Integration doppelt, einmal vom Rand in das Innere des Querschnittes mit dem Zeiger 1 


und dann wieder zurück mit 


Inteeration um. so können beide Integrale 
Sn 

\I7„.,ö0F, +F.0ofF, 

da F und folelieh auclı F? auf beiden Seit 


eration Ist auf dıe Länge: 


dem Zeierer 2 auszuführen ıst. 


! der Gleitlinie auszudehnen. 


Kehrt man den Sinn der zweiten 


vereiniet werden. Man erhält so mit $n: on, 
Fi.’ FF. | 

wi , () | dd 5 v) (24 = 

en der Gleitlinie denselben Wert hat. Die Inte- 


alle 


Kntlanz der Gleitlinien besteht die Nebenbedingune F-—- ks=0, Denken wir uns 
(rößen dureh einen Parameter £ ausgedrückt, welcher für s=0 gleich Null ıst, so können 
wır J /, ya y”’dt=0 schreiben, wo die Striche dıe Derivatıion nach ?f bezeichnen. 
Die Variation dieser Gleichung ergibt 


() F - R- x 


Das erste Glied verschwindet an der unteren Grenze, da dort 92 =dy=U. 


Man erhält 


a oy)d ? 


Den Vektor mit 


den Komponenten d.x, dy können wir durch die Komponenten Ös, dn ausdrücken. Bezeichnet 
o das Krümmungsmaß der Gleitlinie, positiv, wenn der Krümmungsmittelpunkt auf der Seite 


Jr 


der positiven n liegt, so erhält man aus (25) 
N) F (k F)o N F', () 17 ], Ö N) NH da bal. 
wo die Borenlänee als Inteerationsvarlable mit o bezeichnet wurde. Für «dıe andere Seite eilt 
Öf, (k F.os Fr OM | ooando (265). 
u 
Diese Werte ın (24) einzesetzt erhält man 
| F', LP | z 
\ F. (f F\,.)ds . on kıf -Fn:)\vlolon do dAsSs=U r 2 (u) 
Ist die Gleitlinie, wie wir bisher angenommen haben, gerade, so folgt wegen 0 0, F,—/| 
ınd ı F' Ps. ER 
F\, Pi. (28). 
Im Falle, daß o==0, definieren wir die Hilfsgröße do od ndo und erhalten aus (27) 
| | d 
mit Rücksicht auf On —- 008 
() ‘IS 
ER, ’_... 
\ Ve kiF. F 2. ART (Is En (79) 
Io ds 
Durch partielle Integration folgt hieraus 
} F, ER ER, nei 
ATi \ - kiF..+/ \omds: (), 
0 } (ds ID z en 
Das erste Glied verschwindet, da an der unteren Grenze d®=U undan der oberen F„,,—+ F„. = U 
Ist Weiter folet 
d F, F', .) 
(U), 


(is Io 
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er ] An den Nullstellen von o kann mit o nicht dividiert werden. An diesen Stellen muß 
iten 


° ( \ .y r . . Fur . . 
jedoch auch N dc von derselben Größenordnung Null sein wie o, da dön beschränkt ist. Bei 
(IS .. 


der Herleitung des Fundamentallemas der Variationsreehnune nimmt man für do» Funktionen 
an. für welche dies keine Einschränkung bedeutet. so daß Gl]. (30) aueh dann eültie ist. wenn o 
24), eetrennte Nullstellen hat. 


Laut Gl. (30) sind eekrümmte Gleitflächen wohl mörlich. Entweder ist also die an sich 


nte- naheliegende Annahme gerader Gleitlinien falsch, oder ist die Anwendung des Castig- 
lianoschen Prinzipes nicht gestattet, da es sich hier nicht um einen Gleiehgewichtszustand 

alle handelt, sondern um die Ausbreitung einer Gleichgewiehtsstörung. Zur Entscheidung sind 

nen weitere Untersuchungen notwendig. 

nen. 6. Numerisches Beispiel. Wir betrachten den verdrehten Rundstab im Augenblick, wenn 
die Entwicklung der ersten Gleitfläche zum Stillstand gekommen ist. 

Den Halbmesser des Kreisquerschnittes können wir gleich eins annehmen und er 
halten dann 

9% (rı) 
(25). ie ai —r?) 
In den Gleitlinien ist s—= 1 - r. Aus (6) folgt 
. (11) 
7, ks a 11 —r)—kiil—r) . :. .». . u N 
Für den Kreis ist die Greensche Funktion in komplexer Schreibweise bekanntlich 
. ’ 
mıt | l— 2z 
I : log 3 a GE Gr DE (92), 
ınet an 2 —2 
veite : ; ns 
wo z=er-+iy, !2=S+in. Legen wir die x-Achse in die Gleitlinie, so haben wir aus (S) 
zur Bestimmung von v die Integralgleichung 
ba). 1 
51 — ar) kl —)=\r(d),— log | — BE... | : (83). 
gilt 2 | er E_r 
7 
6) Die Inteeration 1st bıs zum Werte A auszudehnen. für den via) (). 

Will man bei vorgegebenen Werten von W, (@, k den Wert von a berechnen, so kann 
man derart verfahren, daß man die Interralgleiehung (33) für verschiedene, willkürlich an- 
genommene Werte von a löst. Für jedes a erhält man eine Belegungsfunktion »„(r). Die 

(27) Lösung der Gleichung v„(a)=0 durch Interpolation liefert den gesuchten Wert von a. 
Da es sich hier nur um ein Beispiel handelt, schlagen wir einen anderen, indirekten, 
1. aber weniger mühsamen Weg ein. Wır wissen, daß » für die beiden Endpunkte der Gleit- 
linie verschwindet und setzen 
(28) ’ | e- 
Br r(E)=- Cl — EZ lE-—-ua) . ... a 
(27) mit einer Konstanten Ü. Nehmen wir a willkürlich an, so können wir 9, G, k so berechnen, 
daß die beiden Seiten von (33) je weniger voneinander abweichen. 
Für «a=02 wurden die Werte von 
’ z a | u 2 a 
(29, / \(] Aaul!- a) '’]oge - d = a (.3.)) 
29). EB _n 
berechnet und in die zweite Spalte von Tabelle 1 eingetragen. Mit diesen Werten wird 
(4 (e'] 
(] a”) (| x) 2 Zr + 186), 
2 k 2akM 
„ . r 4 . >. . n . (11) 
wo (” eine neue Konstante und M=='""loge ist. Wir bestimmen nun die Konstanten 1. 
0 u 
und €” mit Hilfe der Methode der kleinsten Quadrate so, daß (36) womöglich genau erfüllt 
’ e ö (1 1) e ‚ 
sel. Es ergibt sıch 0,94 und ( 3.1062. Wird noch k= 2000 ke/em?, G = SOON ke/em 


(30). 4 >]. 
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en anzenommen, so hat man 9 —= 0,0047, ferner Ü 16950. Mit diesen Werten wurden die Werte 
von u, F*, F ın Tabelle 1 berechnet. Wegen der Symmetrie ist r, 2 Diese Werte sind 
in die vorletzte Spalte eingetragen. Die letzte Spalte enthält 7=(r,?+ r,?)'!2. Der größte Wert 
von r beträgt 2420 kg/em” in der Mitte der Gleitlinie, 3700 kg/em? am Querschnittsumfang. 
An Hand von Tabelle 1 läßt sich unser quadratischer Ansatz für » prüfen. Zwischen 
x—=0.> und » 1.0 sollte der Unterschied zweier benachbarter F-Werte 200 betraeen. In der 
"abelle finden wir zwischen 2£=02 und 03 den Unterschied 177,9, also eine Abweichung 
von Ilv.H. Bei der rohen Annäherung von » ist die Übereinstimmung noch befriedigend. 
Die Kenntnis von » ermöglicht die Berechnung von «u, F*, F für beliebige Stellen. 
Kinige Werte wurden für die y-Achse mit Hilfe der Simpsonschen Regel berechnet und ın 
labelle 2 zusammengestellt. In der letzten Spalte sind auch die x-Komponenten von r be- 
rechnet. 
Tabelle 1 
E I u I F E E, T 
| 1.0 0.OOOU U 0,0 0,0 3700 
4 1.9 0.0010 62 3572 351.0 3320 
i 0.8 0.0022 13.7 676.8 663,1 2930 
En 0,7 0,0035 21,1 958,8 937.1 250 
H 0,6 0.0049 30.4 1203.2 172.8 2160 
r (),.) V,UOU066 11.0 1410.0 1369.0 1770 
0.4 0,0085 52,8 1579,2 1526,4 370 
0.3 0.0108 67.1 1710,8 1643,7 970 
0,2 0.0134 83,2 1804,8 1721.6 70 
0,1 V,.0167 103.7 1561.2 11797.9 150 
0,0 V0208 1292 180,0 170.8 300 
01 ),0262 162.5 1861.2 1698.4 SI 
0.2 0.0544 213.7 IS04.N 1591.1 OO () 2000 
0.3 0,0479 297.6 1710.8 1415.2 OO 590 2090 
(0,4 0.0597 370,9 1579.2 1?083 — 000 1020 IIH0 
0. 1,0064 112.5 1410,0 997,5 — 2000 1270 3370 
06 0,0663 111.9 1?203.> 7913 ZE 1360 2420 
ar, VODSS 365.3 OAS.S 593.5 200 1270 2370 
0.8 0,0443 275.2 676.8 101.6 2000 1020 2250 
(0,0 0,0230 148.5 391.2 208.7 200 590 2090 
1.0 ATI 0), 0,0 0 2000 0) 200 














l[abelle 2. 
F' F T} 
(dt) 128,0 ISSU.0 1751.1 () N, 
1] 125.1 1S61.2 1736.1 310 
0,2 114,8 IS04,8 1690.0 630 Fund —— 
0.3 100.5 1710, 1610.53 IS 
0,4 4,6 1579,2 144.6 1540 
(I. 68.3 1410.0 1341,7 1720 > l \ \ 
Z zZ 

(6) 325 1?03.> 1150.7 2100 T 
0.7 37.7 958,8 921,1 2490 

Q 23.9 676.8 652.9 DISS + 
100 11.4 357.2 345. 3260 [IE “ne 
1.0 RE EAL 0 3650 


Abb. 1 


\uf Grund der Daten von Tabelle 1 und 2 wurde in der Abbildung F mit Hilfe seiner 
Schichtenlinien roh dargestellt. Zur genaueren Darstellung müßte man von einer genaueren 
Bereehnune von » auszeehen und F aueh für Werte außerhalb der Koordinatenachsen be- 
rechnen, was keine prinzipielle Schwierigkeiten verursacht, aber einen beträchtlichen Mehr- 
aufwand an Rechenarbeit fordert 
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8. Verallgemeinerungen. Wir wollen unsere Ergebnisse verallgemeinern. Bei dieser Ver- 
allgemeinerung verlassen wir das Gebiet des exakt Bewiesenen und begeben uns auf das 
(rebiet der Vermutungen. Diese Vermutungen haben den Charakter einer Interpolation zwischen 
den Ergebnissen unserer Berechnungen für den speziellen Belastungsfall eines idealisierten 
Körpers und den experimentellen Beobachtungen. 


Das wichtigste Ergebnis unserer Betrachtungen ist die Erkenntnis, daß die Abnahme 
der Spannung bei zunehmender Formänderung zur Entstehung von Gleit- 
linien führt. Eine Abweichung von den Eigenschaften unseres idealisierten Körpers kann 
wohl Änderungen in der Spannungsverteilung gegenüber der bereehneten hervorrufen, die 
Entstehung der Gleitlinien erfolet aber genau in derselben Weise wie beim idealisierten 
Körper. Die Stelle mit der größten Schubspannuug befindet sich im rein elastischen Zustand 
in den meisten wenn nicht ın allen Belastungsfällen auf der Oberfläche des Körpers. 
Durch die Abnahme der Spannung in der Nähe dieser Stelle setzt die Bildung einer Gleit- 
läche ein, die sich infolge der Kerbwirkung ihres Randes plötzlich ausbreitet. Durch die 
hierdurch hervorgerufene Entlastung entstehen weitere Gleitflächen, die das Bestreben haben, 
sich äquidistant zu verteilen. Es hat also einen Sinn, von der Dichte der Gleitllächen zu 
sprechen, die jedoch nicht räumlich konstant sein muß. Die durch das Auftreten der Gleit- 
linien hervorgerufene Entlastung zeigt sich im Formänderungs-Belastungsdiagramm als eine 
Unregelmäßigkeit an der Fließgrenze. 


Die plastische Formänderung, in die Gleitflächen konzentriert, geht mit einer Verfestigung 
einher. Wir finden statt Gleitllächen Fließschichten. Mit fortschreitender Formänderung 
werden diese Fließschiehten immer dieker, da immer weitere benachbarte, noch unverfestigte 
Schichten ins Fließen geraten, bis sich der ganze Körper im plastischen Zustand befindet. 


Die Behauptung, daß die Abnahme der Spannung bei zunehmender Formänderung zur 
Entstehung von Gleitflächen führt, kann nicht ohne weiteres umgekehrt werden: man kann 
aus dem Auftreten von Gleitflächen nieht auf das Vorhandensein einer oberen und unteren 
Fließerenze schließen, da es wohl vorstellbar ist, daß auch andere Ursachen zur Bildung von 
Gleitflächen führen können. 


Auf den homogenen Spannungszustand können unsere Betrachtungen nicht unmittelbar 
angewandt werden. Auch können wir ihn praktisch nie ganz genau verwirklichen. Sind 
aber Abweichungen vom ideal homogenen Spannungszustand vorhanden, mögen sie noch 
so klein sein, so besteht keine Schwierigkeit mehr. Der homogene Spannungszustand er- 
scheint hier als eine Entartung des allgemeinen, der nur durch einen Grenzübergang erfaßt 
werden kann. Dies ist wohl der Grund. weshalb alle Bemühungen, die Lüders-Hart- 
mannschen Linien vom homogenen Spannungszustand ausgehend zu erklären, scheitern mußten. 


Eine Verallgemeinerung ganz anderer Art ist die Anwendung unserer Gedankengzänge 
auf Kristalle. Das Gleiten erfolgt bei Kristallen in Gleitebenen, die dureh niedrige Indizes 
ausgezeichnet sind. Die Gleitspannung ist beim Gleiten in einer solchen Gleitebene eines 
idealen Kristalles eine periodische Funktion des Gleitweges, deren Periode gleich der Elementar- 
distanz des Gitters in der Gleitrichtung ist. Infolge der Periodizität muß «die Spannung in 
bestimmten Abschnitten der Formänderung abnehmen, ja, sogar das Vorzeichen ändern. 


Die Umstände, die in isotropen Körpern zur Bildurg von Gleitflächen führen, sınd also 
vorhanden. Da es sich hier um Dimensionen der Gitterkonstante handelt, können die Gleit- 
ebenen auch von den inneren Oberflächen der Gitterfehlstellen, die laut den experimentellen 
Beobachtungen immer vorhanden sind, ausgehen. Es besteht aber gegenüber den isotropen 
Körpern ein grundlegender Unterschied: bei den Kristallen sind die Gleitflächen Ebenen, die 
nicht nur durch den Spannungszustand, sondern auch durch die Struktur des Krıstalls be- 
stimmt sind. Von den Gleitebenen sind jene wirksam, in denen die Schubspannung einen 
nur von der Orientierung der Gleitebene und der Gleitriehtung abhängigen, also vom Spannungs- 
zustand unabhängigen kritischen Wert überschreitet. Im allgemeinen sind die Gleitebenen 
Ebenen mit der kleinsten kritischen Schubspannung. Die Kleinheit der kritischen Schub- 
spannung in den Gleitebenen gibt aber keine Erklärung dafür, warum von den parallelen 
Ebenen mit derselben kritischen Schubspannung nur einzelne als Gleitebenen wirksam sind. 
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Gesellschaft für angewandte Mathematik und Mechanik. 


Hauptversammlung in Göttingen. 


Die diesjährige Hauptversammlung der Gesellschaft für angewandte Mathematik und 
Mechanik fand im Anschluß an die Tagung des „Fachausschusses für Strömungsforschung 
beim VDI* in der Zeit vom 20. bis zum 22. Oktober in Göttingen statt. Sie war später als 
sonst üblich zeleet worden, um den Mitgliedern der Gesellschaft, die den internationalen 
Mechanikkongereß in Cambridge (Mass.) besucht hatten, die Möglichkeit zu geben, auch an der 
Göttineer Tazrune teilzunehmen und dort über die Verhandlungen des internationalen Kon- 
eresses zu berichten. Eine Übersicht über die Reise der deutschen Delegation und den all- 
eemeinen Verlauf des Kongresses gab Herr Prandtl-Göttingen; über die einzelnen Ver- 
handlungsgebiete referierten die Herren Tollmien-Dresden (Hydrodynamik, insbesondere 
Turbulenz), Schliehting- Braunschweig (Aerodynamik) und Marguerre- Adlershof (Elastı- 
zität). Kurze Inhaltsangaben der übrigen in Göttinzen gehaltenen Vorträge sind, soweit sie 
von den Vortragenden zur Verfügung gestellt wurden, im folgenden, abgedruckt. Diese 
Berichte werden nebst dem im nächsten Jahrgang dieser Zeitschrift erscheinenden ausführ- 
lichen Vortraesberieht des Herrn Collatz über die Bereehnune von Näheruneswerten füı 
die Eigzenwerte den Mitgliedern der Gesellschaft zugehen. 


Am Freitag fand vor Beginn der Vorträge unter dem Vorsitz von Herrn Prandtl die 
ordentliche Hauptversammlung statt, über deren Verlauf der ausführliche Jahresbericht die 
Mitelieder unterriehten wird. 


Vorträge und Hauptversammlung waren in diesem Jahre so gut wie bisher beı keiner 
Taeunge besucht. Zu manchen Vorträgen hatten sich fast 150 Hörer eingefunden. Ebenso 
war die Beteiligung an den am Freitagnachmittag stattfindenden Besichtigungen der Aerody- 
namischen Versuchsanstalt und des Universitätsinstitutes für angewandte Mechanik überaus 
stark. Es spricht unter anderem für den in jeder Weise erfreulichen Verlauf der diesjährigen 
Tarunge, daß während der Tagung selbst mehr als 25 neue Miteliedermeldungen beim Ge- 
schäftsführer der Gesellschaft einzinzen. 905 


Vorträge aus dem Gebiet der Aero- und Hydrodynamik. 


Über das Messen turbulenter Längs- und Querschwankungen. 
Von H. Reichardt in Göttingen. 
\us dem Kaiser Wilhelm-Institut für Strömungsforschung.) 


I. Über die Bedeutung von Schwankungsmessungen. Kin tieferes Verständnis der tur- 
bulenten Strömungsbewegungen läßt sich aus den Untersuchungen der zeitlichen Mittelwerte 
der Strömungen allein nicht gewinnen. Ein Studium der Schwankungsvorgänge selbst ist 
unerläßlich. Die Turbulenzforschung ist daher in ein neues, erfolgversprechendes Stadium 
zetreten, seitdem man begonnen hat, Schwankungsmessungen auszuführen und auf diese 
\lessungen Theorien zu gründen, 


Die bisher durchgeführten Schwankungsuntersuchungen beziehen sich fast ausschließlich 
auf die sogenannte isotrope Turbulenz. Sie stellt einen Abklingungsvorgang dar und ist der 
einfachste Typ der Turbulenz, bei der man bloß die Längsschwankungen zu messen braucht. 
Das wesentliche Problem der turbulenten Austauschbewegung bzw. der turbulenten Schein- 
reibung steckt aber gerade in der nichtisotropen Turbulenz. Für deren Untersuchung sind 
\Messungen von Längs- und Querschwankungen erforderlich. 


2. Richtungssonden für Schwankungsmessungen. Zuverlässige Messungen quadratischer 
\littelwerte, wie man sie für die Schwankungsforsehung braucht, sind praktisch nur mit Hitz- 
drähten ausführbar. Für das gleichzeitige Messen von Längs- und Querschwankungen ist 
jede Hitzdraht-Richtungssonde verwendbar, die sich in hinreichend geringen Abmessungen 
herstellen läßt. 

RKiehtungssonden bestehen mindestens aus zwei eleichartigen Hitzdrähten, die winkel- 
empfindlich angeordnet sind. Die Differenz der Spannungen über diesen Drähten ist ein Maß 
für die Wıinkelabweichung bzw. für die Quergeschwindigkeit »’ der Strömung und die Änderung 
der Spannungssumme ist ein Maß für die Änderung w’ der Längsgeschwindigkeit. 
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Bei der Anordnung von Simmons und Bailev (Abb. la) beruht die Riehtungs- 
empfindlichkeit der Hitzdrähte auf ihrer: Schräganblasung. Bei den zweı parallelen Hitz- 
drähten von Burgers fin Abb. Ib im Sehnitt dargestellt) ist die Riehtungsempfindliehkeit 


durch wechselseitige Beeinflussung der Temperatur- und Gesehwindiekeitsfelder der Drähte 


bedingt. 
Bei der Dreidraht-Anordnung des Verfassers (Abb. le) werden die Geschwindigkeits- 
und Temperaturgradienten im Nachlauf des vorderen (dritten) Hitzdrahtes benutzt. Der 


vordere Draht kann zum Messen von w’, aber auch zum Messen von Temperatursehwankungen 
T’ dienen. Es ist daher möglich, außer der turbulenten Scheinreibung o w’r’ auch die tur- 
bulente Wärmeübertragung oc, T’»’ direkt aus den Schwankungen zu ermitteln (o == Luft- 
spezifische Wärme). 


dichte, Cy 


. \® 
& PR a 
LU —— -— Eure > . un. Et 3 
” . 
A371 a & 


Apbn.1. 


3. Allgemeine Beziehungen. Ganz unabhängig davon, in welcher Weise die Richtungs- 
empfindlichkeit erzielt wird, gelten für Richtungssonden, die für Schwankungsmessungen Ver- 
wendung finden, folgende allgemeine Betrachtungen. 

Die „wirksame Kühlgeschwindigkeit“ w') an einem Hitzdraht hängt ab vom Winkel « 
(Abb. 1) und von der Geschwindigkeit « der ungestörten Strömung. Die » (a)-Kurven eines 
Hitzdrahtes für die einzelnen Geschwindigkeiten sind symmetrisch zur Achse @=-0, auf der 
die Kleinstwerte von w liegen. Die Größtwerte von w sind im allgemeinen mit der unge- 
störten Geschwindigkeit « identisch. 

Für Schwankungsmessungen ist die Sonde so aufzubauen, daß die Hitzdrähte feste Winkel 
+ a bzw. a mit der Hauptströmungsriehtung x bilden. Dann werden beide Hitzdrähte mit 
der gleichen Geschwindigkeit w angeblasen, wenn der Strömungsvektor in der «-Richtung 
liegt und die Anblasgeschwindigkeit « örtlich konstant ist. 

Ändert die Strömung gegenüber der «-Achse ihre Richtung um den kleinen Winkel 4 


lu, lu, 


und ihre Stärke um Au und besteht ein Geschwindigekeitsgefälle r “in der y-Rich- 
tung senkrecht zu x, so hat man für die Änderungen der Kühlgeschwindigkeiten Aw, bzw. 


I »v, näherunesweise: 


. 


| a” „ar | (l) 
w,: po >, I ee re 
IT ya RT 
Ä Om RT) Ä ' 
WW, DD -P Be. :;; | > ’ : a - . i s r 2 1). 
ö Va’ Out . 


Die Spannung E über einem Hitzdraht ist eine Funktion von mw. Diese Funktion fm) 
hängt im einzelnen von der gewählten Schaltung ab. Allgemein ist zu sagen, daß f mit wachsen- 
dem w abnimmt bis auf den asymptotischen Wert E, des kalten Drahtes. 


0 


Für kleine Anderungen Aw kann man schreiben 


df 


Ben nr :; 2. a 
dm 


E=f(w)—+ 
bzw. für die in der Wheatstoneschen Brücke gemessene Spannungsdifferenz 
„dm ,omdu 


1f 
4 A ea Er 


. —_ dw\ da ou dy 
(19 Abstand der Hitzdrähte.) Führt man noch die Querschwankung "gu ein und 


bezeichnet 7 w/a. so hat man 


(»). 


d \ 


> — „r y | » 
e.= hiı +7,u+t 1 ‚pr 


I), In w sind definitionsegemäß auch die Einflüsse einer veränderten Lufttemperatur mit einbegriffen. 
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Hier bedeutet , ein kleiner Fehlwinkel, der beim Aufstellen der Sonde sehr leicht auftreten 
kann und 


d f ( / 
b==2 (6), 
“ 
d mo «A 
\ 
{ 7 
HM / 
( H ' J 
( r | - ıy (4) 
0) ( / 
) .) l 
= Q (1 0 dA 


Die Größe b ist ein Maß für die Empfindlichkeit der »’-Messung und e gibt den Einfluß des 
(eschwindigkeitsgradienten auf die »’-Messung wieder. 


Die Größe 5 ıst im alleemeinen geeschwindiekeitsabhäneie. Bei der Dreidrahtsonde läßt 
sich b allerdines im Bereiche kleiner Gesehwindiekeiten konstant machen. da bei den kleinen 
Revnoldsschen Zahlen des vorderen Drahtes (von etwa 0,01 mm 8) zwei Effekte gereneinander 


arbeiten: Während mit wachsendem ı® abnimmt. vertieft sieh allmählich die Geschwindie- 


dw 

Q r ; E j 

keitsdelle hinter dem vorderen Draht, wodurch \ ansteigt. Bei Steigerung der Geschwindig- 
7; 

keit weit über I m/s sehreitet aber die Dellenveränderune nieht mehr so stark fort und das 

” df. 
Abklıneen von überwiegt. 
dd I 


Da die turbulenten Gesehwindiekeitsschwankuneen nieht immer klein sind. erfordert die 
(veschwindiekeitsabhäneiekeit von b noch eine Korrektur. Man kann näherungsweise schreiben: 


a, | 
»=bhb(n)+ u WI an FE ET 
da 


wobei # die mittlere Geschwindiekeit und «’ die Längsschwankung bedeutet. Nach Einführung 
der dimensionslosen Zahl 


MH Ad, 


( 
H : &.- 5 1 u N 2 "Te (*)) 
h du 


hat man statt (5) unter Vernachlässieune der Korrekturgelieder zweiten Grades: 


un du 
e,‚b=v -y,u+x ur 0 STE BE re 
7 dy 


4. Bildung verschiedener Mittelwerte. Im allgemeinen mißt man nicht die Brücken- 
spannung selbst, sondern einen Strom, der von der Brückenspannung veranlaßt wird. Ist 
dieser Strom der Brückenspannung proportional, so sind auch die Strommittelwerte den ent- 
sprechenden Spannungsmittelwerten proportional. 


Für das lineare Spannungsmittel hat man 


u’ du 


eb QD,uUTrxX r( (il). 
’ 7 dy 
Diese Gleichung ıst bei der experimentellen Bestimmung von e zu verwenden. 
/ur Ermittlung des quadratischen Mittelwertes ist die Gleichung 
Ä ’ dw (We we 
e,=vV+tp,W+c | I Re? ara DE RZ 
dy 7 7 
zu quadrieren. Man erhält 
URN Ben N Wo? = 
ey» h° „’--2g u Vv I) ( Zu; ) L-2% i - (13). 
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. . 4 . . > V m . . ° 
Im e-Glied dieser Gleichung steckt die Annahme, daß u” __ wegen fehlender Korrelation 
V& 


verschwindet. 
Existiert für die Länesschwankuneen «’ die Gleichung 
ed " ‚ (14), 


so hat man für das gemischte Produkt w’»’, das der turbulenten Scheinreibung proportional 
ist, cemäß (12) 


’ ’ 


a 5 e dw: u’2r 
Z ! 35, » n x 
/ wu rypy,W°Tr59 ar (Im). 
ab 2 dy 7 


Die Glieder mit z in den Gl. (13) und (15) sınd klein und verschwinden außer für 2 = 0 


auch bei symmetrischer w"-»"-Verteilung. 











4 
2| 
4N vi nn / 
700 —— | Um=100 cm/s 
PR 
2! 
je ä 14 
N WERE! En l | LU - 2 
17] 2 r u. 70 72. 0m 0 £ Y 6 y 8 70 12 cm 
4 
Abh. ?. Abb, 3. 


5. Beispiel. Die Abb. 2 und 3 zeigen Y v’?- bzw. uw’ »’-Messungen mit der Dreidraht- 
sonde in einer ausgebildeten turbulenten Kanalströmung als Funktionen des Wandabstandes y. 
Die gerenseitigen Abstände der Hitzdrähte waren von der Größenordnung Ol mm. Der Ab- 
stand der Kanalwände betrug 24,4 em und die mittlere Maximalgeschwindigkeit U, 100 em/s. 


In dem vorliegenden Geschwindiekeitsbereich ist für die Dreidrahtsonde b = const und 
somit <=0. Die #-Glieder wären aber auch ohnehin fortgefallen. da sich die «-»’-Verteilung 
als symmetrisch herausgestellt hatte’). 


Die gestrichelte Kurve der Abb. 2 verbindet die y "’”?-Meßpunkte. Die ausgezogene 
Kurve gibt die unter Berücksichtigung des e-Gliedes verbesserten Y v"?-Werte wieder. 


Die durch den Nullpunkt gehende «’»’-Kurve der Abb. 3 entspricht der turbulenten 
Scheinreibung 


dü ® 
WW V" == Taon u (16). 
” “ ad Y 
un’ E 
Die gesamte Schubspannung Tges entspricht der durch die Punkte y = 12,2 und 10° 1. 25,8 


m 

gehenden Geraden. Die Lage dieser Geraden wurde aus dem Druckabfall bestimmt. Die mit 
der Dreidrahtsonde gewonnenen Meßpunkte liegen also im Mittel mit guter Genauigkeit auf 
der dureh das Druckgefälle und den Geschwindiekeitsgradienten zeforderten «’ »-Kurve. Ohne 
Berücksichtigung des e-Gliedes würden die ınittleren Meßwerte auf die darunter liegende 
schwach ausgezogene Kurve fallen. 397 


2), H. Reichardt: Messungen turbulenter Schwankungen. Naturwissenschaften 26 (1938), S. 404. 
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Frequenzanalyse turbulenter Schwankungen. 
Von H. Motzfeld in Göttingen. 


\us dem Kaiser Wilhelm-Institut für Strömungsforschung. 


Die der mittleren Geschwindiekeit einer turbulenten Strömung überlaeerten Geschwindie- 
keitsschwankungen lassen sich zusammengesetzt denken aus einer sehr großen Anzahl von 
harmonischen Schwingungen verschiedener Frequenz und verschiedener Amplitude. Für viele 
theoretische und auch experimentelle Fragen ist es nun von Bedeutung zu wissen, wie sich 
die gesamte Schwankungsenergie über den vorkommenden Frequenzbereich verteilt, d. h. man 
fragt nach dem Spektrum der turbulenten Schwankungsbewegung. 

Ist «.’ (f) die an einem festen Punkt innerhalb einer turbulenten Strömung der mittleren 
Geschwindiekeit überlagrerte Geschwindiekeitsschwankung, so ist der zeitliche Mittelwert der 
gesamten Schwankungsenergie «'?. Es sei nun flo) (»® =2rn Kreisfrequenz) die Verteilungs- 
funktion der Energie, für welche eilt 


WAC LE A500: GE re EG ee © 


ı) 
Dann stellt 
u? flo)Ao Beet: Ale. ©. r Beh, „1 EDER zutun. 0 SE | FZEEE BE Eee FE 


den Beitrag dar, den das Frequenzband von der Breite Io zur Gesamtenergie liefert, und 


\ 
v\flo)deo F(N) . ; (3) 


den Beitrag, den das Frequenzband zwischen O0 und X liefert. 

Ks ist nun gleichgültig, welehe der beiden Funktionen Fi(o) oder fi») experimentell 
ermittelt wird. Die Aufnahme der Geschwindigkeitsschwankungen geschieht am besten mit 
einer Hitzdrahtbrückenschaltung, wie sie des öfteren beschrieben worden ist. Die Ge- 
schwindigkeitssehwankungen « werden dureh eine solehe Schaltung in Spannungsschwankungen 
überführt. Durch Röhrenverstärker können diese Spannungsschwankungen verstärkt werden 
und «dureli geeignete Wahl des Arbeitspunktes der Röhren und durch Anwendung von Kompen- 
sationsschaltungen ist erreichbar, daß die Spannung am Ausgang des Verstärkers E'=k, «' 
wird (A, = const). 

legt man jetzt diese Spannung an die beiden Spulen eines Elektrodynamometers, so 
zeigt dieses einen Wert an, der proportional zu «’? ist. Schaltet man aber vor das Dynamo- 
meter eine Siebkette, die nur die Frequenzen & = N gleichmäßig hindurchläßt, so ist der 
von dem Dynamometer angezeigte Wert 

\ 
viflo)do=w"F(N). 
0 
Durch Änderung von N an der Siebkette läßt sich die Kurve «’? F(N) punktweise bestimmen. 
Dieses Verfahren wurde von Simmons und Salter') angewendet, um die Spektralkurven 
der turbulenten Schwankungen hinter einem Maschengitter in einem Windkanal zu ermitteln. 

Nach einem Vorschlag von H. Reichardt, von dem auch die Anregung zu den vor- 
liegenden Untersuchungen stammt, kann die Ermittelung von f(») mittels einer Reihe von 
mechanischen Schwingungssystemen von verschiedenen Eigenfrequenzen @, erfolgen, die durch 
einen zu E’ proportionalen Strom ./" elektromagnetisch erregt werden. Aus den Schwingungen 
dieser Systeme lassen sich die Werte von f(®) an den Stellen der Eigenfrequenzen w, be- 
rechnen. Eine kurze Betrachtung soll dies zeigen. 

Da bei den Versuchen, über die unten noch berichtet werden soll, als Schwingungs- 
systeme Drehschwinger verwendet werden, so wollen wir die Formeln auch für Drehschwinger 
anschreiben. Bezeichnen wir mit 


) das Trägheitsmoment des Schwingers, 
D die Direktionskraft des Torsionsfadens, 
R die Reibungsziffer der Dämpfung, 

\/ı(t) das erregende Moment, 





so Jautet die Gleiehunge der Drehscehwinzung 
BB + Dr HN . : >. 3 ua ar 


1), Simmons und Salter: An experimental determination of speetrum of turbulence. Proe. Roy. Soc. A 920, 
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Nehmen wir zunächst einmal an, daß die turbulenten Schwankungen sich aus S diskreten 
harmonischen Schwingungen zusammensetzen, deren Frequenzen Vielfache von Aw sind. 
So können wir schreiben 

N 


> Cy sın iv | of AR E’ k, - ; . (»)). 
] 


k,w=h, 
Ist nun die Meßanordnunge so eewählt. daß das errerende Moment 
\ 


Mtl)=k,J'=K N ec,sin(rdAwot-a,. . | (6) 


) 
| 


ist. dann erhalten wir als Lösung der Differentialeleichung 


N 
o=®Pe *!sıno,t- > 0, sin (v Ioati— a ur . MM 
| 
mit 
I er 
I) x\° Pr 
LI PRE- | | A ‘ıb) 
0 (-) ) 
und 
| Ke, ü 
( „ = i FE4) 


) s 
| a, (y Im?) T z(v | co»)? 


Den ersten Term der Lösung (7) brauchen wir nicht weiıter zu betrachten. Er klingt nach 
einiger Zeit ab. In dem zweiten Term stören die Phasenwinkel y. Um sie herauszuschaffen, 
erheben wir diesen Term ins Quadrat und bilden den zeitlichen Mittelwert. Es ergibt sich 
dann der einfache Ausdruck 


y : .) a . . . . . . . . R . ; ’ r (S) 


Von dem Quadrat der Gl]. (5) bilden wir ebenfalls den zeitlichen Mittelwert und erhalten 


- | , un 
x 1, FOR. a, DR 7er % -, dr 


Die Größen ec,” sind nun Teilbeträge der gesamten Schwankungsenergie. Nach Gl. (2) und (9) 


können wir setzen 
„=2u’flo)do. ... ei ER EI (10). 
Ist nun 5 eine sehr große Zahl (S>x) und A sehr klein (A&®>do), so läßt sich in Gl. (8) 


die Summe als Integral schreiben. Sie lautet dann, wenn wir noch die Werte für C,’ und e,? 
einführen, 


—.— a (m) 
Ti .. D ; u er re re 
| (-) w”| > x° “ 
( 


Um auf einfache Weise hieraus die Spektralkurve f(») zu ermitteln, bedenken wir, daß 
sich (11) auch schreiben läßt 


es 


do 


o=K,u’f(w,) D RN Er re (12), 
() a") 7 


« / 
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worin f(o,) der Wert der Spektralkurve an einer Zwischenstelle &, bedeutet. Da ferner gilt 


F 


do 7 
a, „10.c 2, ot2 Wave CE Tre 


\ | D 2 ! a 20, 
(-) ‘') | 


so ergibt sich, wenn wir nach f(o;) auflösen, 


fi.) N, 0. U N ee 2 Pe 

Die Stelle », ist uns nicht bekannt. Für kleine x ist aber die Resonanzbreite des 
Schwingers selır schmal. Wir machen also keinen großen Fehler, wenn wir als Näherunge 
setzen: ®- ” @,. Um die Größe dieses Fehlers zu ermitteln, wurden bei den Messungen 
(siehe unten) dem Schwinger verschiedene Dämpfungen gegeben. Es zeigte sich hierbei, daß 
für Frequenzen » — 4 Hz., der Fehler innerhalb der Mesgenauigkeit lag. Für tiefere Frequenzen 
dareren war ein Fehler festzustellen. Es ließ sich aber die durch Versuche ermittelte Ab- 
hängiekeit der Größe 9°” von # durch eine Potenz darstellen. Wir erhalten somit für sehr 


tiefe Frequenzen die empirische Formel 
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Abb. 1. 


Die Apparatur, mit der nach diesem Verfahren Frequenzanalysen durchgeführt wurden, ist 
in der Abb. I schematisch dargestellt. Es wurden als Schwingungssysteme Drehschwinger 
verwendet. An einem dünnen Torsionsfaden ist als Achse des Schwingers ein steifer Draht 
aufeehänet, der dureh ein Glasröhrehen in der -Mitte geteilt ist, so daß zwischen dem oberen 
und dem unteren Stück keine elektrisch leitende Verbindung besteht. An dem unteren Ende 
trägt die Achse einen Schwingkörper, an dem oberen einen Spiegel und in der Mitte eine 
Spule. Ein Ende dieser Spule ist mit dem oberen Stück der Schwingerachse, das andere 
mit dem unteren verlötet. Die Stromzuführung geschieht durch den Torsionsfaden und durch 
einen dünnen Platindraht, der in ein Gefäß mit Quecksilber taucht. In dem Gefäß befindet 
sich über dem, Quecksilber eine Dämpfungsflüssigkeit. In diese hinein taucht der Schwing- 
körper. Bei Kigenfrequenzen über 20 Hz. zeigte es sich, daß es vorteilhafter ist, dem Dreh- 
schwingeer einen etwas anderen Aufbau zu geben (Abb. 1b). Da das Trägheitsmoment der 
Spule und des Spiegels groß genug ist, erübrigt sich die Anbringung eines zusätzlichen Träg- 
heitsmomentes. Dagegen ist das Gewicht der Spule zu gering, um den Torsionsfaden zu 
straffen. Es befindet sich daher an dem unteren Ende der Schwingerachse ein zweiter Tor- 
sionsfaden, an dem ein Gewieht hängt, das beide Fäden strafft. Durch dieses Gewicht 
wird weiter erreicht, daß der Drehschwinger stets senkrecht hängt, wodurch überlagerte 
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Schwingungen, hervorgerufen durch das Schwerefeld, vermieden werden. Das Gewicht taucht 
in ein zähes Öl und bleibt bei den schnellen Oszillationen des Schwingers praktisch in Ruhe. 

In der zur Verfügung stehenden Apparatur befinden sich 6 Schwinger der ersten und 
t Schwinger der zweiten Art. Die Eigenfrequenzen. liegen zwischen 0,2 und 43 Hz. Die 
Schwinger sind so aufgehängt, daß sich die Spulen zwischen den Polen N und S eines 
Elektromagneten befinden. Für jeden Schwinger ist eine Beleuchtungsvorrichtung vorhanden. 
Der von dieser Beleuchtungsvorrichtung ausgehende Lichtstrahl wird von dem Spiegel reflek- 
tiert und fällt auf einen Filmstreifen. Letzterer kann mittels eines Motors mit verschiedenen 
Geschwindigkeiten vorbeigezogen werden. Die ganze Apparatur ist in einem lichtdiehten 
Kasten untergebracht. 

An die Spulen der Schwinger wird nun die Ausgangsspannung E’ des Verstärkers gelegt. 
Durch die jetzt auf die Spulen wirkenden Drehmomente werden die Schwinger angeregt und 
ihre Schwingungen werden auf den Filmstreifen aufgezeichnet. Die Versuche haben gezeigt, 
daß zur Erlangung guter Mittelwerte Beobachtungszeiten von 12 bis 15 Minuten erforderlich sind. 

Wie aus der Formel (15) hervorgeht, sind für die die Berechnung von f(w,) verschiedene 
Größen erforderlich. Die Konstante A, wird durch Eichung für jeden einzelnen Schwinger 
ermittelt. «’° wird mit dem Dynamometer gemessen. Zur Bestimmung der Eigenfrequenzen 
werden während der Messungen auf den Filmstreifen in gewissen Zeitabständen Lichtmarken 
gegeben. Das logarithmische Dekrement z ergibt sich aus der Abklingkurve der Schwinger. 
Der quadratische Mittelwert der Amplitude 9? wird aus den aufgezeichneten Schwingungs- 
kurven mittels eines Quadratplanimeters bestimmt. Zur Ermittelung der Abhängigkeit des 
Wertes 92 von x ist es notwendig, verschiedene Meßreihen aufzunehmen mit geänderter 
Dämpfung der Schwinger. Aus den zu einem Schwinger gehörenden Werten ist durch 
logarithmische Auftragung der Exponent p und der Faktor A, für diesen Schwinger zu finden. 
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In der Abb. 2 sind nun die Ergebnisse der Messungen wiedergegeben, die der Verfasser 
in der turbulenten Strömung eines rechteckigen Kanals von 24 em Höhe und 100 em Breite 
bei einer Mittengeschwindigkeit von 100 em/s nach diesem Verfahren durchgeführt hat. Es 
zeigte sich dabei, daß innerhalb der Meßgenauigkeit die Spektralkurven für die verschiedenen 
Wandabstände zusammenfallen. 

In der Abb. 3 werden nun diese Ergebnisse mit denen von Simmons und Salter ver- 

fin) . n W; Sa a 1722 
glichen. Als Ordinate ist hier aufgetragen L, über der Abszisse Mn Für die Turbulenz 
hinter einem Gitter ist U die mittlere Geschwindiekeit im Windkanal und Z die Maschen- 
weite des Gitters. Für die Wandturbulenz ist U die mittlere Geschwindiekeit im Kanal und 
I, die Kanalhöhe. Die ausgezogene Kurve gibt die Ergebnisse von Simmons und Salter 
wieder. Man sieht, daß sie mit den unseren gut vergleichbar sind. S3S 
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Drei Filme über konforme Abbildung. 
Von Helmut Heinrich in Breslau. 


Aufgabe eines mathematischen Unterrichtstfilms. der als Anschauunesmittel bei mathe- 
matischen Vorlesungen verwendet werden soll. kann es m. E. nur sein, eeometrische Be- 


Ä 


ziehungeen durch Bewegungesvorgänge darzustellen. So ıst der Grundgedanke, dem 
die voreeführten drei Filme über konforme Abbildung ıhre Entstehung verdanken. der. daß 
bei jeder Abbildung zwischen zwei Ebenen der geometrische Zusammenhang zwischen den 
Punkten der einen Ebene und ihren Bildpunkten in der anderen Ebene besonders deutlich 
wird. wenn man zueeordnete Bewerungen der beiden Punkte darstellt, d.h. wenn man einen 
Punkt der ersten Ebene eine bestimmte Bahn durchlaufen und den zugehörigen Bildpunkt 
die dadurch zwanesläufie gerebene Bewegung gleichzeitig ausführen läßt. 


Es ist einleuechtend, daß derartize Filme Triekfilme sein müssen. Das bisher übliche 
außerordentlich zeitraubende und darum auch sehr kostspielige Verfahren der Herstellung 
soleher Filme, bei dem für jede Bewegungsphase eine besondere Zeichnung notwendig war 
und die aufeinanderfolgenden Zeichnungen genauestens Justiert werden mußten, mag in erster 
Linie schuld daran sein, daß die Zahl der mathematischen Unterrichtsfilme verschwindend 
eerine ist. Um auf diesem Gebiete schneller vorwärts zu kommen, mußte ein anderes wesent- 
lich einfacheres Aufnahmeverfahren entwickelt werden. Es ist in der Reichsstelle für den 
Unterriehtsfilm (R. f.d. U.) in Berlin zum ersten Male im September 1937 bei der Herstellung 
des ersten der drei vorgeführten Filme erprobt worden. Es ist außerordentlich einfach und 
kann kurz folgendermaßen umrissen werden. 


Der zu verfilmende Bewegungsvorgang wird auf geeignetem Pauspapier mit Bleistift 
derart vorgezeichnet, daß die Kurven, deren Entstehung gezeigt werden soll, mit der not- 
wendigen Einteilung in kurze Schritte versehen werden. Diese Vorzeichnung wird auf eine 
von unten stark beleuchtete waagerechte Glasplatte gebracht, über der sich der Aufnalıme- 
apparat befindet. Die einzelnen Felder, in die die vorgezeichnete Kurve durch die Schritt- 
einteilunz zerfällt, werden nunmehr nacheinander mit schwarzer Tusche von Hand ausgefüllt, 
und nach jedem Schritt wird eine Aufnahme gemacht. Die Belichtung und Entwicklung des 
Films kann so vorgenommen werden, daß im fertigen Film von der Vorzeiechnung nichts zu 


sehen Ist. 


Diese Art des Triekverfahrens hat ersichtlich ihre Grenzen. Innerhalb dieser Grenzen 
aber ist es, was Einfachheit und Aufwand an Zeit und Kosten angeht, nicht zu übertreffen, 
und es führt sehr schnell zum Ziel. Die Vorzeichnung erfordert in Anbetracht der Ansprüche, 
die an Sauberkeit und Genauigkeit vor allem auch der Schritteinteilung eestellt werden, 
selbstverständlich ihre Zeit. die Aufnahme selbst aber ist bei einiger Übung der daran 
Arbeitenden in wenieen Taren erledigt. Da es sich bei den vorgeführten dreı Filmen um 
erste Versuche gehandelt hat, sind verhältnismäßig einfache Gegenstände aus dem Gebiete 
der konformen Abbildung für die Verfilmung ausgesucht worden. Der erste Film zeigt die 
unter der Bezeichnung Kreisverwandtschaft bekannten Beziehungen zwischen den Geraden 
und Kreisen der z- und der m-Ebene, wie sie bei der konformen Abbildung durch die Funktion 


—=-— auftreten. Der zweite Film stellt für die Funktion »=z-+-— die Abbildung der 


Nullpunktsgeraden und -kreise auf konfokale Ellıpsen und Hyperbeln dar. Der dritte Film 
schließlieh veranschaulicht für dieselbe Abbildungsfunktion die bekannte Konstruktion von 
Joukowskv-Profilen aus ihren Bildkreisen durch einen Bewegungsvorgang und zeigt noch 
die Zuordnung der zugehörigen Stromlinien bei zirkulationsfreier Strömung um Bildkreis 
und Profil. 

Kinzelheiten über den Aufbau der Filme, auf die hier nicht eingegangen werden kann, 
können leicht aus den Veröffentlichungen der R.f.d. U., die für jeden einzelnen Film heraus- 
gegeben werden, entnommen werden. 02 
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Zur Berechnung der turbulenten Ablösung. 


Von F. Sehultz-Grunow ın Göttineen. 


Um zur näherungsweisen Berechnung der turbulenten Ablösung an der ebenen Wand 
einen Beitraege zu hefern. wird an bekanntem Versuchsmaterial!) über die turbulente verzögerte 
Grenzschieht die Frare untersucht, ob in der Grenzschiehtbeweruneseleichung 


7; O7 OU Iodr 


Ü x OÖ „ 0x 00 „ 


(u, v, x, y Geschwindigkeitskomponenten und Koordinaten parallel und senkrecht zur Wand; 
o Dichte; r Schubspannung; U Potentialgeschwindigkeit am Rande der Grenzschicht) die unter- 
striehenen Glieder in einem Streifen ın der Grenzschicht, der durch gewisse Linien ein- 
gegrenzt ist, von gleicher Größe sind. Würde das Experiment die Frage bejahen, so ergäbe 
sich eine wesentliche Vereinfachung der Rechnung, indem dann in dem betreffenden Streifen 
die Bewegrungsgleichung 

O7; OU 


u U 
Tr 0x 


lautete, welche sich sofort integrieren ließe zu 


But et EEE ie 3 


p- l 2 l 


Index 2 beziehe sich auf den Querschnitt, in dem die Ablösung einsetzt, Index 1 bezeichne 
einen Querschnitt in der turbulenten Verzögerungsstrecke. Die Ablösung sei dadurch gekenn- 
zeichnet, daß «,—=0 ist, so daß dann aus (1) sich 


Bu io bee ana rs MR 


ergibt. Diese Kennzeichnung der Ablösung verlangt es, die Beziehung (2) in einem Streifen 
zu bestätigen, der sich im Verlaufe der Verzögerung nicht zu weit von der Wand entfernt. 
da erfahrungsgemäß die Ablösung ın unmittelbarer Wandnähe einsetzt. Deshalb wurde der 
Streifen durch Linien s —= konst. (d = Grenzschichtdicke) eingegrenzt und nicht dureh Linien, 
4 
deren Wandabstände der Verdrängungsdicke oder Impulsdieke proportional sind. Letztere 
Länzen hätten den Vorteil, sich wesentlich genauer als die Grenzschichtdieke aus dem Versuch 
ermitteln zu lassen. 
Die Gültigkeit der Beziehung (2) wurde am Versuch in der Weise geprüft, daß jeweils 
y r a - er 
auf Linien 5; 0,1: 0,15; 0,2 in Jedem Meßquerschnitt die Geschwindigkeit «, abgelesen, mit 
( 
der Potentialgeschwindigkeit U, des gleichen Querschnittes in (2) eingeführt und U, daraus 
berechnet wurde, das ist Jene Potentialgeschwindigkeit, bei der die Ablösung einsetzen soll. 
. ö > Ki u. Fe 2: ’ - ; 
Es ergab sich, daß dıe für s —(),15 erhaltenen U,-Werte mit Ausnahme des Beginns der Ver- 
( “ Ä 
zögerungsstrecke tatsächlich in guter Näherung mit der Potentialgeschwindiekeit des der 
Ablösung am nächsten gelegenen Meßquerschnittes übereinstimmen und daß die in den übrigen 
Abständen erhaltenen Werte eine Übereinstimmung zeigen, daß sich (2) auch in der Umgebung 


7 RR 27 2 rn ur 2: 
von -0,15 in guter Näherung bestätigt, womit sich die unvermeidbare Unrzenauiekeit in 


der Ablesung der Grenzschichtdieke als wenig einflußreich erweist. 

Der Vorteil bei der Anwendung der Beziehüng (2) zur Vorausberechnung der Ablösung 
liegt darin, daß nur der Verlauf der Pontentialgeschwindiekeit am Rande der Grenzschicht 
bekannt zu sein braucht. Ein Nachteil besteht in der wıllkürlichen Annahme von «,, da das 
Geschwindigkeitsprofil im Querschnitt 1 nicht bekannt ist. Am besten ist es, im Beginn der 
Verzögerungsstrecke nach Gruseh witz'), im Rest nach dem hier entwickelten, wesentlich 
einfacheren Verfahren zu rechnen, da das erste Verfahren im Verlauf der Verzögerungsstrecke 
gute Übereinstimmung mit dem Versuch ergeben hat, während sich das zweite Verfahren für 
den Rest der Verzögerungsstrecke als brauchbar erweist. Die Auswertung hat ergeben, daß 
man mit dem ersten Verfahren höchstens bis zu dem Wert 0,65 des in diesem Verfahren 
auftretenden Formparameters zu rechnen braucht. Der diesem Wert entsprechende, in (2) 
einzuführende Wert von u, läßt sich aus den Versuchen ermitteln, da die Form der Ge- 
schwindigkeitsprofile in der verzögerten Grenzschicht in guter Näherung nur von diesem 
Formparameter abhängt. S99 

1) E. Grusehwitz: Die turbulente Reibungsschiceht in ebener Strömung bei Druckabfall und Druckanstieg. 
Ing.-Arch., Band II (1931), S. 321. 
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Zur Berechnung von Kreiselrädern. 
Von A. Betz in Göttingen und I. Flügge-Lotz in Berlin. 


Die Berechnung der Schaufelform von Kreiselrädern erfolgt in vielen Fällen mit Hilfe 
der Eulersechen Theorie unter der Annahme, daß die Strömung nicht nur auf der Oberfläche 
der Schaufeln der Schaufelform folgt, sondern auch zwischen den Schaufeln annähernd den 
eleiehhen Verlauf hat. Bei schnellaufenden Rädern mit geringer Schaufelzahl, bei denen diese 
Voraussetzung nieht mehr zutrifft, hat man Verfahren, die in der Flugtechnik für den einzelnen 
"lügel entwickelt sind, zur Berechnung herangezogen und den Einfluß der Nachbarflügel als 
Störune hinzugenommen. Ein wesentliches Hilfsmittel dieser Methode ist die konforme 
Abbildung. Diese Berechnungsart ist um so bequemer, je weiter der Abstand der Schaufeln 
ist. Das vorlieeende Verfahren soll die Lücke zwischen den bestehenden Methoden ausfüllen. 
Man denkt sich die einzelnen zylindrischen Schaufeln durch tragende Wirbelflächen ersetzt 
und bereehnet die Strömungsfelder bei Rädern mit konstanter Radbreite und dünnen Schaufeln 
mit Hilfe periodischer Wirbelfelder, bei Rädern mit abnehmender Radbreite mit Hilfe 
periodischer Wirbel- und Quellfelder. Die endliche Dicke der Schaufeln kann ebenfalls durch 
hinzutretende @Quellfelder berücksichtigt werden. Es hat sıch gezeigt, daß die Abnalıme der 
Radbreite für Räder mit mehr als zwei Schaufeln nur Einfluß hat auf die mittlere radiale 
Durchfllußgeschwindigkeit, so daß man praktisch mit Wirbelfeldern allein auskommt. Der 
Aufbau der Rechnung gestattet es, den Einfluß der endlichen Schaufelzahl auf die Wirkungs- 
weise des Rades abzuschätzen. 6 





Messungen an Kreistragflächen und Vergleich mit der Theorie der 
tragenden Fläche. 


Von M. Hansen in Göttingen. 


Vor 2 Jahren berichtete Kinner?) auf der Tagung der GAMM in Dresden über seine 
theoretischen Rechnungen zur Aerodynamik der Kreistragfläche. Inzwischen wurden zur 
Nachprüfung der Theorie in der Aerodynamischen Versuchsanstalt, Göttingen, Messungen an 
ebenen und gzewölbten Kreistragflächen vorgenommen, die in der Arbeit von Kinner theore- 
tisch behandelt waren. Im folgenden werden die gefundenen Meßergebnisse mit der Theorie 
verglichen. 

Aus der Arbeit von Kinner sei hier zum besseren Verständnis der Meßergebnisse nur 
ein Punkt näher behandelt, nämlich die Saugkraft. Sie entsteht bekanntlich beim Umströmen 
des Vorderrandes, also bei nicht stoßfreiem Eintritt, und äußert sich in einer Erniedrigung des 
Widerstandes. In der Theorie der tragenden Linie ist die Saugkraft im Widerstand enthalten. 
Der so errechnete Widerstand sei mit Linienwiderstand W7, bezeichnet. In seiner Arbeit 
berechnet Kinner für die Kreistragfläche auch den Widerstand olıne die Saugkraft. Dieser 
Widerstand heiße Flächenwiderstand Wr. Der Unterschied der beiden Widerstände ergibt 
also die theoretische Sauekraft bei Umströmung der Vorderkante. Man bekommt somit 
iheoretisch zwei Polaren, indem man den Auftriebsbeiwert c„ einmal zusammen mit dem 
Beiwert des induzierten Widerstandes cp und das andere Mal mit ec»; zusammen aufträgt. 











Abt 
Kine ausführliche Veröffentliehung erfolgt demnächst im Ingenieur-Archiv. 
®, Kinner, W.: Die kreisförmige Tragfläche auf potentialtheoretischer Grundlage. Ingenieur-Archiv. Bd.S, 


H. ı (1837). Vel. ZAMM, Bd. 16 (1936), S. 349 bis 352. 





Band 18, Heft 6 



















































4 Desamber 1928 Vorträge der Hauptversammlung in Göttingen 369 

Bei der wirklichen Tragfläche ist eine Umströmung der Vorderkante im Sinne der 
Theorie der tragenden Linie, d. h. unter vollständiger Erhaltung der Sauekraft, nieht zu 
erwarten. Man weiß andererseits aus einer genaueren Untersuchung der Strömunesverhältnisse 

fe an der Vorderkante bei verschiedener Abrundune von Wolff?) daß selbst bei der feinsten 
E möglichen Zuschärfung immer noch eine geringe Sauekraft vorhanden ist. Um einen gewissen 
Ren Einblick in die Verhältnisse an der Anströmkante zu bekommen, wurden bei den Messungen 
a an der ebenen Kreistragfläche diese Kanten verschieden ausgebildet. 
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Abb. 3. 

Es wurden 3 ebene Flächen mit den auf Abb. 2 und 3 skizzierten Vorderrandausbildungen 
und 4 gewölbte untersucht. Von den gewölbten waren 2 Kugelkalotten mit den Pfeilhöhen 
h/D=-0,025 und 0,05, die restlichen beiden S-Flächen mit den Werten a/D = 0,067 bzw. 0,134. 
Die S-Flächen (Abb. 1) wurden nach einer von Kinner angegebenen räumlichen Funktion 
gerechnet. Diese Funktion ergab sich aus 2 Bedingungen: Die Tragfläche sollte druckpunktfest 
sein wie die ebene Kreisfläche. Sıe sollte dabei aber eine geringere theoretische Sauekraft 
besitzen als die ebene Kreistragfläche. 

Für die Kraftmessungen waren die ebenen Kreisscheiben aus Eisenblech und die ee- 
wölbten zwecks leichterer Herstellung aus Leichtmetall von 4 mm Dicke und 400 mm © 
angefertigt. Die Messungen selbst wurden im Kanal IV der Aerodynamischen Versuchsanstalt 
bei einer Windgeschwindigkeit von 4O m/sece durchgeführt. Die Kräfte wurden in der üblichen 

de 8, 3) Wolff, E.: Einfluß der Abrundung scharfer Eintrittskanten auf den Widerstand von Flügeln. Ingenieur- 


Archiv, Bd. IV, 193, S. 521. 
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Weise ermittelt. Für den Vergleich mit der Theorie wurden nur die Kräfte bei kleinen 
\nstellwinkeln herangezogen. 

\us den Ergebnissen werden im folgenden drei Abbildungen herausgegriffen. Abb. 2 
zeiet die Auftriebsbeiwerte sämtlicher untersuchter Kreistraeflächen aufeetraegen über den 
\nstellwinkel a. Die theoretischen Kurven sind zu den Meßpunkten eingezeichnet worden. 
s ergibt sieh für die ebenen Kreistragflächen befriedigende Übereinstimmung von Theorie 
und Messung bis zu etwa @a=-3" bei den ebenen Flächen. Dann tritt rasch ansteigend mit « 
eine Abweichung ein. Der gemessene Auftrieb ist erößer als die Theorie ergibt. Bei den 
ewölbten Kreisflächen ist die Übereinstimmung mit der Theorie befriedigend etwa in der 
(ezend, wo die Saurkraft der Theorie nach Null ıst. Die zugehörigen Winkel sind bei den 
kalotten 2,45" bzw. 49°, bei den S-Flächen 3,855 ° bzw. 7,63 °. Oberhalb und unterhalb dieser 
Stelle beeinnt bald eine mit a@ stark zunehmende Abweichung wie bei den ebenen Flächen. 
Bei den S-Flächen mußte eine Vergrößerung des gemessenen Anstellwinkels «a um 0,5° bzw. 
0,5" vorgenommen werden, um Messung und Theorie einigermaßen zur Deckung zu bringen. 
Die Erklärung für diese Abweichung kann gesucht werden in der Reibungsgrenzschicht, 
welche die Wirkung der S-Krümmung abschwächt. Ein ähnlicher, wenn auch geringerer 
Kinfluß zeigt sich bei den Kalotten. Für diese wurden Jedoch die Meßpunkte nicht korrigiert. 

Abb. 3 gibt die theoretischen Polaren für die ebenen Kreistragflächen mit den zugehörigen 
Meßpunkten. Zwecks besseren Vergleichs mit der Theorie wurde der Widerstand cw, bei 
ce, 0, der die Oberflächenreibung darstellt, abgezogen. 

Die Meßerzebnisse der ebenen Kreistragflächen mit verschiedener Ausbildung der An- 
strömkante zeieten deutliche Unterschiede. Die größte Einbuße an Sauekraft hat erwartungs- 
eemäß die Scheibe mit zugeeschärftem Rande. Etwas mehr Sauekraft hat die Scheibe, beı 
welcher der Rand abgerundet war. Praktisch vollständig erhalten hat sich die Saugkraft beı 
der vorne mit Wulst versehenen ebenen Kreisfläche. 
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Abb. 4. 


Die Abb. 4 zeiet die Polaren einer gewölbten Fläche und zwar der Fläche mit starkem 
S-Schlae. Der Wert der Oberflächenreibung wurde hierbei ermittelt durch Verschieben der 
heillen theoretischen Polaren unter Beibehaltung der Abszissenachse. bis die Meßpunkte 
zwischen ep und Cw7 lagen. Der auf diese Weise gefundene Wert von c.„, stimmte praktisch 
mit dem e,., der ebenen Fläche von gleicher Abrundung der Vorderkante überein. Infolge 
der gewählten Abrundung ist zu erwarten, daß ein Teil der Saugkraft verlorengeht, was 
dureh die Messung bestätigt wird. 

An einer ebenen Fläche wurde außerdem eine Messung der Druckverteilung vorgenommen. 
/u diesem Zwecke wurde eine ebene Fläche von 4 mm Dicke mit abeerundeter Kante her- 
eestellt. Die mit dieser Platte zemessenen Drücke wurden mit den theoretischen Druck- 


verteilungen bei “Anstellwinkel 4° bzw. 8" verglichen. Bis 4° stimmen die Mittelschnitte 
verhältnismäßie gut mit der Theorie überein. Am seitlichen Rande zeigen sich stärkere 


Abweichungen an der Stelle, wo der Wirbelzopf sich bildet. Bei S°’ erkennt man bereits in 
der Nähe des Vorderrandes stärkere Abweichungen. Es zeichnet sich hier deutlich die 
Stelle des Wirbels ab, welcher sich bei der Umströmung der Vorderkante auf der Saugseite 
bildet. Im übrigen zeigen sich am seitlichen Rande die gleichen Abweichungen wie bei 4° 
Anstellwinkel. YOU 
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Zur Theorie der Leitapparate für Propeller’). 
Von A. Betz VDI in Göttingen. 

In dem Strahl hinter einem Schraubenpropeller sind zusätzliche Geschwindigkeiten 
enthalten, deren Energie einen Verlust bedeutet. Die Axialkomponente dieser zusätzlichen 
Geschwindigkeiten hängt nach dem Impulssatz mit dem Schub des Propellers zusammen, läßt 
sich also nieht vermeiden. Dagegen läßt sich die Energie, welche in den mit dem Drehmomen! 
zusammenhängenden Tangentialgeschwindigkeiten, dem Scehraubendrall, enthalten ist, durch 
flügelartige Flächen hinter dem Propeller, sogenannte Leitapparate, zurückgewinnen. Sieht 
man von den Verlusten, die dureh den Widerstand des Leitapparates entstehen, ab, so läßt 
sich die Drallenergie durch einen Leitapparat mit unendlich vielen Sehaufeln vollständig 
wiedergewinnen, bei endlicher Schaufelzahl dagegen nur zum Teil. Man kann sich nun die 
Frage vorlegen, wieviel Energie man mit einer bestimmten Schaufelzahl wiedergewinnen kann. 
und wie man die Zirkulation um die Schaufeln verteilen muß, um dieses Optimum zu erreichen. 
Diese Frage läßt sich mit den auch sonst in der Tragflügel- und Propellertheorie üblichen 
Verfahren behandeln. 

Zunächst wurde die günstigste Zirkulationsverteilunge bei endlieher Schaufelzahl des 
leitapparates untersucht, wobei für die Schaufeln in radialer Richtung keine Begrenzung 
vorgeschrieben war. Auf Grund der gefundenen Zirkulationsverteilung ergibt sieh dann der 
Energie- bzw. Schubgewinn durch den Leitapparat für verschiedene Schaufelzahlen. 

Bei dieser Rechnung zeigte sich, daß die außerhalb des Strahles liegenden Teile des 
Leitapparates noch verhältnismäßig viel Zirkulation haben mußten, um den Optimuns- 
bedingungen zu genügen. Um zu erkennen, ob diese äußeren Teile auch für den Schubzewinn 
wesentlich sind, wurden die Rechnungen dann auch für Leitapparate mit begrenzter Schaufel- 
länge durchgeführt. Es ergab sich, daß bei geringer Schaufelzahl die außerhalb des Strahles 
liegenden Schaufelteile auch noch merklich zum Schubgewinn beitragen. 

Dieses Ergebnis ist aber für die praktische Verwendung irreführend, da der Widerstand 
der Schaufeln vernachlässigt ist. Wenn man diesen Widerstand berücksichtigt, der bei den 
äußeren Schaufelteilen verhältnismäßig groß im Vergleich mit dem theoretischen Schubgewinn ist, 
so ergeben sich je nach dem Drall des Schraubenstrahles und je nach der Gleitzahl der 
Leitschaufeln Schubgewinne, welche bei einer bestimmten Schaufellänge ein Maximum haben 

Die Nabe des Leitapparates kann einen anderen Durchmesser als die Nabe des Propellers 
haben. Auf das kegelige Übergangseebiet zwischen Propellernabe und Leitapparatnabe wirken 
infolge der vom Drall herrührenden Zentrifugalkräfte Unterdrücke, die eine Kraftkomponente 
in axialer Richtung ergeben. Ist die Leitapparatnabe größer als die Propellernabe, so ergibt 
sich eine Kraft in Richtung des Schubes, also ein Gewinn, welcher aus der Drallenergie 
entnommen wird. Umgekehrt entsteht, wenn die Leitapparatnabe kleiner als die Propellernabe 
ist, ein zusätzlicher Widerstand, der in einer Erhöhung der Drallenergie zum Ausdruck kommt. 
Es ist deshalb nicht günstig die Leitapparatnabe kleiner als die Propellernabe zu machen. 
Wenn sich hinter dem Propeller noch Körper mit größerem Durchmesser befinden (Flugzeug- 
rumpf, Motorgondeln), so ist es zweckmäßig, die Leitschaufeln an der dieksten Stelle dieses 
Körpers anzubringen. so 


Vorträge aus dem Gebiet der Elastizität. 


Beitrag zur näherungsweisen Ermittlung der Biegeeigenfrequenzen 
mehrfach abgesetzter Wellen. 


Von M. Sehilhansi in Heidenheim/Brenz. 

Die überwiegende Mehrzahl der Untersuchungen über Biegeschwingeungen von Wellen, 
welche das Problem analytisch behandeln, bezieht sich auf die „glatte“ Welle, d. h. solche 
Wellen, die durchlaufend einen und denselben Durchmesser besitzen. Mit den graphischen 
Verfahren dagegen ist es ohne weiteres möglich, auch die abgesetzte Welle zu berechnen: 
bei mehrfach gelagerten Wellen wird jedoch dieses Verfahren langwierig und bei der Über- 
lagerung der einzelnen Biegelinien, die zur Ermittlung der statisch unbestimmten Auflager- 
kräfte notwendig sind, auch ungenau. In der „Hütte* (3. Abschnitt, III. Kapitel, Unter- 
abschnitt IIea auf Seite 441 und 442 der 26. Auflage) ist ein analvtisches Verfahren mit- 
geteilt, das auf einer Reduktion auf die „glatte* Welle beruht, aber lediglich die niedrigste 
Eigenfrequenz einigermaßen zuverlässig liefern kann. 

Ich möchte nun ein neues Verfahren in Vorschlag bringen, das in diesem Auszug nur 
an dem Beispiel einer masselosen Welle mit einer Einzelmasse in der Mitte aufgezeigt werden 


1, Kine ausführliche Darstellung erscheint im „Ing.-Areh.“* 9, Heft 6. 
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soll. Die Welle seı zweifach also statisch bestimmt velagert und setze sıch aus mehreren 
Stiieken zusammen, deren Durchmesser voneinander verschieden sınd. Beı dieser Problem- 
stellung genügt für die Berechnung der kigenfrequenz die Kenntnis der Masse und der 
elastischen Konstanten. Letztere ist aber nichts anderes als diejenige Kraft, die an der 
Stelle der Masse angebracht werden muß, um die Masse um die Längeneinheit zu verschieben 
oc p anders ıuseedrückt . der reziproke Wert der elastischen Konstanten Ist die Verschiebung f 
des Massenpunktes, wenn an der Stelle der Masse die l,.asteinheit angreift Um dıe Ver- 
ehiehunz f zu berechnen, gehe ich von der Differentialgleichung der elastischen Linie aus: 
d’y 


Zr M (1). 


Mit Rücksicht darauf. daß für die abresetzte Welle J nicht konstant, sondern eine Funktion 


der Abszisse .r ıst, schreibe ıch foleendermaßen 


4 ‘ ‘) 
tdx ’ 1) 
llıer ıst ./. ireend eın Bezuesträrheitsmoment. Um auf der rechten Seıte eıne mit den 


üblichen und jedem Ingenieur zugänglichen Methoden ıintegrierbaren Ausdruck zu erhalten, 
entwiekle ich J,/J, in eine Reihe, von der ich verlangen muß, daß sıe den Verlauf der J,/J, 
nicht nur in der Umgebung eines beliebigen Punktes, sondern über die ganze Wellenläng« 
out annähert. Es kommen also Fourier-Reihen und Reihen mit Kuzelfunktionen ın Betracht. 


letztere sınd empfehlenswerter Man wählt als Länzeneinheit den halben Lagerabstand: dann 


hat das linke Auflaeer die Abszısse | und das rechte dıe Abszısse j )as Morent vom 
linken Auflaeer bis zur Last ist 1/2 (1 + ») und von der Last bıs zum rechten Auflager 1/2 (] +) 
Die Reihenentwieklung selbst liefert, wobeı P,„ die Kugelfunktion »n-ter Ordnung bedeutet: 
Ar > a„P., (.) 

() 


Bezüglich der numerischen Berechnung der Koeffizienten a,, die aus der allgemeinen Vorschrift 
(dl 3,3 a dr ; (4) 


(oleen. wird auf die Funktionentafeln von Jahnke und KEmde (Verlage Teubner S. 79 ff.) 


verwiesen. Setzen wir in Gl. (3) die expliziten Ausdrücke für dıe Kugelfunktionen 


|-3:9...-(2 mn 1) nn 1) nn lin - 3) B 
pP (| gr r" 4 = ge j l (.)\ 
n! Fl |) 3.4.2 n lıl?!n 3) 


lo 


ein, multiplizieren aus und ordnen nach Potenzen von x. so erhalten wir 
JlJ »’ u.” (6), 


wobei man allerdines in Kauf nehmen muß, daß die Koeffizienten b, selbst wieder unendliche 
Summen bestimmter Vielfacher von a, sind und ihr Genauigkeitsgrad davon abhängt, wie 


nt 


viele Glieder der Kurelfunktionenreihe man berücksichtigt hat. 


Die Gleichung der elastischen Linie lautet also für den Ast von | bıs O 
dd 7 | ‚ re 
F. — (] I r) \ I), a (71 
dr ) 
oder 
dl? 7 v 
2EJ - \ u As; | (S), 
od: 
wobeı 
Cm b, pe h,, Wr i . ; 5 : : : j : N 4 (9) 


ist. Bei der Integration treten für jeden Ast 2 Integrationsfreiwerte auf, mit deren Hilfe 
man den Verlauf der elastischen Linie den Randbedingungen ( Auflagerbedingungen) anpassen kann. 

Kontinuierlieh mit Masse belerte Wellen mit variablem Trägheitsmoment lassen sich 
mit diesem Verfahren bereehnen, indem man die kontinuierliche Massenbelegung in eine 
Reihe von gleichwertigen Einzelmassen zerlegt. Ebenso lassen sich mehrfach gelagerte 
Wellen damit untersuchen. Dabei zeiren sich erst die Vorteile des Verfahrens, weil die Ent- 
wieklune des Trärheitsmomentenverlaufs nur einmal gemacht zu werden braucht und bei 
allen Lastfällen. die man zur Ermittlung der statisch unbestimmten Lagerkräfte braucht, 
wieder verwendet werden kann. 3 
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Knickfestigkeit eines ebenen Systems beliebig vieler Stäbe, die in einem 
Knotenpunkte zusammenstoßen. 
Von Konrad Ludwig ın Hannover 


In einer Ebene liegen beliebig viele Stäbe, die in einem gemeinschaftlichen Knotenpunkte 
starr miteinander verbunden sind. Der andere Endpunkt jedes Stabes ıst ın einem Gelenk 
eelarert. Jeder Stab soll geren eın Ausknicken ın der Ebene fest sem Daher ıst nur eın 
Ausknieken senkreeht zur Ebene möglich, so daß jeder Stab zu einer ebenen Kurve verformt 
wird. Wählt man ın der Ebene durch den ausgeknieckten J-Stab das Gelenk als O-Punkt 
eines kartesischen Koordinatensystems, bezeichnet das 'Trägheitsmoment seines Querschnitte: 
in bezug auf die durch den (uersehnittsschwerpunkt gehende Achse parallel zur Kbene dei 
nicht ausgeknickten Stäbe mit 7, und seine Knickkraft mit 7’, so ıst die Differentialgleichung 
(des auseeknickten Stabes 


Mıt der Abkürzung 


wırd das alleemeine Integral der Differentialeleichung 


yj= A;sın A500; + DB; cos /;a 


Die Randbedineune im Gelenk ergibt 

yı= A;sınd;,&; 
Bezeichnet man die Länge des j-Stabes mit 1; und die Verschiebung des knotenpunktes ımıt f, 
so wird die Gleichung des ausgreknieckten Stabes 


/ 


Y; / Sin A; .# 


sın 7; 1, 


Hiernach ıst die Auskniekung proportional der Knotenpunktsverschiebung Daher ıst olım 
Knotenpunktsverschiebung kein Ausknieken möglich. Die Steigung im Knotenpunkte ıst 


d „H, 


’ / 
dar ;)® j 


/ 


| 











\bh. |]. Abb 


Wegen der starren Verbindung der Stäbe in ihrem gemeinsechaftlichen Knotenpunkte liegen 
die Tangenten aller ausgekniekten Stäbe im Knotenpunkte in einer Ebene, d. Ih. jeder aus 
geknickte Stab steht im Knotenpunkte auf der Normalen dieser Ebene senkrecht. Die Kben« 
der nicht ausgeknickten Stäbe wird als (z, y)- Ebene gewählt und die positive z. Achse so (Abh. 2) 
daß sich auf ihr der Knotenpunkt verschiebt. In diesem Koordinatensystem werden dis 
Riechtungscos jener Normalen mit «, # und w bezeichnet. Um die Riehtungseos des j- aus 
zekniekten Stabes im Knotenpunkte zu berechnen, wird seine Tangente im Knotenpunkte so 
verschoben, daß sie durch den O-Punkt geht, ın dem der Knotenpunkt vor der Auskniekung 


liegt. Der Winkel zwischen der Tangente und der z;-Achse wird mit a; bezeichnet 
ter a; fi;ctg 4;l;. 





LITE 
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Das 3-Kant aus der verschobenen Tangente, der »;-Achse und der x#-Achse ist rechtwinklig. 
Der Riehtungseos der Tangente gegen die x-Achse ist nach dem cos-Satz der Seiten 


CO)S Pj cos qa;j. 
Entsprechend ist der Riehtungscos der Tangente gegen die s- Achse 
sın Tj COS dj. 
Der Riehtuneseos der Tangente geren die z-Achse ıst 
sın d;. 


Wird in der (#;, 9;-Ebene das Quadrat der Steigung gegen -1 vernachlässigt, so vereinfachen 
sich die Rıehtungscos der Tangente zu 
COS qQ ; sın 1; to a;. 


} 


\us der Orthoeonalitätsbedinzung folet 


tv a; mw v). 


Wird die Anzahl der Stäbe mit 


COS 9 j HM 


sınyjr 


» bezeichnet. so ist der Range der Matrix 


COS sın pP to A, 
COS 5, sın (Vo tg As 
COS, sın Pr [g a, 
>» Demnach erzeben sieh die an — 2 Gleichungen 

COS q sin y, te a, 
COS G sing 0 Ig a, | == U) kB, H, us Wi 
COS gQ sın 77 to a, 

Ilierin werden für te a; die Werte eingesetzt 
COS P, sing ),etgi,l, | 
COS qQ, sıng R ), etg | | 0 
COS 77 sın war 7} ( ig Al; I} | 

Die Determinante wird ausgerechnet 


sin (4, p,)4,ctg 4,1, —sin(pE — 9.) 4, ctg /,l, + sin (9, — Y,)/r etg /, IK =V. 
Diese n > Stahilitätsbedineuneen ereeben zusammen mit den beiden Gleichzewichts- 


bedineuneen 


n Gleiehuneen für 


Der Fall u 


0. Band, Berlin 1 
Knickfestiekeit. Bi 


lie 


) 


» P COS 9 
— 


1 


n Knicekkräfte 


3, In dem nur eine Stabiılitätsbedineune auftritt. wurde von L. Vıanello. 
Kniekfestigkeit eines dreiarmigen ebenen Systems, Zeitschrift des Vereins Deutscher Ingenieure, 


YO, 
rlın 


falsch. da sie eeeen die Permutation 


f; 


| l: 


Seite 1753 behandelt. 
1921. Seite 109 versucht. 


Der 


Fall» 
Doch ıst dort auf Seite 110 die Gl. (9) 


3 wurde von R. Mayer, Die 


In 
f; | 


invarlant sein müßte’). Auch bestimmen die auseekniekten Stäbe mit ihrer Ebene vor dem 
\usknieken keine Pyramide. Weiter folgt dort statt n — 2 Stabilitätsbedingungen nur die 


eine Gl. (19). 


I Hinweis a 


o 


nke ich Herrn Prof. 


Ss 


Dr.-Ing. E. Kohl, Braunschweig. 
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Veranschaulichung und Anwendung der Minimalsätze 
der Elastizitätstheorie. 


Von Ü. Weber in Dresden. 


Der Satz vom Minimum der potentiellen Energie ist unabhängig vom Elastizitätsgesetze 
und auch für große Verschiebungen eültie'. Der Satz von Castieliano vom Minimunı 
der Spannungsenergie gilt nur, wenn die Gleichungen zwischen den Verschiebungen, Zerrungen 
(Dehnungen und Schiebungen) linearisiert sind, also 1. für kleine Verschiebungen (w, », ıw) und 
Zerrungen (e,y) und 2. nur für das Hookesche Gesetz. 

Gelten diese Einschränkungen, so lassen sich beide Sätze, da dann Überlagerungen 
möglıch sind, veranschaulichen. 

Abb. 1 zeigt als Beispiel einen unbelasteten Körper, der links an eine starre Wand 
befestigt ist. Die Querstücke a und 5b seien starr; Wand und a sind mit drei elastischen 
Bändern verbunden, die hier als Federn dargestellt sind, ebenso a und b. Die beiden oberen 
Bänder haben dieselben elastischen Eigenschaften wie die unteren. 


| 
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Abb. 3. Abh. 4 

















\bh. 5. 


Der Körper seı nun verformt und belastet. An den Punkten 1. 2 und 3 schreiben wir 
die Verschiebungen 


u 0 u, Ad, TA () 


1 , 2 e 3 
vor, am Punkte 4 bringen wir die Kraft X, = 9 an (Abb. 2). An allen übrigen Punkten sind 
ebenfalls die Kräfte (gleich Null) gegeben. Für alle Oberflächenpunkte sind hiermit entweder 
die Verschiebungen oder die Randkräfte bekannt. 
Zu den vorgeschriebenen Randbedingungen gehört ein bestimmter Verschiebungszustand 
aller Punkte; aus ihm lassen sich die Zerrungen und auf Grund des Hookeschen Gesetzes 


1) Vorausgesetzt wird, daß der belastete elastische Körper im stabilen Gleichgewicht ist. 
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die Spannungen und hieraus die Formänderungsarbeit A,,, ausgedrückt durch die Ver- 
sehiebuneen. aneereben. Die Formänderungsarbeit ist gleich der Arbeit der äußeren Kräfte. hier 


er l 2 | 


X, at 0-u 


(der Index »# bezieht sieh auf die „wirkliche Lösung“ ) 
7/u den vorzeschriebenen Randverschiebungen 


HM, 0, H. A. u, = 0 


lassen sieh noch beliebie andere Verschiebungszustände angeben. Berechnet man aus ihnen 
die Zerrungen und dann die Spannungen, so sind letztere weder untereinander noch mit den 
Randkräften im Gleichgewicht (Abb. 3). 

Um diesen Zustand aus dem wirklichen zu erhalten, sind zusätzliche Kräfte anzubringen 
und zwar innere, P’, um die inneren und äußeren, P”, um die äußeren Gleichgewichts- 
bedinzungen zu befriedigen. Den zusätzlichen Verschiebungen d u = u — u, allein entsprechen 
Zerruneen und eine Formänderungsarbeit Ö A, die gleich der Arbeit dieser zusätzlichen Kräfte 
und damit größer als Null ist; also 


| ’ [3 | ’ Zi 
VA Yu P PS a 


-) 


Nıin stellen wir uns vor. daß der neue Zustand aus dem wirklichen entstanden ist. be- 
rüeksichtiren die Arbeit der äußeren Kräfte und erhalten die Formänderunesarbeit des neuen 
Zustandes 


| 


A=Aut5>9wWP +z0du”P"”+0Q(u, —u,w) 


hieraus 

A Yu, Au -- Yu w 
und zwar, weil stets die positive Arbeit der Zusatzkräfte hinzukommt, wenn man aus der 
reehten Seite die linke bildet. 

Bezeichnen wir A - Q@-u, als potentielle Energie des Systems, so sehen wir, daß für 
den wirklichen Verschiebungszustand die potentielle Energie zum Minimum wird. 

Zur Darstellung des Satzes von Castigliano gehen wir von den Spannungszuständen 
aus. Bei Verformung und Belastung nach Abb. 2 entstehen die Spannungen 9,,, T,,, aus denen 
die Zerruneen mit Hilfe des Hookeschen Gesetzes berechnet werden können: aus diesen 
müssen sich Verschiebungen ergeben, die mit den vorgeschriebenen Randbedingungen in Ein- 
klang stehen. Aus den Spannungen bestimmen wir die aufgespeicherte Spannungsenergie 

Br tn, Ka) ar . X,.a- € Osin. 

Wir betrachten nun wieder den unbelasteten und unverformten Körper; das heißt, die 
voreeschriebenen äußeren Kräfte sollen gleich null und die vorgeschriebenen Randver- 
schiebungen gleich null sein. Bei diesem Zustande nimmt man an, daß alle Spannungen 
null sind. Es könnten aber auch Anfangsspannungen do und ör vorhanden sein. Diese 
müssen die Gleiehgewiehtsbedingungen erfüllen sowohl ım Innern, als auch an den Stellen, 
wo die Randkräfte null sind. An den Stellen, wo die Randverschiebungen null sind, könnten 
jedoch Spannungen auftreten. Zerlegt man den Körper dann in kleine spannungslose Teile, 
so passen sie nicht mehr zusammen. Aus den Anfangsspannungen ö60 und Ör läßt sich die 
Spannungsenergie Öd B berechnen. 

Uberlagern wir nun den Spannungs- und Formänderungszustand, der den Spannungen 
On, 7,, entspricht, so erhält man eine neue Lösung von Spannungen 6=0,+do und r=r,-+ Ör, 
die im Inneren im Gleiehgewiehte sind und am Rande den äußeren Kräften das Gleichgewicht 
halten. Umgekehrt kann man auch ein solehes neues Spannungssystem o,r annehmen und 
dann 96 = 0 —0,. dr=r --r,. berechnen. 

Abb. 5 veranschaulicht den Körper mit dem neuen Spannungszustande, Abb. 4 den 
/ustand, der den Anfangsspannungen entspricht. Eigentlich kann man den Verschiebungs- 
zustand gar nieht darstellen, da beim Zerlegen in kleine spannungslose Teile diese nicht mehr 
zusammenpassen. Hier ist das Vorhandensein von Anfangsspannungen dadurch angegeben, 
daß dureh Ansätze am Querstück a die Bänder vorgespannt werden. 

Der Spannungszustand nach Abb. 5 entsteht dureh Überlagerung der Spannungszustände 
nach Abb. 2 und Abb. 4. Beim Zustande nach Abb. 4 ist die Spannungsenergie dB >O 


vorhanden. 
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Nehmen wir allmählieh die Verformune und Belastune nach Abb. 2 vor, so wird ım 


Punkte 4 die äußere Arbeit 5 0 u, „ und im Punkte 2 die äußere Arbeit ZA,» a + (X,  A,n) a 


dem Körper zugeführt, so daß zum Schlusse die Spannungsenergie B(o,r) im Körper auf 
gespeichert ist: 


| Bi 
RBR=) R 5 1, () U r 5 Au a—+(X, RT 7 
Da 
| | 


Ju 


ist. so folet 


B—-PF,& Bas —- P,»a. 


Ist insbesondere a =0, so folgt, daß für den wirklichen Spannungszustand die Spannungs- 
energie B(o,r) zum Minimum wird. 


Anwendungen. 1. Es lassen sich aus den Minimalsätzen die jeweils bei der Herleitung 
nicht benutzten Grundgleichungen der Elastizitätslehre zurückzewinnen, also aus Ad(u, v, m) 
die Gleichgewichtsbedingungen und aus B(o,r) die Beziehungen zwischen &,y und u,e, m. 

2. Es lassen sich zweekmäßige Näherungslösungen aufstellen. 

3. Für Stabilitätsprobleme kann nur der Satz vom Minimum der potentiellen Energie 
verwendet werden, wobei die Gleichungen nicht linearisiert werden dürfen. 

t. Mit Hilfe beider Sätze läßt sich die Verschiebungsgröße des Angriffspunktes einer 
KEinzelkraft, bzw. der Angriffsebene eines Kräftepaares eingrenzen. Dieses Verfahren läßt 
sich erweitern auf die Bestimmung der Verschiebung eines beliebigen Punktes. 


Beides werde an Beispielen gezeigt. 


Abh. 6 zeigt einen zweimal gestützten Balken, der in der Mitte durch die Kraft P 
belastet ist, von der Höhe h und der Breite b<h. In den Stützquerschnitten und im Be- 
lastungsquerschnitt seien versteifende Wände vorhanden, die in Richtung der äußeren Kräfte 
starr sind, so daß alle Punkte dieselben Verschiebungen ® haben. In Riehtung der Balken- 
achse kann die Verwölbung w beliebig sein; die Verschiebungen « senkrecht zur Bild- 
ebene sind nur gering und werden so bemessen, daß o,—=0 ist. Durch Einführung der ver- 
steifenden Wände wird die Frage nach den Verschiebungen sinnvoll. 








Abh. 6%. Abb. 7. 


iz; zZ 


Der Verschiebungsansatz besteht aus folgenden zwei Teilen w’, »’, w’ und u”, »"’. w 


l. Die Mittellinie der linken Balkenhälfte hat die Verschiebungen nach unten: 


3 


> r m 

a e : | +.) w | | 4 | 

m .) ») F 

5 2\L/ I 

In gleicher Weise mögen sich alle anderen Punkte eines jeden @Querschnittes verschieben. 
Die Querschnitte stellen sich senkrecht zur Mittellinie ein, hieraus 


0 mm’ Q Pi A) ’ 


’ 
n u —. .q 
7] > 27 
(y von der Mittellinie aus gerechnet). 0, =0, hieraus die geringfügigen «'. 
2, Die Querschnitte stellen sich schräg zu der Mittellinie nach Abb. 7. 


m 


’ 2 - 
" " u —=Q, w =Ü,. 


In / 


Im Gesamtansatze für die Formänderune befinden sieh die beiden noch wählbaren Bei- 
werte %,, und vun. 
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Man erhält die Formänderunesarbeit 


ni 3JE m ur ) 
| ’ -F@GlIl Int F h.h: J = h* 


r 3.) F. RR FF fı ] : er ] | ] 
Di . Pa. -# (tm + 9%. A, Fa On u m 
(] It" / Ä ) 7 
die wirkliche Durchbiegung des mittleren (@uersehnittes). Die günstigste Lösung 


in Hu 


erhalten wir, wenn wir die linke Seite möglichst klein machen. 
können Be und 4. bestimmt und eingesetzt werden. Man erhält 


Aus dieser Minimalbedineunge 


| | | Bi 11” , | 
DR? JE > FG) 3 rm 
Hieraus 
| Pa | 1”) 4 
LiZur JE 2FG 


Spannungen. Die Querkraft werde durch 


Nun machen wir einen Ansatz für die 
Die Normalspannungen o, werden nach 


parabolisch verteilte Schubspannungen übertragen. 
der gewöhnlichen Elastizitätslehre bestimmt. 

Für die einzelnen Balkenteile läßt sich hieraus die Formänderung bestimmen. In der 
Mitte passen Jedoch wegen der Verwölbung die Querschnittsflächen nicht zusammen und wegen 
der Querzusammenziehung nieht mit der eingefügten Querwand. 


Der Ansatz eıbt 
ir 1.35 
} P 


mi 6 JE >2FG 
Dıe Ansätze lassen sıch selbstverständlieh beliebige verfeinern. 
", ıst hiermit zwisehen zweı Werten eineeerenzt. 




















Nun soll ein Balken nach Abb. S durch zwei Lasten P, und P, belastet sein. 
Durch die Kraft 7’, verschiebt sich der Angriffspunkt von P’, um e,,P, 


von 7 mr... 


Kraft 7, verschiebt sich der Angriffspunkt von 7’, um e,,F, 


Durch die 
von P, um e,,JP,. 


Die gesamte Verschiebung ®, ,. Ist also ce, P,+ec,FP;., 
ebenso ou - ,ı. 
(sresucht wird der Wert e 


Der Satz vom Mınımum der potentiellen Energie gibt mit einem Ansatze ähnlich dem 


vorherieen: 


(l e)l’(l, +1)”  1(,—+1,) 


ai? 5 GEJI >GFl 
I e)l,i, ?—I1?’—1?) 4 
Bi IP EJI DGr! 


1—«)(l, +11? (d,-+1,)-7 
2 6GEJI  2GFl 
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Der Satz vom Minimum der Spannungsenergie gibt 
| >) >) | » pP? | \ PP PP’ | » P? \ pP: ı»Dp PP I’ ) PP .) 
) P, vw ro9 P,v,w > %ıı ı Tea tıde ra an ko Pııı T=Pı ) ud}; (-) 


wobei sich die 5 von den a dadureh unterscheiden, daß bei den ersten Gliedern (1 €) fort- 
fällt. bei den zweiten der Faktor 1, 2 hinzutritt, 


Aus (1) und (2) ist e,, so zu finden, daß es zwischen zwei Werten eingegrenzt ist. 


In Ungleichung (2) setzen wir einmal P, gleich Null, dann 7’, gleich Null und ziehen 
die so erhaltenen Ungleichungen von Ungleichung (1) ab, so daß nur das Glied e,, P, P, links 
nachbleibt und das Ungleichungszeichen bestehen bleibt. 


ee 
Auf dieselbe Weise erhält man dureh Vertauschung der Ungleichungen (1) und (2) 


2 (P IP FE ai. Hr 


11 11 1 ı nu 
Wir teilen die Ungleiechungen (3) und (4) dureh P, P,, wobei die Ungleiehungszeichen bestehen 
bleiben, wenn wir für P, und P, dasselbe Vorzeichen annehmen, andernfalls sind sie zu 
ändern. Beide Fälle sind getrennt zu untersuchen. Der erste Fall zıbt 


) > 1 
’ ) { £ m ( ) a TE '\ 
Ca (Pr kıı IP l’22 22 P ’ 
> > l 
> ) Pr .) } Al ( N (; 
C 19 (Pi dp = Dis T (Pas Ass ıp ä 


Das Verhältnis P/P, ist so zu wählen, daß die rechten Seiten jeweils zum Extremum 
werden. Das gibt für beide Ungleichungen 


P, Pa "2 


n Pıi A 
wobei die Wurzel für Fall 1 mit dem positiven Vorzeichen zu nehmen ist. Man erhält 
20 Ey Pr - u) Pa a) El Ef r2V ld, — ar) (Das Mr)» - » ID). 
Wählt man für P, und 7’, verschiedene Vorzeichen, so erhält man 

2 PP 2) (Pi: A) (Pa: (yo) Cs 2 Ga j 27 (Pi A) (Pa: Ayo) - . (6). 


Von beiden Ungleichungssystemen ist dasjenige zu wählen, das den Wert e,, stärker 


einerenzt. 


r 
Für unser Beispiel ist 5, größer als a,,, so daß Ungleichungssvstem (6) in Betracht 
kommt. 


Für den symmetrischen Belastungsfall mit !,=1!, erhält man 
1? —21°)—-2uW1?’(l—1)” 1?—-021 (1-21) 
GEJI GHFI 


,’(?—21’) + ul? —21)% 1?+021,(1-—1) 
GEJI G Fl 


Auf diese Weise kann man sowohl die Verschiebung des Angriffspunktes einer Kraft 
als auch eines anderen Punktes eingrenzen. Die Eingrenzungsspanne läßt sich durch Änderung 
der Ansätze beliebig verfeinern, wobei man die Verbesserung der Ansätze mit Hilfe der An- 
schauung vornehmen kann. 393 
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BUCHBESPRECHUNGEN 


ERNST BITTNER, Momententafeln 


, 


ınd Kinflußflächen für kreuzweise be 
wel 4 ‚Ist nbetonplatten. Mit 16 Textabb. 
ı, 81 Zahlentaf. Wien 1938. Verlae von Julius 
Springer. Preis geh. 
Während die Bereehnunz von Platten für Vollast 
mit Hilfe der Mareusschen DBeiwerte » für prak 
tische Bedürfnisse genügend genau und zenügzenid 
einfach möglich ist, ist ein ähnlich einfaches Ver 
fahren für Einzellasten bisher nieht zereben worden. 
Hier wird dem Ingenieur ein Tafelwerk zur Ver 
fürung gestellt. das bei hinreichend engen Tafel 
hri Platte unter wandernden Einzellasten 
ı berechnen zestattet: für «die Herstellung der 
lafeln sind die zenauen Erzebnisse («der Platten- 
theorie verwertet, für den Gebrauch der Tafeln ist 
es aber nicht erforderlich. die Theorie selbst zu 
kennen. Die Kinzellasten werden als Vollasten über 
kleine Reehtecke von bestimmten Abmessungen 
eineeführt. Am Schluß des Buches wird die Platten 
rie im Anschluß an Nädai) so weit entwickelt. 
laß eine Nachprüfung und allfällige Ergänzung der 
lafelwerte jederzeit möglich ist. Die Lösung für 
len eingespannten unendlichen Plattenstreifen und 
für die unendliche Kragplatte ist vom Verfasser 
selbst entwickelt worden. Das Werk wird für viele 
Zwecke der Praxis von Nutzen sein, die eine über 
siehtliche Verwertunze der exakten Theorie ver 
angt, ohne in jedem Einzelfalle auf deren Entwick 
luneen Bezur zu nehme: 


scehritten ine 


tr’) 


Karlsruhe. T 

H. ERTEL, Methoden und Probleme 
ler dvnamischen Meteorolorie (.Er- 
rebnisse der Mathematik und ihrer Grenzzebiete“. 
herauszereben v, d. Schriftleitune des .Zentralblatt 
für Mathematik“, Bd.5. Heft 3). 122 S. m. 14 Fie. 
Verlag Julius Springer. Preis brosch. 


Berlin 1938, 
I4 M. 

Die dynamische Meteorologie will an Hand von 
(sedlankenmodellen das physikalisch Wesentliche des 


(seschehens in der Erdatmosphäre herausarbeiten. 
dab (diese Aufeabe im alleremeinen schwerer ist als 


lie entsprechende in manchen 
chen Anwenduneseebieten, ist 
man bedenkt. daß die Versuchsbedineuneen in der 
\tmosphäre der Beeinflussung ent 
ven sind. Es ist sehr dankenswert. daß H. Ertel 
jetzt einen Bericht vorlegt, der in gedrängter Kürze 
nen Überblick über dies entwieklunesfähire Ge- 


biet zibt. 


andern physikali- 
sofort klar, wenn 


menschlichen 


Im ersten Kapitel wird zunächst die 
ivnamik trockener und feuchter Luft. 
| dem benötirten Umfane 
die hydrodynamischen Gleichungen 
in Eulerscher und Lazraneescher Form aufzestellt. 
Im zweiten Kapitel über die allgemeine Dynamik 
ler Atmosphäre werden insonderheit die beiden 
Tatsachen in die Betrachtunr einzeführt. daß die 
heinungeen in einem rotieren 
en Koordinatensvstem beobachtet werden und daß 
ırch die vertikale Koordinate alus- 
rezeichnet ist. Dieses Kapitel gzipfelt in der Be- 
les Zirkulationstheorems von Bjerknes 
Enereiebilanz. Bei der 


Probleme im dritten Kapitel 


Thermo- 
sowie die 
Q dareelert. 


| 
Danach werden 


eltOrTOlOPIS«ı en Frse 


Schwere die 


atmosphärischen 


« 
1) I, . y Naar 
behandlung spezielle 


ist in der atmosphärischen Statik auf die thermo- 
Ivnamische Ableitune- der Bedinrune für stabile 
Schiehtune hinzuweisen. Natürlich finden in der 


atmosphärischen Dynamik die Windverteilune über 


lem Erdboden, die eng mit dem Turbulenzproblem 
‚usammenhängt, ferner die Theorie der Diskontinui 
täten sowie die atmosphärischen Störungsgleichun 
een und das Zyklonenproblem eine ausführliche 


Besprechung. 


Der Verfasser hat es ausgezeichnet verstanden, 
einen Überblick über die Erzrebnisse und Probleme 
eines weitverzweigten Forschungsgebietes zu geben. 
Besonders wertvoll sind die ausgedehnten Literatur 
aneaben. 

Dresden. W. Tollmien. 81 

Grimsehls Lehrbuch der Physik, neubearb. von 
Professor Dr. R. TOMASCHEK, Direktor d. Phys. 
Inst. d. Techn. Hochsch. Dresden, II. Bd., 2. Teil, 
Materie und Äther, 8 Aufl.. VIII + 456 >. 
m. 339 Abb. Leipzig und Berlin 1938, Verlag 
PB. G. Teubner. Preis geb. 14 M. 

Der vorlierende Schlußband des Grimsehlschen 
Lehrbuches ist eine völlige Neuschöpfung des 
Herauszebers: die darin behandelten Probleme 
standen zu Lebzeiten Grimsehls, der ja im Ok- 
tober 1914 als Vierundfünfzigjähriger in Flandern 
vefallen ist, bestenfalls in den ersten Anfängen 
ihrer Entwieklung. Die Eigenart des alten „Grim- 
sehl‘“ ist aber auch in diesem Bande zewahrt ze- 
Hauptgewicht wird auf die Er- 
fahrunzstatsachen zelert und erst im Anschluß 
daran Theorie zebracht. Wenn freilich der Ver- 
fasser im Vorwort betont. er habe sich bemüht. 
„niemals Tatsachen als Bestätizunzen der Theorie 
auszureben“. so ist das wohl eine etwas mißver- 
ständliche Redeweise: Richtiekeit und Geltungs- 
bereich einer Theorie können schließlich doch nur 
an Tatsachen geprüft und durch sie erwiesen 
werden. Auch aus der Darstellunz des Verf. geht 
übrieens hervor, daßer dieses Verhältnis zwischen 
Theorie und Erfahrung durchaus anerkennt. In 
liesem Zusammenhanze sei auf zwei kleine In 
korrektheiten hingewiesen. Die S. 200 zerebene 
I>arstellunz der Abhängigkeit der Masse von der 
(eschwindiekeit kann dahin mißverstanden wer- 
den. als ob jene Formel induktiv aus Versuchen 
oefoleert worden sei, ] 


hlieben: das 


während ja in Wirklichkeit 
die von der Relativitätstheorie zelieferte Formel 
an den Experimenten zeprüft und durch sie be- 
stätiet wurde. LUmeekehrt wird auf S. 281 der 
lLandesche Aufspaltungsfaktor als Ergebnis von 
Reehnunsen zewonnen: hier wäre der Hinweis 
anzebracht wzewesen. daß Lande, dessen Name 
ibrirens in diesem Zusammenhang gar nicht er- 
wähnt wird. den Faktor durch eine zlückliche In- 
duktion aus Experimenten erschlossen hat. 

Nun zum Inhalt des Buches. Man muß ihm vor 
allem nachriühmen, daß ein ungemein reichhaltiges 
Tatsachenmaterial in die Darstellung hineinge- 
arbeitet worden ist. Auf dem verhältnismäßig 
beschränkten Raume von 450 Seiten wird über 
eine solehe Fülle von Einzeluntersuchungen be- 
richtet, daß. wenn das Buch auch kein Lernbuch 
im eirentlichen Sinne des Wortes ist, es doch. 
namentlich auch dem älteren Studierenden. eine 
auseezeichnete ÖOrientierunz und ein tieferes Ein- 
drineen in die darin behandelten Gebiete er- 
mösrticht. Vortreffliche Abbildungen und an- 
schauliche Skizzen. zahlreiche Diagramme und 
reichhaltires Tabellenmaterial — das alles be- 
sonders auch deshalb eindrucksvoll, weil vielfach 
den ÖOririnalabhandlunzen entnommen —- unter- 
stützen den Leser beim Studium aufs wirksamste. 
Iazu kommt, daß neben den experimentellen 
Forsehunzen auch die neueren Theorien durchaus 
zur Geltunz eekommen sind. Es wird nicht nur 
das klassische Bohrsche Atommodell — in. wie 
mir scheint. etwas übertriebener Pietät als Lenard- 
Rutherford-Bohrsches Modell bezeichnet — mit 
aller Gründlichkeit und Klarheit behandelt. son- 
dern die schöne. in sich zeschlossene Darstellung 
des Abschnitts „Welle und Korpuskel“ dringt bis 
zur Schrödingergleichung und ihren wichtigsten 
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Anwendunren vor und vermittelt auch dem Leser, 
der nicht über tiefere mathematische Vorkennt- 
nisse verfürt, eine auszezeichnete Vorstellung von 
der wellenmechanischen Deutung der Vorgänge in 
Atom und Molekül sowie der radioaktiven Vor 
eänee, wobei auch die Grundzüre der Gamow 
schen Theorie des ß-Zerfalls gezeben werden. Ge 
lewentlich wirkt sich m. E. die sorgfältige Ver- 
meidune der Relativitätstheorie zum Nachteil einer 
sinnvollen Deutung der Erscheinungen aus. Dab 
Verfasser auf dem Boden der Äthervorstellung 
steht, was ja auch der Untertitel andeutet, ist für 
die meisten Kapitel ohne Belang, da dieser Äther, 
der zwar Substanz (d. h. Träger von Eigenschaf- 
ten). aber nicht materielle Substanz sein soll, der 
aber andrerseits „körnige Struktur“ besitzen soll, 
die sieh im Auftreten des Planckschen Wirkungs- 
quantums äußere, zur Deutung der krscheinungen 
ebensowenig herangezogen wird wie der an zwei 
Stellen erwähnte Lenardsche „Uräther“. Wie in 
den vorhergehenden Bänden finden sich auch in 
diesem Band biographische Angaben über eine 
Reihe von Forschern: die Auswahl scheint etwas 
zufällirer Natur gewesen zu sein. Leider sind die 
Angaben wesentlich kürzer und farbloser gehalten 
als namentlich im ersten Bande, wo sie meist 
recht gute Vorstellungen vom Wesen und Werk 
der einzelnen Forscher vermitteln. Angaben 
über Ernst Grimsehl, den Schöpfer des nunmehr 
abzeschlossenen Werkes, fand ich in keinem der 
drei Bände. 

;erlin. E. Mosch. 857 


s 


Dr. phil. HARALD GEPPERT, o. Prof. d. Math. 
a.d. Univ. Gießen, und Dr. phil. habil. SIEGFRIED 
KOLLER, Vorstand d. Statistischen Abt.d. Kerck- 
hoff-Instituts Bad Nauheim, Erbmathematik. 
Theorie der Vererbune in Bbevölke 
rune und Sippe. VII + 236 S. m. Abb. Leip 
zie 1938, Verlag Quelle & Meyer. Preis geb. 18 M. 


In diesem Buch. das aus der gemeinsamen 
Arbeit eines Mathematikers und eines Biologen 
hervorgegangen ist, werden die wahrscheinlich 
keitstheoretischen Zusammenhänge im biologi- 
schen Geschehen von Familie, Sippe, Volk und 
Rasse dargestellt; die statistischen Methoden der 
Biologie bleiben also außer Betracht. Die „Erb- 
mathematik“ kann den Fachleuten beider Diszi 
plinen bestens empfohlen werden: Der Biologe 
tindet darin eine ausführliche Entwicklung der für 
seine Arbeit wiehtigen Verfahren der Wahrschein 
lichkeitsrechnung, und dem Mathematiker bietet 
das Buch von Geppert und Koller vielerlei An 
rezung, an der Lösung dringlicher und zugleich 
mathematisch reizvoller Aufgaben der Biologie 
mitzuarbeiten. Allerdings bin ich nicht davon 
überzeugt, daß ein Leser, der nur die üblichen 
mathematischen Schulkenntnisse mitbringt, der 
Darstellung überall zu folgen vermag (wie die 
Verfasser meinen); er muß zum mindesten die 
Grundlagen der Wahrscheinlichkeitsrechnung be 
herrschen und darüber hinaus mit dem Rechnen 
mit Zahlenfoleen etwas vertraut sein. Für den. 
der den Vormarsch des Mathematikers in die Be 
reiche der Biologie mit Mißtrauen beobachtet, sei 
ausdrücklich bemerkt, daß die Verfasser der 
„Erbmathematik“ keineswegs eine abstrakte 
mathematische Theorie gezüchtet haben, sondern 
iiberall in engstem Kontakt mit der biologischen 
Wirklichkeit bleiben; das zeigen die gewissenhafte 
Untersuchung der Grenzen, innerhalb deren die 
Theorie den tatsächlichen Verhältnissen gerecht 
zu werden vermag, die vielfache Prüfung .der 
Theorie durch Gegenüberstellung mit Beobach- 
tungsergebnissen und die Anwendung auf viele 
Fragen von unmittelbar praktischem Interesse 
(Verhütung erbkranken Nachwuchses, Inzucht, 


Rassenmischung und -entmischung, Abstammungs 
prüfung). Die „Erbmathematik“ ist folgender- 
maßen zegliedert: I. Die Mendelschen Regeln. 
II. Das Erbeefüre einer Bevölkerung bei völliger 
Durchmischune. III. Natürliche und künstliche 
Auslese in einer Bevölkerung. IV. Das Erbgeefüre 
einer Sippe in einer beständireen Bevölkerung. 
V, Erbbegutachtunge. Die Beiträge des Mathe 
matikers und die des Bioloren sind so eut auf 
einander abgestimmt, daß eine „Abstammungsprü 
fung“ hier gewiß nicht zum Ziel führen würde! 
Freiberg (Sa.). Grüß. 878 


Dr. SIEGFRIED VALENTINER, Prof. f, Physik 
a. d. Bergakademie Clausthal, Vektorana 
Iysis (Sammlung Göschen Bd. 354). 5. erneut 
durehgzesehene Aufl., 136 S. m. 13 Fir. Berlin 1938. 
Verlag Walter de Gruyter & Co. Preis ceb, 1,62 M. 

Die vorliegende Auflage des als erste Einführun» 
in das Gebiet der Vektoranalysis durchaus zu emp 
fehlenden Bändehens ist im wesentlichen ein un 
veränderter Abdruck der in Band 10 (1930) S. 309 
dieser Zeitschrift besprochenen vierten Auflage. 

Willers, 89 


CZUBER-BURKHARDT, Die statistischen 
Forscehungsmethoden. 3. erweit. Aufl. 
Herausgeg. von F. Burkhardt, Prof. a. d. Univ. 
leipzig. XVI-+330 S. m. 38 Fig. i. Text. Wien 
1958, Verlage IL. W. Seidel u. Sohn. Preis reb. 
13,50 M. | 

Die deutsche Literatur ist nicht reich an guten 
lwehrbüchern der mathematischen Statistik. Um so 
mehr ist es zu berrüßen, daß das vielbenutzte 
Werk von E. CUzuber nunmehr in dritter erweiter 
ter Auflage, bearbeitet von F. Burkhardt-Leipzig, 
vorliegt. Die Anlage des Werkes sowie die mathe 
matische Darstellune ist im wesentlichen dieselbe 
seblieben (vgl. die Besprechune Zamm 3, 1923, 
S. 827). Die Erweiterungen betreffen u. a. die Be 
völkerungsstatistik, die Zerlegung von Zeitreihen 
und geldliche Ausgleichsprobleme. Wünschenswert 
wäre eine Berücksichtigung der Anwendungzen der 
Statistik in der Technik und Industrie (Probleme 
des Fernsprechverkehrs, Verfahren der Fabrikati 
onsüberwachung). 

Berlin. Günther Schulz. 844 


Ferner sind bei der Schriftleitune folgende 
Bücher einzeraneen ausführliche Besprechung 
bleibt vorbehalten): - 

Professor PAUL B. FISCHER, Studienrat am 
Gymnasium Berlin-Steglitz, Arithmetik (Samm 
lunze Göschen, Nr. 47). 152 8. m. 19 Abb. Berlin 
1938, Verlag. Walter de Gruyter & Co, Preis 1,62M. 


_— 


Der feste Körper, Vorträge a. d., Tagung 
der Physik. Gesellsch. Zürich anläßlich der Feier 
ihres SVjährigen Bestehens, m. einem Vorwort von 
RR. SANGER, Präsident der Physikal, Gesellsch. Zü 
rich. V + 154 S. m. 30 Fig. Leipzig 1938, Verlag 
S. Hirzel. Preis kart. 7 M. | 


Dr. GERHARD KOWALEWSKI, 0. Prof. a. d. 
Techn. Hochschule Dresden. Grundbeeriffe 
und Hauptsätze der höheren Mathe 
matik insbesondere für Inzenieure 
und Naturforscher. 156 S. m. 40 Fig. Ber 
lin 1938, Verlag Walter de Gruvter & Co. Preis 
seh. 5M. 


Dr. Inze. FRIEDRICH TÖLKE VDlI, o. Prof, a. d. 
Techn. Hochschule Berlin, Talsperren, Stau- 
dämme und Staumauern (Wasserkraftan 
laren 2. Hälfte, I. Teil der Handbibliothek für Bau 
inzenieure, Ill. Teil. 9 Bd.) XI+ 734 S. m. 1189 
Abb. im Text. Berlin 1938, Verlag Julius Springer. 
Preis geb. 78 M. 
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NACHRICHTEN 


Internationaler Kongreß für angewandte Mechanik Hauptversammlung der Lilienthal-Gesellschaft für 


ER RÄT 


in Cambridge (Mass.). Luitfahrtiorschung. 
PR ie internationalen Kongresse für Angewandte Die diesjährige Hauptversammlung der Lilienthal- 
Ri Mechanik. «ie jetzt in vierjährigem Abstand statt Gesellschaft für Luftfahrtforschung fand vom 12. 
" finden, sind seit ihrer Beeründunge im Jahre 1924 bis 15. Oktober in Berlin statt. Die Versammlung 
3 nicht mehr aus der Entwicklung dieser Wissen- war wesentlich stärker besucht als die vorangehen- 
4 schaft fortzudenken. Der fünfte dieser Kongresse den. Von den allgemeinen Vorträgen werden den 
& tagte in diesem Jahre vom 12. bis 16. September in l,eserkreis dieser Zeitschrift interessieren: 1. I. Si- 
3 (ambridee bei Boston im S>Staate Massachusetts korskv: Das Großflueboot. E. Heinkel: Erhö 
$ ler Vereinieten Staaten von Amerika. (Gasteeber hunze der Geschwindigkeit (der Flugzeuge in den 
£ waren «das Massachusetts-Institut of Technology letzten Jahren, W. Findeisen und L. Ritz: Verei- 
€ und die Harvard Universitv, In der sehon traditio- sunz, X. Ferri: Untersuchuneen und Versuche im 
4 nellen Weise fanden weeen der Überfülle der Vor- Ubersehallwindkanal zu Guidonia. Von den Sek- 
träre eleichzeitie drei Sektionssitzungeen statt. Die tionsvorträgen seien die foleenden besonders hervor 
erste Sektion beschäftiete sich mit Hlastizität und eehoben: ©. Minelli: Das energetische Verfahren in 
Festiekeitslehre. die zweite mit Hydrodynamik, der Flugzeugstatik, F. Bollenrath: Zeit- und Dauer- 
während in der dritten Sektion sonstige Fragen der festirkeit der Werkstoffe, ©. Nissen: Festiekeits 
Mechanik. wie Schwingungen, Schmierung, Flug- fragen bei der Gestaltung neuzeitlicher Flugzeuge, 
‚eugedvnamik, Wärmeübertragung, Flußbaufragen W. G. A. Perring: Neuere Luftschraubenuntersu- 


und manches andere behandelt wurden. Außerdem chungen in England mit besonderer Berücksichtigung 
fanden gemeinsam für alle Sektionen sechs allge des Startproblems, F, Haus: Aerodynamische Grunid- 


meine Vorträge statt, deren Themen angedeutet lagen der selbsttätigen Stabilisierung, H. W. Kanl: 


werden mögen. Sverdrup sprach über Ozeanzirku- Statistische Erhebungen über Betriebsbeanspruchun 
lation. Chatfiell über Aneewandte Mechanik und ven von Flugzeugflügeln, E. Schmidt: Die gra 
Juezeuebau. Tavlor über neuere Entwieklungen im phische Bereehnung der Vergleichsprozesse von Ver- 
"urbulenzproblem, Körber über das Verhalten me brennungsmotoren unter Berücksichtigung der Tem- 
tallischer Werkstoffe bei kleinen Verformungen, peraturabhängirkeit der spezifischen Wärmen. 
erös über die Verwendune von Gleichnissen (z.B. C, Cremona: Bericht über Untersuchungen und Ver- 
elektrischer Art) in der Angzewandten Mechanik und suche im Schleppkanal zu Guidonia, H. N. Garner: 
schließlieh Hardv Uross über Probleme der Bau- Neuere Versuche an Flugbooten und Flugbootmo- 
konstruktionen. "In den Sitzuneen ler ersten Sek dellen, >. Truseott: Mitteilungen über den ver 


erößerten NAUA-Schleppkanal und hier ausgeführte 


tion standen die elastischen Stabilitätsprobleme und 
Versuche, W, Sottorf: Start und Landung im Mo 


lamit «die Ansätze für endliche Verformungen sehr 


im Vordererund. In der zweiten Sektion wurde ddellversuch. Die Vorträge werden z„esammelt 
a las Turbulenzproblem besonders ausgiebig erörtert, herauszerreben. 
ER wofür eigens eine mehrtägire Aussprache in dem Berlin. l. Flügge-Lotz. 912 
Di „Turbulenee Symposium“ vorgesehen war. Erfreu- 
| licherweise ist es diesmal wieder möglich, längere Der a. o. Prof. a. d. Technischen Hochschule 
= \uszüge der Vorträge und sogar einen Teil der Darmstadt. Dr. Kurt Schmieden, wurde zum 
> Diskussionsbemerkungen durch den Druck der wei o. Prof. ler Mathematik an dieser Hochschule er 
FE teren Fachwelt zugänglich zu machen. nannt. 
r ? \n (len Kongreß schlossen sich ausgedehnte Be- Der Dozent für Mathematik und Grundlagen der 
th sichtigungen an. In New York wurden bemerkens "heoretischen Physik an der Universität Köln, Dr. 
2. werte Inzenieurbauten besiehtiet. in Washington Hubert Cremer, wurde zum a. 0. Prof. (daselbst 
Be las National Bureau of Standards, das etwa der ernannt. 

E ee ” a 
; j MURBODEN Ph sikalisch ons ee era Der ord. Prof. der Darstellenden Geometrie an der 
N u ee Pin nn Technischen Hochschule Braunschweig, Dr. Heit.- 
4 Committee of \eronauties | rich Fimerdin 53 wurde Wegen Lrreichens 
3 " der Altersgrenze von seinen amtlichen Verpflich- 
Wie man aus diesem kurzen Bericht sieht. ist tuneen entbunden. 

- DOSE SEEN kun Bramagigen Gastfreund- Der Dozent für Mathematik und für mathema- 

schaft der amerikanischen Fachgenossen, die sich tigeh-technische Grenzgebiete Dr.-Ing. Wilhelm 

natürlich auch in geselligen Veranstaltungen kund- Quade wurde zum a. 0. Professor ernannt. 

tat. zur vollen Befriedieunge aller Teilnehmer ıim % 

Der ord. Professor für anzewandte Mathematik 


Geiste einer echten. sachlich beeründeten inter- 


nstionsien Zusammenarbeit verlaufen. an der Universität Jena, Dr. Max Winkelmann, 


wurde von seinen amtlichen Verpflichtungen ent- 
Dresilen. W. Tollmien. 911 bunden (emeritiert). 
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Einbanddecken für den Jahrgang 1938. 


Erst durch das Einbinden wird der beendete Jahrgang zu einem handlichen und übersichtlichen 
Nachschlagewerk, so daß man seinen Inhalt jederzeit bequem auswerten kann. Wir haben deshalb auch 
für den Jahrgang 1938 der Zeitschrift für angewandte Mathematik und Mechanik wieder Einbanddecken 
herstellen lassen, die zum Preise von 2.25 RM (für VDI-Mitglieder 2.—) durch jede Buchhandlung 


VDI-Verlag G. m. b. H., Berlin NW 7. 


bezoren werden können. 
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Für den Textteil verantwortlich: Professor Dr. Fr. A.Willers, Dresden-A.20, Dorotheenstr. 12. — Printed in Germany. — 
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Wir suchen zum baldigen 
Eintritt: 


Ingenieure 


die sich als Statiker im Flugzeugbau 
betätigen und vorwärtskommen wollen. 


ee Aerodynamiker 

e Physiker 

oe Maihematiker, 
die sich in der statischen Abteilung durch 
Bearbeitung der Lastannahmen eine ent- 
wicklungs- und ausbaufähige Stellung ver- 
schaffen wollen. 
Hauptbedingung: 
Gute mathematisch-technische Veranlagung. 
Günstige Gelegenheit zur Umschulung auf das 
Sondergebiet des Fiugzeugbaues wird vom 
Werk geboten. 
Bewerbungen mit Lebenslauf, Lichtbild, 
Zeugnisabschriften, Angabe von Referenzen 
und Gehaltsansprüchen sowie des frühest 
möglichen Eintrittstages sind zu richten an 


Flieseler-Flugzeugbau 
G.m.b.H. 
Kassel-B 





u Die „Zeitschrift für angewandte Mathematik und Meohanik‘‘ 
A. Bezugsbedingungen = erscheint jeden zweiten Monat (Februar, April, Juni usw.) 
Die Bezugszeit ist halbjährig (Ausland ganzjährig). 


Der Bezugspreis beträgt 1/sjährlich 15.— RM.; ganzjährig 30.— RM. 
für Mitglieder ” 13.50 RH.; m) 27.— BM. 
Bestellung kann bei jeder Buchhandlung oder beim Verlag aufgegeben werden. 


Abbestellung ist nur am Schluß eines Halbjahres bezw. Jahres zulässig. Nicht befristete Bestellungen laufen 
von Halbjahr zu Halbjahr weiter (Ausland von Jahr zu Jahr). 


Bezugsgebühren sind stets zu Beginn eines neuen Halbjahres oder Jahres fällig und werden bei Ausbleiben der 
Zahlung durch Nachnahme eingezogen. 


Lieferung unter Streifband. | 
Einzelheftpreis: 6.— RM., für Mitglieder 5.40 RM, zuzüglich 15 Rpf. Postgeld. 


VDI-Literaturkartei: Auf Wunsch erhalten unsere Bezieher unberechnet Karteikarten über den Inhalt sämt- 
licher Zeitschriften des VDI in monatlichen Zusammensiöllungen ‚gegen Erstattung der jährlichen Versand- 
kosten von 1.— RM. 


ZURBRBEACHTUNG& 


eg rer werden am schnellsten beseitigt, wenn das Ausbleiben der Zeitschrift umgehend mit- | 
geteilt wird. 


Um- und Abbestellungen sind zweckmäßig stets dorthin zu richten, wo die Bestellung aufgegeben wurde, 


Bei Zuschriften und Zahlungen wird um genaue Angabe der vollständigen Anschrift und des Verwendungs- 
zweckes gebeten, da nur dadurch Verwechslungen und Lieferstörungen vermieden werden. 


Satzspiegel der i/ı Seite 171/350 mm. 
B. Anzeigenpreise m site anteiliz. 
Nachlaß: bei 3maliger Aufnahme im Jahr 3 vH, bei 6maliger Aufnahme im Jahr 5 vH. 
Aufschläge für Vorzugsplätze nach vorheriger Vereinbarung laut Tarif. Erfüllungsort für beide Teile Berlin-Mitte. 


_VDI-VERLAG G. M. B. H., BERLIN NW 7, DOROTHEENSTR. 40 


Fernsprecher: Sammel-Nr. 116171: Postscheckkonto: Berlin 102373, Wien: 174439, Prag 77329, Budapest 59951, 

Warschau 194372. Bankverbindungen: Dresdner Bank, Depositen-Kasse 65, Berlin C 2, Wallstr. 5 bis 8. — Deutsche 

Bank, Stadtzentrale Abt. A, Berlin W 8, Mauerstraße 26/27. Drahtanschrift: Ingenieurverlag.. Geschäftszeit: 
Montag bis Freitag 8 bis 17 Uhr, Sonnabend 8 bis 1230 Uhr. 
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